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SUMMARY 


PROJECT  OBJECTIVES,  METHODOLOGY,  AND  OUTLINE  OF  THE 

REPORT 


OBJECTIVES  AND  METHODOLOGY 

This  technical  report  serves  as  the  deliverable  for  a  project  that  involves  (1)  a 
review  of  recent  studies  of  the  mechanical  behavior  of  polycrystalline  silicon  used  to 
fabricate  microelectromechanical  systems  (MEMS)  devices,  (2)  the  development  of 
original  theoretical  models  and  experimental  techniques  that  will  lead  to  improved 
understanding  of  such  behavior,  and  (3)  suggestions  for  future  research.  It  is  hoped  that 
this  report  will  serve  as  a  document  that  will  introduce  those  interested  in  the  subject  to 
recent  developments  in  the  field. 

The  scientist’s  original  proposal  defined  the  following  objectives: 

1 .  Provide  an  assessment  of  the  state  of  the  art  of  theoretical  and  experimental 
mechanics  in  the  field  of  MEMS  technology. 

2.  Develop  analytical  models  to  treat  outstanding  issues,  such  as,  the  appropriate 
scale  for  fracture  mechanics  analyses. 

3.  Develop  fracture  mechanics  specimens  with  characteristic  dimensions  of  the 
same  order  of  magnitude  as  those  of  MEMS  devices. 

4.  Define  the  research  needs  to  predict  the  life  and  strength  of  MEMS  devices. 

5.  Report  on  any  novel  MEMS  applications  that  are  important  to  the  Air  Force, 
•such  as,  sensors  for  the  detection  of  small  cracks  under  fatigue  or  regions 
degraded  by  corrosive  agents. 

The  methodology  was  comprised  of  the  following  tasks: 

1.  Site  visits  to  university  and  industrial  laboratories  involved  in  MEMS 
engineering  and  science. 

2.  Literature  review  and  critique  of  theoretical  and  experimental  studies  of  the 
mechanical  behavior  of  materials  used  to  fabricate  MEMS  devices. 

3.  Search  of  evidence  of  failures  due  to  residual  stresses,  fatigue,  delamination, 
etc.. 

4.  Formulation  of  analytical  models. 

5.  Development  of  fracture  mechanics  specimens. 

6.  Identification  of  research  needs  based  on  Tasks  1-5. 

Task  1  was  accomplished  through  visits  to  the  University  of  Michigan  Center  for 
Integrated  Sensors  and  Circuits,  the  University  of  California  at  Berkeley  Sensor  and 
Actuator  Center,  Microcosm,  and  Analog  Devices.  As  a  result  of  the  discussions  he  had 
with  researchers  at  these  organizations,  the  scientist  quickly  learned  the  state-of-the-art  of 
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MEMS  technology,  especially  the  issues  related  to  the  objectives  of  the  project.  Trip 
(oral)  reports  were  provided  to  Dr.  N.  Pagano  of  Wright  Patterson  Air  Force  Base,  and 
are  not  included  in  this  report. 

Task  2  was  realized  as  follows.  On  the  experimental  side,  the  scientist  decided 
early  on  to  concentrate  solely  on  studies  of  elastic  properties,  tensile  strength,  and 
fracture  toughness  of  polycrystalline  silicon,  the  most  commonly  used  structural  material 
for  MEMS  devices.  Furthermore,  the  scientist  decided  to  limit  this  review  to  work 
published  after  1995.  While  it  is  clear  that  additional  topics  are  important  in  MEMS 
(stiction,  wear,  etc.),  the  properties  treated  in  this  report  are  the  most  fundamental  for 
engineering  design.  Concomitant  to  the  literature  review,  the  scientist  and  his  coworkers 
continued  the  development  of  an  electrostatically  loaded  fracture  mechanics  specimen 
whose  characteristic  lengths  are  comparable  to  those  of  typical  MEMS  devices.  The 
motivation  for  the  electrostatic  loading,  which  eliminates  the  need  of  external  loading 
devices  (as  are  used  in  all  other  test  methods  for  polysilicon  described  in  this  report),  is 
the  development  of  on-chip  test  structures  that  will  allow  in-situ  measurement  of 
strength/fracture/fatigue  of  actual  devices.  This  on-chip  philosophy  was  well  received  by 
the  MEMS  device  developers  that  the  scientist  visited  as  part  of  Task  1. 

It  is  important  to  note  that  the  literature  review  in  this  report  is  not  exhaustive. 
Rather  it  is  comprised  of  a  representative  cross-section  of  the  techniques  available  to 
measure  the  mechanical  properties  of  polycrystalline  silicon.  Furthermore,  most  of  the 
research  groups  that  were  chosen  for  review  consider  mechanical  testing  and  theoretical 
modeling  their  “bread  and  butter”,  not  simply  a  necessary  burden  that  is  part  of  the 
development  and  design  of  specific  MEMS  devices. 

On  the  theoretical  side,  two  issues  have  been  addressed  through  literature  review 
and  model  development.  The  first  is  the  prediction  of  the  scatter  in  the  nominal  elastic 
constants  of  polycrystalline  films.  This  information  is  needed  in  applications  that  rely  on 
the  knowledge  of  stiffness  and  its  uncertainty.  The  two  models  reviewed  are  stochastic, 
the  first  one  developed  by  D.  Mirfendereski  and  coworkers  at  University  of  California  at 
Berkeley,  the  other  by  the  scientist  and  coworkers  at  Case  Western  Reserve  University. 

The  second  theoretical  issue  involves  quantification  of  the  effects  of  heterogeneity 
and  anisotropy  on  crack  tip  parameters  in  cracked  poly  crystalline  films.  This  information 
provides  insight  into  the  appropriate  scale  for  fracture  mechanics  analyses  of  MEMS 
devices.  The  two  models  reviewed  estimate  the  stress  intensity  factors  and  energy  release 
rate  of  a  relatively  long  crack  in  a  polycrystalline  film,  and  provide  useful  information 
about  cracks  in  components  that  contain  a  relatively  large  number  of  anisotropic  grains. 
The  first,  developed  by  K.  Abdel-Tawab  and  coworkers  at  University  of  Texas  at  Austin 
in  their  study  of  fracture  of  columnar  sea-ice,  relies  on  an  approximate  model  for  the 
interaction  of  the  crystals  with  each  other  and  with  the  crack  tip.  The  second,  developed 
by  the  scientist  and  co workers  during  the  course  of  this  project,  relies  on  a  more  accurate 
finite  element  method  based  Monte  Carlo  model.  Both  models  predict  that  the  crack 
driving  forces  are  sensitive  only  to  the  details  of  the  grains  in  the  immediate  vicinity  of 
the  crack  tip.  The  study  of  cracks  whose  lengths  are  comparable  to  the  grain  size  is  one  of 
the  outstanding  problems  that  will  be  discussed  in  Part  4  of  this  report,  as  part  of  Task  4. 
The  information  is  essential  if  one  is  to  perform  fracture  mechanics  analyses  of  very  (in¬ 
plane)  thin  components  of  MEMS  devices,  such  as  connecting  beams  in  resonators. 
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Evidence  of  field  failures  due  to  residual  stress,  fatigue,  delamination  cracking, 
etc.  was  sought,  but  very  little  was  found.  This  is  not  surprising;  the  technology  is  new, 
the  number  of  applications  of  MEMS  is  still  relatively  small,  mechanical  reliability  data 
is  not  available,  and  the  scientist  suspects  there  is  inertia  to  admit  the  possibility  or 
experience  of  failures.  However,  one  laboratory  failure  under  service  conditions  was 
found  and  is  reported,  and  suggests  that  long  term  reliability  of  MEMS  devices 
necessitates  a  basic  understanding  of  the  fracture  mechanics  of  polycrystalline  silicon. 
Additional  fracture  mechanics  phenomena,  including  interfacial  cracking,  will  eventually 
have  to  be  addressed,  since  MEMS  structures  comprised  of  layered  films  are  being 
developed.  An  example  described  briefly  in  Part  4  is  SiC  coated  polycrystalline  silicon 
developed  by  Mehregany  and  coworkers  at  Case  Western  Reserve  University  (K.N. 
Vinod,  C.A.  Zorman,  A.  Yasseen  and  M.  Mehregany,  Fabrication  of  low  defect 
density  3C-SiC  on  Si02  structures  using  wafer  bonding  techniques,  ”  Journal  of 
Electronic  Materials,  in  press). 

MEMS  applications  important  to  the  Air  Force  are  not  described  in  this  report. 

OUTLINE 

The  report  is  divided  into  four  parts,  each  consisting  of  chapters.  Part  1  (Chapters 

1  and  2)  serves  as  an  introduction.  Chapter  1  summarizes  data  for  elastic  modulus, 
strength,  and  fracture  toughness  of  polycrystalline  silicon  films.  The  variability  of  the 
data  clearly  illustrates  the  need  for  robust  measurement  techniques  and  data  bases  for 
mechanical  properties  as  functions  of  processing  procedures,  specimen  size,  etc..  Chapter 

2  presents  a  device  that  failed  in  the  laboratory  as  a  result  of  static  fatigue  (or  stress 
corrosion)  under  the  action  of  applied  loads  and  a  humid  environment.  This  example 
suggests  the  need  to  study  the  fracture  mechanics  of  polycrystalline  silicon,  including 
fracture  toughness  and  subcritical  crack  growth. 

The  experimental  work  performed  by  six  research  groups  is  summarized  in  Part  2, 
which  contains  Chapters  3-8.  In  this  part  of  the  report  the  scientist  has  essentially 
rewritten  the  papers  published  by  the  following  research  groups;  T.  Tsuchiya  and 
coworkers  at  Toyota  Central  Research  and  Development  Laboratories,  W.N.  Sharpe  and 
coworkers  at  Johns  Hopkins  University,  G.  Johnson  and  coworkers  at  the  Berkeley 
Sensor  and  Actuator  Center,  S.  Brown  and  coworkers  at  Exponent  (formerly  Failure 
Analysis  Associates),  the  scientist  and  coworkers  at  Case  Western  Reserve  University, 
and  M.T.A.  Saif  and  coworkers  at  University  of  Illinois  at  Urbana  and  Cornell 
University.  The  work  of  this  last  group  is  associated  with  single  crystal  silicon.  It  is 
summarized  in  this  report  because  it  presents  an  on-chip  strength  test  that  does  not  rely  on 
the  assumption  of  an  elastic  modulus.  The  scientist  chose  this  method  of  writing  Part  2  in 
an  effort  to  eliminate  the  interested  readers’  need  to  search  for  the  original  papers  (in  the 
MEMS  community  many  papers  appear  in  conference  proceedings,  which  as  a  general 
rule  are  more  difficult  to  find  in  libraries). 

The  theoretical  models  are  presented  in  Part  3,  which  is  comprised  of  Chapters  9 
and  10.  These  chapters  summarize,  respectively,  the  aforementioned  predictive  models 
for  elastic  constants  and  crack  driving  forces. 

Part  4  provides  conclusions  and  suggestions  for  future  research. 
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In  this  report  references  are  provided  within  the  text  in  bold  italic  fonts.  This 
nonstandard  format  eliminates  the  need  to  turn  to  a  reference  section  while  reading  the 
text.  The  following  abbreviations  are  used  throughout  the  text. 


B -doping 

BSAC 

COD 

HF 

JHU 

LPCVD 

LTO 

MCNC 

MEMS 

MUMPs 

NSG 

P-doping 

poly-Si 

PSG 

RIE 

SEM 

SNL 

TEM 

TCRDL 

UCB 

UM 

UI 

WPAFB 


boron  doping 

Berkeley  Sensor  and  Actuator  Center 

crack  opening  displacement 

hydrofluoric  acid 

Johns  Hopkins  University 

low  pressure  chemical  vapor  deposition 

low  temperature  oxidation 

Microelectronics  Center  of  North  Carolina 

microelectromechanical  systems 

Multi-User  MEMS  Process 

non-doped  silicate  glass 

phosphorous  doping 

polycrystalline  silicon 

phosphosilicate  glass 

reactive  ion  etching 

scanning  electron  microscope 

Sandia  National  Laboratories 

transmission  electron  microscope 

Toyota  Central  Research  and  Development  Laboratories 

University  of  California  at  Berkeley 

University  of  Michigan 

University  of  Illinois 

Wright  Patterson  Air  Force  Base 
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CHAPTER  1 


VARIATIONS  IN  MECHANICAL  PROPERTIES  OF  POLYSILICON 

Until  recently  the  MEMS  community  has  focused  on  the  development  of 
fabrication  processes  and  clever  devices.  Based  on  his  discussions  with  the  individuals 
acknowledged  in  the  previous  section,  the  scientist  concluded  that  in  order  for  MEMS 
technology  to  reach  its  full  potential,  durability  and  reliability  issues  of  designed  products 
need  to  be  addressed.  Therefore  basic  understanding  of  the  mechanical  behavior  of  the 
materials  used  to  make  devices  is  paramount,  especially  for  devices  that  are  expected  to 
survive  overloads,  a  large  number  of  cycles  of  load,  or  harsh  environments.  As  mentioned 
previously,  in  this  report  attention  is  limited  to  poly-Si.  The  most  basic  static  structural 
design  of  a  poly-Si  device  requires  the  knowledge  of  Young’s  modulus  (and  perhaps 
Poisson’s  ratio)  and  tensile  strength.  The  endurance  limit,  as  obtained  from  an  S-N  curve, 
is  needed  to  ensure  that  the  structure  will  not  fracture  under  cyclic  loads.  Flaw-tolerant 
designs  involve  fracture  toughness,  and  long  term  reliability  necessitates  quantification  of 
subcritical  crack  growth  rates  (as  discussed  in  Chapter  2). 

While  efforts  at  measuring  Young’s  modulus  and  tensile  strength  of  poly-Si  for 
MEMS  devices  have  been  made  over  the  last  decade,  the  data  base  is  sparse,  and 
properties  vary  significantly  as  a  result  of  differences  in  specimen  configuration  and  size, 
property  measurement  techniques,  and  device  fabrication  procedures.  Fracture  toughness, 
stress  corrosion  behavior,  and  fatigue  life  data  is  even  more  scarce.  The  variations  in 
stiffness  and  strength  are  clearly  illustrated  through  Table  1,  which  presents  in  (less 
detail)  the  data  reported  by  W.N.  Sharpe  (W.N.  Sharpe,  Jr.,  B.  Yuan  and  R.L.  Edwards, 
“ Variations  in  mechanical  properties  of  polysilicon”,  Proceedings  of  43rd 
International  Symposium  of  the  Instrumentation  Society  of  America,  Orlando,  FL,  pp. 
179-188,  1997;  W.N  Sharpe,  Jr.,  B.  Yuan,  R.L.  Edwards  and  R.  Vaidyanathan, 
“Measurements  of  Young’s  modulus,  Poisson’s  ratio,  and  Tensile  Strength  of 
Polysilicon”,  Proceedings  of  the  Tenth  IEEE  International  Workshop  on 
Microelectromechanical  Systems,  Nagoya,  Japan  pp.  424-429,  1997)  for  Young’s 
modulus  and  tensile  strength  of  polycrystalline  silicon,  measured  since  1989.  The  last  two 
entries  in  this  table  were  made  by  the  scientist. 
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Table  1-  Variations  in  Young’s  modulus  and  tensile  strength  of  poly-Si. 


Reference 


Tabata  [SI] 


Ta 


Jones [SI 3] 


Modulus  Strength 

(GPa)  (GPa) 


Walker  [S3] 

190 

- 

Koskinen 

[S4] 

164-176 

2.69-3.37 

Maier- 

Schneider 

[S5] 

151-162 

Michalicek 

[S6] 

163 

Biebl  [S7] 

2.11-2.84 

Biebl  [S8] 

130-168 

Biebl  [S9] 

170 

Read  [S10] 

140 

0.7 

Kahn  [SI  1] 

150 

Gupta  [SI 2] 

142-147 

Tsuchiya 
[SI  4] 

167 

2.0-2.7 

Sharpe  [SI 5] 

163-173  . 

1.14-1.27 

Marxer 
[SI  6] 

3.89 

Kahn  [SI 7] 

4.0-4. 9 

The  above  referenced  paper  provides  a  much  more  detailed  description  of  (1)  the 
various  types  and  number  of  tests  used  to  derive  the  data,  (2)  material  processing 
procedures  (type  of  doping,  annealing  temperature,  etc.),  (3)  standard  deviations  of  the 
strength  and  stiffness,  and  (4)  measurement  error  estimates.  As  presented,  this  table 
simply  illustrates  the  relatively  large  variations  in  the  most  fundamental  mechanical 
properties;  Young’s  modulus  varies  from  120-190  GPa,  and  the  tensile  strength  varies 
from  0.7-4. 9  GPa  (the  smaller  values  correspond  to  direct  tension,  while  the  large  values 
are  characteristic  of  flexure).  These  variations  in  turn  provide  a  powerful  argument  for 
studying  the  mechanical  behavior  of  poly-Si. 

Indeed,  Sharpe  suggested  and  organized  a  round-robin  program  for  measuring  the 
stiffness  and  strength  of  poly-Si.  This  involved  having  (1)  MCNC  fabricate  (during  the 
same  fabrication  run)  the  test  structures  used  by  four  different  research  groups,  (2)  each 
group  make  their  measurements  of  Young’s  modulus  and  tensile  strength,  and  (3)  the 
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results  presented  to  the  MEMS  community.  The  results  of  this  very  successful  exercise 
are  presented  in  Part  4. 

The  references  quoted  in  the  Table  above  are  listed  next  for  the  reader  that  may 
want  to  read  the  original  papers. 

[51]  O.  K.  Tabata,  Kawahata,  S.  Sugiyama  and  I.  Igarashi,  “ Mechanical 
property  measurements  of  thin  films  using  load-deflection  of  composite  rectangular 
membranes,”  Sensor  and  Actuators,  20,  pp.  135-141, 1989. 

[52]  Y.-C.  Tai  and  R.S.  Muller,  “ Measurement  of  Young's  modulus  on 
microfabricated  structures  using  a  surface  profiler,  ”  IEEE  Micro  Electro  Mechanical 
Systems,  pp.  147-152, 1990. 

[53]  J.A.  Walker,  K.J.  Gabriel  and  M.  Mehregany,  “ Mechanical  Integrity  of 
Polysilicon  Films  Exposed  to  Hydrofluoric  Acid  Solutions,  ”  Journal  of  Electronic 
Materials,  20,  pp.  665-670, 1991. 

[54]  J.  Koskinen,  E.  Steinwall,  R.  Soave  and  H.H.  Johnson,  “Microtensile 
testing  of  free-standing  polysilicon  fibers  of  various  grain  sizes,”  Journal  of 
Micromechanics  and  Microengineering,  Vol.  25,  pp.  13-17, 1993. 

[55]  D.  Maier-Schneider,  J.  Maibach,  E.  Obermeir  and  D.  Schneider, 
“ Variations  in  Young's  modulus  and  intrinsic  stress  of  LPCVD-polysilicon  due  to 
high-temperature  annealing,  ”  Journal  of  Micromechanics  and  Microengineering,  VoL 
5,  pp.  121-124, 1995. 

[56]  M.A.  Michalicek,  D.E.  Sene  and  V.M  Bright,  “ Advanced  modeling  of 
micromirror  devices,”  Proceedings  of  the  International  Conference  on  Integrated 
Micro/Nanotechnology  for  Space  Applications,  NASA  &  Aerospace  Corp.  Publication, 
pp.  214-229, 1995. 

[57]  M.  Biebl  and  H.V.  Philipsborn,  “. Fracture  strength  of  doped  and 
undoped  polysilicon,”  Transducers  {95-Eurosensors  IX,  Proceedings  of  the  8th 
International  Conference  on  Solid-State  Sensors  and  Actuators,  and  Eurosensors  IX, 
Stockholm,  Sweden,  June  '95,  pp.  72-75, 1995. 

[58]  M.  Biebl,  G.  Brandi  and  R.T.  Howe,  “Young's  modulus  of  in  situ 
phosporous-doped  polysilicon,  ”  Transducers  ‘95-Eurosensors  IX,  Proceedings  of  the 
8"'  International  Conference  on  Solid-State  Sensors  and  Actuators,  and  Eurosensors 
IX,  Stockholm,  Sweden,  June  '95,  pp.  80-83, 1995. 

[59]  M.  Biebl,  T.  Scheiter,  C.  Hierold,  H.V.  Philipsborn  and  H.  Klose, 
“Micromechanics  compatible  with  an  0.8  mm  CMOS  process,”  Sensors  and  Actuators 
A,  Vol  46-47,  pp. 593-597, 1995. 
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[510]  D.T.  Read  and  J.C.  Marshall,  “ Measurements  of  fracture  strength  and 
Young’s  modulus  of  surface-micromachined  polysilicon,  ”  Microlitography  and 
Metrology  in  Micromachining  II,  SPIE  Vol.  2880,  Austin,  Texas,  pp.  56-63, 1996. 

[511]  H.  Kahn,  S.  Stemmer,  K.  Nandakumar,  A.H.  Heuer,  R.L.  Mullen,  R. 
Ballarini  and  M.A.  Huff,  Mechanical  properties  of  relatively  thick,  surface 
micromachined  polysilicon  films,”  the  Ninth  Annual  International  Workshop  on 
Micro  Electro  Mechanical  Systems,  San  Diego,  California,  USA,  pp.  343-353, 1996. 

[512]  R.K.  Gupta,  P.M.  Osterberg  and  S.D.  Senturia,  “ Material  properties 
measurements  of  micromechanical  polysilicon  beams,”  Microlitography  and 
Metrology  in  Micromachining  II,  SPIE  Vol.  2880,  Austin,  Texas,  pp.  39-45, 1996. 

[513]  P.T.  Jones  and  G.C.  Johnson,  “Micromechanical  structures  for  fracture 
testing  of  brittle  thin  films,”  Micro-Electro-Mechanical  Systems,  ASME  DSC-Vol.  59, 
pp.  325-330, 1996. 

[514]  T.  Tsuchiya,  O.  Obata,  J.  Sakata  and  Y.  Taga,  “Specimen  size  effect  on 
tensile  strength  of  surface  micromachined  polycrystalline  silicon  thin  films,” 
Proceedings  of  the  Tenth  IEEE  International  Workshop  on  Microelectromechanical 
Systems,  Nagoya,  Japan,  pp.  529-534, 1997. 

[515]  W.N.  Sharpe,  Jr.,  B.  Yuan,  R.  Vaidyanathan  and  R.L.  Edwards, 
“ Measurement  of  Young’s  modulus,  Poisson’s  ratio,  and  tensile  strength  of 
polysilicon,”  Proceedings  of  the  Tenth  IEEE  International  Workshop  on 
Microelectromechanical  Systems,  Nagoya,  Japan,  pp.  424-429, 1997. 

[516]  C.  Marxer,  M.-A.  Gretillat,  N.F.  de  Rooij,  R.  Battig,  O.  Anthamatten,  B. 
Valk  and  P.  Vogel,  “Reliability  considerations  for  electrostatic  polysilicon  actuators 
using  as  an  example  the  REMO  component,  Sensors  and  Actuators  A  61,  pp.  449-454, 
1997. 


[SI 7]  H.  Kahn,  R.  Ballarini,  R.L.  Mullen  and  A.H.  Heuer,  “Electrostatically 
actuated  failure  of  microfabricated  polysilicon  fracture  mechanics  specimens,” 
submitted  for  publication,  1998. 
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CHAPTER 2 


LABORATORY  FAILURE  OF  A  MEMS  DEVICE 

This  chapter  provides  evidence  that  poly-Si  is  associated  with  static  fatigue.  This 
time  dependent  failure  mechanism,  defined  as  environmentally  assisted  subcritical  crack 
growth  under  an  applied  stress,  must  therefore  also  be  studied  if  long-term  reliability  of 
MEMS  devices  operating  in  an  aggressive  environment  is  to  be  achieved. 

Figure  1  shows  a  fractured  REMO  (a  light  modulator)  component  developed  in 
Switzerland  by  the  Institute  of  Microtechnology,  University  of  Neuchatel,  and  Ascom  (a 
telecommunications  company)  (C.  Marxer,  M.-A.  Gretillat,  N.F.  de  Rooij,  R.  Battig,  O. 
Anthamatten,  B.  Valk  and  P.  Vogel,  Reliability  considerations  for  electrostatic 
polysilicon  actuators  using  as  an  example  the  REMO  component,  Sensors  and 
Actuators  A  61,  pp.  449-454,  1997;  C.R.  Marxer,  “Silicon  micromechanics  for 
applications  in  fiber  optic  communication,  ”  Doctor  of  Science  Dissertation,  University 
of  Neuchatel,  1997). 


Figure  1-  A  MEMS  device  that  fractured  under  the  combination  of  cyclic 
mechanical  stress,  high  electric  field,  and  humid  air  (Reference  provided  in  text). 


This  device  is  an  optical  transceiver  system  with  a  minimum  of  components, 
achieved  by  integrating  a  photodiode  on  the  modulator  chip.  The  poly-Si  actuator 
includes  a  freestanding  membrane  with  a  nominal  thickness  of  464  nm,  separated  from 
the  underlying  electrode  by  an  air  gap  of  1.6  yum.  In  order  for  the  device  to  perform  its 
function,  the  membrane  must  deform  300  nm  under  the  action  of  electrostatic  forces.  The 
electric  field  associated  with  this  displacement  is  as  high  as  60  V  yum'1. 

The  authors  performed  a  standard  reliability  analysis  of  the  device,  which 
consisted  of  the  following  steps. 

1 .  A  finite  element  analysis  of  the  device  was  performed,  yielding  a  maximum 
stress  equal  to  0.13  GPa. 

2.  The  average  and  standard  deviation  of  the  tensile  strength,  evaluated 
experimentally  using  200  pm  long  and  2.8  yum  wide  beams,  were  3.89  GPa 
and  0.53  GPa,  respectively.  The  data  was  fitted  to  the  Weibull  reliability 
function  R(cr)  =  exp[-(cr/cr0)'"],  where  cr  is  the  applied  stress.  This 
produced  a  Weibull  modulus  m=1.5  and  <r0  =  4.15 GPa . 

3.  Using  these  results,  they  calculated  a  practically  zero  probability  of  failure  for 
the  maximum  stress  of  0.13  GPa. 

The  predictions  of  the  analysis  were  confirmed  experimentally.  In  dry  air,  no 
failure  was  observed.  However,  the  rupture  shown  in  Figure  1  was  achieved  in  one 
week’s  time  through  the  application  of  a  field  of  50  V  yum'1  (corresponding  to  a  stress  of 
100  MPa)  and  80%  humidity.  The  cyclic  mechanical  loads  were  produced  by 
superimposing  a  30  V  peak-to-peak  AC  square  signal  of  500  kHz  on  a  DC  bias.  Two 
failure  mechanisms  are  observed;  rupture  of  the  free-standing  membrane,  and 
delamination  of  the  fixed  poly-Si  electrode. 

This  is  not  a  surprising  result;  Connally  and  Brown  {J.A.  Connally  and  S.B. 
Brown,  “ Slow  crack  growth  in  single-crystal  silicon,”  Science,  Vol.  256,  pp.  1537-1539, 
1992 )  demonstrated  that  slow  crack  growth  can  occur  in  single  crystal  silicon. 

This  work  clearly  illustrates  the  need  for  a  basic  understanding  of  environmentally 
assisted  subcritical  crack  growth  in  poly-Si  MEMS  devices.  As  will  be  shown  in  Chapter 
7,  preliminary  data  on  tension-compression  fatigue  loading  of  poly-Si  in  dry  air  suggests 
that  subcritical  crack  growth  can  also  occur  with  very  little  (and  perhaps  no)  humidity. 


10 


PARTI 


REVIEW  OF  EXPERIMENTAL  WORK 


CHAPTER  3 

CHAPTER  4 
CHAPTER  5 
CHAPTER  6 
CHAPTER  7 
CHAPTER  8 


TOYOTA  CENTRAL  RESEARCH  AND  DEVELOPMENT 
LABORATORIES 

JOHNS  HOPKINS  UNIVERSITY 

BERKELEY  SENSOR  AND  ACTUATOR  CENTER 

EXPONENT 

CASE  WESTERN  RESERVE  UNIVERSITY 

CORNELL  UNIVERSITY  AND  UNIVERSITY  OF  ILLINOIS  AT 
URBANA 


11 


CHAPTER  3 


TOYOTA  CENTRAL  RESEARCH  AND  DEVELOPMENT 
LABORATORIES 


1.  INTRODUCTION 

This  chapter,  which  summarizes  the  work  conducted  at  Toyota  Central  Research  and 
Development  Laboratories  (Nagakute  Aichi  480-11  Japan)  by  Toshiyuki  Tsuchiya  and  co¬ 
workers,  is  based  on  the  following  three  references. 

[1]  T.  Tsuchiya ,  O.  Tabata,  J.  Sakata,  and  Y.  Taga,  “ Tensile  testing  of 
poly  crystalline  silicon  thin  films  using  electrostatic  force  grip,”  Transactions  of 
Institute  of  Electrical  Engineers  of  Japan,  Vol.  116-E,  No.10,  pp.  441-446  (1996). 

[2]  T.  Tsuchiya,  O.  Tabata,  J.  Sakata,  and  Y.  Taga,  “ Specimen  size  effect  on 
tensile  strength  of  surface  micromachined  polycrystalline  silicon  thin  films,” 
Journal  of  Microelectromechanical  Systems,  Vol.7,  No.l,pp.  106-113  (1998). 

[3]  T.  Tsuchiya,  J.  Sakata,  and  Y.  Taga,  “ Tensile  strength  and  fracture 
toughness  of  surface  micromachined  polycrystalline  silicon  thin  films  prepared 
under  various  conditions,”  Materials  Research  Society  Proceedings,  Vol.  505  (in 
press). 


This  group  developed  a  testing  procedure  for  measuring  the  tensile  strength  of  poly- 
Si  thin  films  through  direct  tension.  They  subsequently  introduced  a  smooth  edge  notch  in 
their  specimen  to  measure  the  fracture  toughness  (actually  the  bending  strength  in  the  vicinity 
of  a  stress  concentrator,  since  their  notch  is  of  finite  radius).  As  mentioned  by  the  authors,  the 
dimensions  of  the  specimen  are  comparable  to  those  of  connecting  beams  in  typical  MEMS 
devices  (surface  micromachined  resonators,  for  example);  thus  the  experimentally  derived  data 
is  directly  applicable  to  design  and  reliability  of  MEMS  devices. 

Reference  [1]  describes  the  experimental  procedure,  as  well  as  preliminary  results  on 
the  effects  of  specimen  length  on  the  tensile  strength  of  undoped  films.  The  results  of  a  more 
detailed  and  complete  study  of  size  effects  and  doping  on  the  tensile  strength  are  presented  in 
Reference  [2];  the  experimental  data  is  presented  through  Weibull  plots  and  Weibull 
parameters,  with  separate  normalizations  by  volume,  surface  area,  and  side  surface  area.  The 
effects  of  grain  size  (which  depends  on  annealing  temperatures)  on  the  tensile  strength  and 
fracture  toughness  are  studied  in  Reference  [3], 

2.  GENERAL  DESCRIPTION  OF  EXPERIMENTAL  TECHNIQUE 

Mechanical  testing  of  thin  films  poses  many  challenges,  including  how  to  grip  the 
specimen  to  a  loading  device.  This  groups’s  ingenious  design  takes  advantage  of  relatively 
weak  electrostatic  forces  to  grip  a  thin  film  (2pm  thick)  specimen  to  a  probe,  thus  eliminating 
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the  need  for  mechanical  grips  or  adhesives.  Moreover,  because  the  silicon  wafer  that  contains 
the  specimen  can  be  placed  into  and  removed  from  the  testing  chamber  without  being 
touched,  there  is  very  little  danger  of  damaging  the  specimen  before  testing.  The  work  is  truly 
seminal,  as  evidenced  by  the  fact  that  at  least  one  research  group  (JHU,  whose  work  is 
summarized  in  the  next  chapter)  has  substituted  their  adhesive  based  gripping  system  with  the 
electrostatic  system. 

Figure  1  is  a  schematic  of  the  test  set-up.  One  end  of  the  specimen  is  fixed  to  a 
(silicon  wafer)  substrate.  The  other  end,  which  is  free-standing,  is  fixed  to  a  probe  by  an 
electrostatic  force.  The  probe  is  a  conductive  silicon  wafer  covered  with  an  insulating  film  of 
silicon  nitride  (Si3N4).  Figure  2  shows  how  the  tensile  loading  is  applied  to  the  specimen.  A 
sufficiently  high  voltage  is  applied  between  the  specimen  (which  has  the  same  electrical 
potential  as  the  substrate)  and  the  probe.  The  charge  on  their  surfaces  generates  an 
electrostatic  force  that  fixes  them  to  each  other.  As  this  voltage  is  maintained,  the  tensile  load 
is  introduced  by  moving  the  probe  with  a  piezoelectric  actuator.  One  nice  feature  of  this 
gripping  device  is  the  stable  release  procedure  subsequent  to  fracture.  By  applying  the 
opposite  pole  voltage  between  the  probe  and  the  substrate,  the  charge  remaining  on  the 
specimen  generates  a  repulsive  force  that  releases  it  from  the  probe,  and  an  attractive  force 
that  fixes  it  to  the  substrate.  The  unattractive  feature  of  the  loading  device  is  that  the  applied 
force,  as  seen  in  Figure  2,  is  eccentric  with  respect  to  the  neutral  axis  of  the  specimen. 
Therefore  the  specimen  develops  bending  stresses.  These  have  not  been  calculated  nor 
considered  in  the  data  reduction. 
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Figure  1-  Schematic  drawing  of  tensile  testing  using  an  electrostatic  force 
gripping  system  (Reference  [2]). 

Figure  3  shows  the  geometry  and  nominal  dimensions  of  the  uniform  tension  and 
fracture  toughness  specimens.  It  is  important  to  note  that  MEMS  processing  procedures  can 
result  in  significant  differences  between  design  and  as  fabricated  dimensions.  In  the  current 
design,  the  design  values  of  gage  length  L  and  width  w  vary,  respectively,  from  30-300/um  and 
2-5 fim.  As  discussed  subsequently,  this  allows  the  strength  and  toughness  to  be  measured  as 
functions  of  specimen  size.  The  notch  depth  c,  which  is  targeted  for  1/zm  or  2/im,  was 
measured  as  0.81yum  or  1.87/zm.  The  notch  tip  radius  was  reported  in  Reference  [3]  as  0.23 
jim.  A  private  communication  with  T.  Tsuchiya  puts  this  value  at  0.6  fim,  and  finite  element 
derived  stress  concentration  factors  of  5.8  for  vv= 2,  c=l,  a= 2,  and  8.0  for  w=5,  c=2,  a= 4. 
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Figure  2-  Tensile  test  procedure  using  electrostatic  force  (Reference  [2]). 
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Figure  3-  Thin  film  specimen  for  (a)  uniform  tension  (Reference  [2]),  (b)  fracture 

toughness  (Reference  [3]). 


3.  CONCEPTUAL  DESIGN  OF  GRIPPING  SYSTEM 


To  better  understand  the  types  and  magnitudes  of  forces  that  make  this  experiment 
work  (and  the  one  described  in  Chapter  4),  it  is  instructive  to  review  the  calculations  involved 
in  its  conceptual  design.  The  gripping  force  Fv  is  calculated  as  the  electrostatic  force  acting 
perpendicularly  to  two  parallel  electrodes 


F  =  V 2 

v  2d 2 


0) 


where  S  is  the  area  of  the  free  end  of  the  specimen,  d  is  the  gap  between  the  probe  and  the 
free  end,  V  is  the  applied  voltage,  and  £0  =  8.854xl0-12  F  /  m  is  the  permittivity  constant  (of 
air  in  this  case).  In  order  for  the  specimen  to  stick  to  the  probe,  its  weight  should  be  less  than 
Fv.  This  provides 


V 2  -mg  =  “  StPS  >  0  (2) 

where  g  is  the  acceleration  due  to  gravity,  and  m,  t  and  p  are  the  mass,  thickness,  and  density 
of  the  specimen,  respectively.  The  applied  voltage  associated  with  electrical  breakdown  of  the 
insulating  film  (several  MV/cm)  is  100V.  For  this  voltage  a  several  micron  thick  polysilicon 
film  (p  =  2.33xl03 kg  /  m3 )  will  stick  at  a  gap  <i=lmm. 

In  order  for  the  specimen  not  to  slip  relative  to  the  probe,  a  force  parallel  to  the 
tensile  force  is  required.  This  force  is  comprised  of  a  horizontal  electrostatic  force  Fh  and  a 
friction  force  Ff  ,  given  by 


£  £  W 

0  r  8  y2 

2d ' 

(3) 

II 

(4) 

where  ws  is  the  width  of  the  free  end,  and  d]  and  £r  are,  respectively,  the  gap  and  dielectric 
constant  between  electrodes.  When  the  specimen  is  fixed,  d  is  equal  to  the  thickness  of  the 
Si3N4  insulating  film  (d= 0.2  pm  and  £r  =8).  The  resultant  of  these  two  forces  must  be  greater 
than  the  fracture  tensile  force  F, ,  which  is  given  in  terms  of  the  tensile  strength  of  ,  thickness 
t,  and  width  w,  as  <7fwt .  The  following  sensitivity  analysis  shows  that  the  friction  force  is 

much  larger  than  the  horizontal  electrostatic  force.  Assume  that  the  coefficient  of  friction 
p  =  0.3.  Since  the  tensile  strength  is  on  the  order  of  2-5GPa,  the  largest  value  of  F,  is  about 
0.1N.  Assuming,  for  the  sake  of  argument,  that  the  free  end  of  the  specimen  is  square, 
S  -  w2,  for  100V  the  necessary  dimensions  of  the  electrode  are  200pm  x200pm.  As  a  result. 
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Ff  is  three  orders  of  magnitude  larger  than  Fh. 


4.  SPECIMEN  FABRICATION 

The  specimens  described  in  References  [l]-[3]  were  fabricated  using  standard 
micromachining  technology.  They  are  poly-Si  films  crystallized  from  an  amorphous  film.  The 
authors  favor  crystallized  polycrystalline  silicon  films  because  they  are  associated  with  a 
tensile  residual  stress,  making  them  suitable  for  micromachined  devices,  and  because  their 
grain  size  can  be  controlled  by  the  annealing  temperature.  Crystallization  occurs  above  600°C, 
the  grain  size  decreasing  with  increasing  annealing  temperature.  In  Reference  [1],  the 
annealing  temperature  was  1000°C.  In  Reference  [3],  the  effects  of  annealing  temperature  and 
associated  grain  size  were  investigated  by  testing  films  that  were  amorphous,  or  annealed  at 
600°C,  700°C,  or  1000°C;  in  this  paper  the  effects  of  additional  annealing  at  1000°C  were  also 
investigated.  The  effects  of  P-doping  were  studied  in  Reference  [2]. 

The  specimen  is  fabricated  using  the  surface  micro-machining  process  shown  in 
Figure  4.  The  following  key  steps  were  reported  in  Reference  [1]. 

1.  A  LPCVD  amorphous  silicon  film  is  deposited  at  520°C  on  the  substrate,  which 
is  comprised  of  a  0.2^m  thick  LPCVD  Si3N4  insulating  film  and  a  2pm  plasma 
CVD  NSG  sacrificial  film. 

2.  The  amorphous  film  is  crystallized  by  annealing  in  nitrogen  at  1000°C. 

3.  In  Reference  [2]  P-doping  of  the  poly-Si  film  was  performed  using  POCL3  at 
1000°C.  The  grain  size  of  the  poly-Si  film  is  0.3^m,  irrespective  of  the  doping. 

4.  The  film  is  patterned  into  the  shape  of  the  specimen  by  photolitography  and  RIE. 
An  aluminum  pad  electrode  is  fabricated  on  the  fixed  end  of  the  specimen.  The 
free  end  of  the  specimen  has  many  etching  holes  which  shorten  the  time  needed 
for  sacrificial  etching  time. 

5.  The  sacrificial  layer  is  removed  by  wet  etching  with  HF. 

6.  The  specimen  is  rinsed  with  distilled  water,  dipped  with  hot  p-dichlorobenzene, 
and  then  dried.  This  last  step  prevents  the  specimen  from  sticking  to  the 
substrate. 

Figure  5  shows  a  scanning  electron  micrograph  of  three  typical  specimens.  The 
properties  of  the  non-doped  and  P-doped  poly-Si  films,  as  reported  in  Reference  [2],  are 
shown  in  Table  1.  The  internal  stress  and  the  Young’s  modulus  were  measured  on  a 
rectangular  membrane  structure  using  the  pressure-deflection  technique  (O.  Obata,  K. 
Kawahara,  S.  Sugiyama,  and  I.  Igarashi,  “ Mechanical  property  measurement  of  thin 
films  using  load-deflection  of  composite  rectangular  membranes,”  Sensors  and  Actuators, 
Vol.  20,  pp.  135-141, 1989). 

In  Reference  [3]  films  with  different  average  size  grains  were  fabricated  by  annealing 
at  600°C,  700°C  and  1000°C.  Scanning  electron  microscope  (SEM)  micrographs  of  the 
different  microstructures  produced  by  these  temperatures  are  shown  in  Figure  6;  the  600°C 
and  1000°C  annealed  films  are  comprised  of  spherical  grains  of  approximate  size  2pm  and 
0.3/im,  respectively,  while  the  700°C  annealed  film  contains  oval  grains  of  about  0.1pm  size. 
As  discussed  previously,  some  of  the  films  went  through  an  additional  annealing  at  1000°C, 
which  had  no  effect  on  the  grain  size. 
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Table  1  Properties  of  poly-Si  thin  film  (Reference  [2]). 


non-doped 

P-doped 

Grain  size 

0.3  pm 

0.3pm 

Thickness 

1.8pm 

2.4pm 

Internal  Stress 

76MPa 

-7.0MPa 

Young’s  modulus 

163GPa 

167GPa 

Figure  4-  Fabrication  process  (a)  deposition  of  films,  (b)  patterning  of  poly-Si  and  A1 
films,  (c)  sacrificial  etching  to  release  poly-Si  (Reference  [2]). 


Figure  6-  SEM  micrograph  of  poly-Si  film  crystallized  at  (a)  600°C,  (b)  700°C, 
and  (c)  1000°C  (Reference  [3]).  Note  the  decrease  in  grain  size  with  increasing 

temperature. 


5.  TESTING  PROCEDURE 

As  shown  schematically  in  Figure  7,  the  tension  test  is  conducted  (in  a  vacuum) 
inside  an  SEM  chamber.  This  figure  shows  the  locations  of  the  piezoelectric  actuator,  the 
strain  gauge  that  measures  the  force,  and  the  strain  gauge  that  measures  the  displacement. 
Figure  8  is  an  SEM  micrograph  of  the  probe  gripping  a  specimen.  This  figure  clearly  shows 
the  eccentric  nature  of  the  loading  device,  and  the  resulting  bending  of  the  specimen.  The  ratio 
of  maximum  bending  stress  to  axial  stress  is  6 e  / 1 ,  where  e  is  the  distance  from  the  line  of 
action  of  the  force  to  the  neutral  surface  of  the  specimen  within  the  gauge  section. 


X-Y  Recorder  SEM  Chamber  ^—Attached 

Mechanism 


Figure  7-  SEM  chamber  testing  environment  (Reference  [3]). 
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Figure  8-  SEM  micrograph  of  tensile  testing  (Reference  [2]).  Note  that  the  specimen 
is  subjected  to  eccentric  loading,  as  evidenced  by  its  bent  shape. 

6.  RESULTS 

This  section  is  divided  into  three  parts.  The  effects  of  specimen  size  and  grain  size  on 
the  tensile  strength  are  presented  first.  This  is  followed  by  the  effects  of  grain  size  on  fracture 
toughness.  The  final  section  discusses  fracture  initiation  site. 

Tensile  strength 

Effect  of  specimen  size  and  doping 

Fracture  is  always  brittle,  and  initiates  at  random  positions  within  the  tested  length. 
The  tensile  strength  for  various  gage  lengths,  as  reported  in  Reference  [1],  is  presented  in 
Figure  9.  In  this  figure,  n  represents  the  number  of  specimens  tested.  As  expected  for  this 
brittle  material,  the  average  strength  decreases  with  increasing  gage  length.  This  size  effect  is 
studied  in  a  more  systematic  fashion  in  Reference  [2],  as  discussed  next 

The  mean  strength,  standard  deviation,  and  Weibull  modulus  for  the  un-notched 
specimens,  taken  from  Reference  [2],  are  summarized  in  Table  2.  The  mean  and  the  range  of 
the  tensile  strength  are  shown  in  Figure  10. 

Because  the  authors  assumed  that  the  strength  data  follows  a  Weibull  distribution, 


the  following  results  for  this  extreme  value  distribution  are  stated  for  completeness.  The 
probability  of  failure  under  an  applied  stress  cr  is  given  as 


F  =  1  -  exp 


(5) 


where  m  is  the  Weibull  modulus,  VE  is  the  effective  volume  of  the  tested  part,  and  a  is  a 
constant.  The  mean  tensile  strength  fi  and  variance  Var  are  given  by 
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where  TO)  is  the  Gamma  function.  It  follows  that  the  mean  strengths  &  and  of 
specimens  with  effective  volumes  VE1  and  VE2  are  related  by 
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For  a  fixed  a  and  VE ,  a  plot  of  Equation  (7)  would  show  that  as  m  decreases,  the  scatter 
increases. 

The  fracture  probability  functions  from  Reference  [2]  are  plotted  in  Figures  11. 
These  were  fitted  to  the  Weibull  distribution  function  using  a  least  square  method,  the 
resulting  parameter  m  tabulated  in  Table  2.  It  is  first  noted  that  the  fitted  distributions  do  not 
fit  the  data  very  well  at  the  tail  ends;  this  is  not  surprising  because  many  more  data  points  are 
required  to  accurately  fit  such  distributions. 

The  strengths  of  the  un-doped  specimens  reported  in  Reference  [1]  are  very  close  to 
those  reported  in  Reference  [2];  this  gives  confidence  to  the  reproducibility  of  the  results.  It  is 
observed  that  the  mean  strength  of  the  P-doped  films  is  slightly  lower  than  of  the  un-doped 
films.  However,  the  strength  of  the  un-doped  films  exhibits,  for  the  same  size  specimen,  a 
larger  scatter  (Weibull  modulus  5.3-6.9)  than  that  of  the  P-doped  films  (Weibull  modulus  8.7- 
12.6). 

The  most  important  result  of  this  work  is  perhaps  the  implication  that  if  poly-Si  is  to 
be  used  to  its  full  potential  (in  terms  of  strength),  processing  techniques  must  be  adjusted  to 
mitigate  the  variability  in  strength.  This  is  demonstrated  clearly  by  considering  the  ratio  of 
allowable  stress  crfl//  to  mean  strength  for  a  prescribed  reliability.  Assume  for  the  sake  of 
argument  that  the  data  is  indeed  Weibull.  Combining  Equations  (5)  and  (6)  this  ratio  becomes 
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(9) 


_[-ln(l-f)P 

6r) 

If  reliability  considerations  demand  that  only  one  part  in  a  million  fail  (F  - 10~6),  then  for 
m= 5,  6.5  and  10,  the  allowable  stress  is  only  6.8%,  12.8%  and  26%  of  the  mean  strength.  In 
other  words,  these  levels  of  variability  in  strength  require  a  very  conservative  design  approach. 

As  discussed  in  the  final  section  of  this  chapter,  there  is  evidence  that  in  these  films, 
the  fracture  initiates  at  surface  flaws  formed  by  the  sacrificial  etching  procedure.  Therefore  the 
most  effective  way  of  reducing  scatter  in  strength  is  through  improved  processing  procedures 
that  mitigate  the  development  of  these  stress  concentrators. 


Figure  9-  Tensile  strength  for  different  gage  lengths  (Reference  [1]) 
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Table  2  Tensile  strength  of  poly-Si  thin  film  against  the  specimen  size 


(Reference  [2]). 


non-doped 

P-doped 

Width  (pm) 

2 

5 

2 

5 

Length  (pm) 

30 

100 

300 

30 

100 

300 

30 

100 

300 

30 

100 

300 

Sample  Number 

11 

14 

17 

13 

18 

17 

11 

12 

11 

13 

12 

15 

Mean  Strength  (GPa) 

2.8 

2.3 

2.2 

2.7 

2.4 

2.0 

2.5 

2.2 

2.0 

2.4 

2.1 

2.1 

Std.  Dev.  (GPa) 

0.61 

0.52 

0.45 

0.48 

0.44 

0.45 

0.37 

0.20 

0.20 

0.25 

0.19 

0.26 

Weibull  modulus 

5.3 

6.0 

6.0 

6.9 

6.3 

5.4 

8.7 

12.0 

9.9 

10.1 

12.6 

9.3 

0.0  11  . . |  ,  [i  ^|  |  i  i  i  t-—- 

w  (pm)  2  2  2  5  5  5 

L(pm)  30  100  300  30  100  300 

(a) 
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L  (jj.m)  30  100  300  30  100  300 

(b) 

Figure  10-  Tensile  strength  of  poly-Si  thin  film;  (a)  un-doped  film,  (b)  P-doped  film 

(Reference  [2]). 
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Figure  11-  Measured  fracture  probability  as  a  function  of  tensile  stress.  The  curve  is 
fitted  to  Weibull  distribution  function  (Reference  [2]). 
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Effect  of  grain  size 

In  Reference  [3]  the  effects  on  tensile  strength  of  the  previously  described  different 
annealing  conditions  and  their  associated  grain  size  were  studied.  The  results  are  presented  in 
Figure  12.  It  is  observed  that  for  all  specimen  sizes,  amorphous  films  have  the  lowest  average 
tensile  strength,  and  that  the  strength  increases  with  increasing  annealing  temperature 
(decreasing  grain  size).  The  authors  attributed  the  lower  strength  in  the  amorphous  films  to 
the  presence  of  0.4%  hydrogen.  For  the  specimens  that  were  subjected  to  additional 
annealing,  there  seems  to  be  a  mixed  pattern;  the  additional  annealing  slightly  strengthened  the 
600°C  specimens,  slightly  weakened  the  700°C  specimens,  and  had  no  effect  on  the  1000°C 
specimens. 
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Figure  12-  Tensile  strength  of  poly-Si  films  for  various  annealing  conditions  (and  grain 

size  (Reference  [3])). 


Fracture  toughness 

The  notched  beam  shown  in  Figure  3(b)  was  used  to  measure  the  fracture  toughness 
(Klc).  Assuming  that  the  notch  is  a  sharp  crack  of  equal  depth,  available  handbook  formulas 
were  used,  together  with  the  fracture  load,  to  derive  the  nominal  fracture  toughness  values  in 
Figure  13.  It  is  observed  that  these  values  range  from  1.9  to  4.5  MPa-m1/2. 

Because  the  notch  tip  has  a  radius  comparable  to  other  characteristic  dimensions  (the 
most  important  being  the  lengths  of  the  notch  and  uncracked  ligament),  these  results  do  not 
represent  the  true  toughness  of  the  material;  this  parameter  requires  propagation  from  an 
atomically  sharp  crack.  The  results  can,  at  best,  be  used  to  determine  the  bending  strength  in 
the  region  of  a  stress  concentrator  (5.8  for  w=2,  c=l,  a= 2,  and  8.0  for  vv=5,  c=2,  a= 4),  or  the 
critical  J-integral  for  the  given  notch-tip  radius.  This  method  of  presenting  results  obtained 
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using  a  finite  radius  notch  is  favored  and  recommended  by  the  scientist  and  his  coworkers,  as 
discussed  in  Chapter  7. 

In  Reference  [3]  the  authors  also  back  calculate  an  initial  flaw  size  for  the  tension 
specimens,  and  concluded  that  while  the  critical  flaw  sizes  did  not  correlate  well  with  grain 
size,  they  did  correlate  well  with  the  surface  pits  produced  by  the  etching  procedure.  While 
the  conclusion  may  be  correct,  the  procedure  they  used  is  questionable.  The  critical  stress  is 
proportional  to  the  fracture  toughness  and  the  inverse  of  the  square  root  of  the  critical  flaw 
size  (knowledge  of  two  parameters  provide  the  third).  The  authors  used  the  critical  stress 
from  the  tension  test,  the  fracture  toughness  from  the  notch  tension  specimens,  and 
determined  the  flaw  size.  This  procedure  is  questionable  because  the  notch  tests,  as  discussed 
previously,  do  not  provide  the  true  toughness. 
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Figure  13-  Fracture  toughness  of  poly-Si  films  for  various  annealing  conditions  (and 

grain  size  (Reference  [3])). 


Determination  of  fracture  initiation  site 

The  ideal  method  of  determining  where  fractures  initiate  is  inspection  of  the  fracture 
surfaces.  Figure  14  taken  from  Reference  [2],  shows  a  typical  fracture  surface,  and  indicates 
that  the  fracture  origin  is  on  the  surface  of  the  specimen.  According  to  the  authors,  this  was 
typical  of  all  observed  fracture  surfaces.  Similar  fracture  surfaces  were  observed  by  the 
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scientist  and  coworkers,  as  described  in  Chapter  7. 

The  authors  also  used  an  unconventional  statistical  method  to  determine  whether  the 
fracture  initiated  within  the  volume,  on  the  total  surface,  or  on  the  side  surface.  They  used  the 
maximum  likelihood  method  to  fit  the  data  to  the  Weibull  distribution  function  using  these 
geometric  parameters  in  Equation  (5).  The  results  are  shown  in  Table  3  and  Figure  15.  Table 
3  presents  the  Weibull  moduli  m,  the  likelihood  functions  L,  and  reliability  coefficients  R2 
for  each  of  these  geometric  parameters.  R2 ,  calculated  by  EXCEL,  is  a  normalized  residual  of 
the  sum  of  the  squared  difference  between  the  measured  values  and  the  fitted  values.  The 
functions  with  the  length  and  the  side  surface  area  are  associated  with  (very  slightly  and 
statistically  insignificant)  larger  values  of  the  likelihood  function  and  reliability  coefficients. 
The  authors  concluded  that  these  slight  differences  suggest  surface  flaw  driven  fracture. 


Table  3  Fitting  results  using  the  maximum  likelihood  method 

(Reference  [2]). 


fitting  parameter 

length 

l 

width 

w 

jail  surface  sidesur- 
area  face  area 
i  S  Ss 

volume 

v 

non-doped 

Weibull  modulus  j  m 

6.3 

!  6.5 

6.3 

6.6 

likelihood  function  llogL 

-1927 

I  -1927 

]  -1927 

-1930 

reliability  coefficient  j  R2 

0.818 

|  0.803 

0.818 

|  0.763 

P-doped 

Weibull  modulus  j  m 

10.0 

7.9 

j  10.2 

10.3 

likelihood  function  j  log  L 

-1541 

|  -1544 

-1541  | 

-1548 

reliability  coefficient  i  R2 

0.790 

|  0.775  | 

0.790  I 

|  0.732 
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Fracture  origin 


Figure  14-  Fracture  surface  of  poly-Si  thin  film  (Reference  [2]) 


Figure  15-  Weibull  plots  normalized  by  the  volume,  the  surface  area,  and  side  surface 

area  (Reference  [2]). 


7.  CRITIQUE 

The  electrostatic  loading  system  developed  by  the  research  group  at  TCRDL  is  truly 
path-breaking;  it  offers  an  effective  method  for  loading  thin  films  with  characteristic 
dimensions  comparable  to  those  of  typical  MEMS  devices.  The  group  has  also  started  on  the 
correct  path  as  they  use  their  system  to  (1)  gain  a  better  understanding  of  the  effects  of 
processing  procedures  on  the  mechanical  properties  of  poly-Si,  (2)  reducing  the  data  using 
standard  statistical  analyses,  which  together  with  stress  analysis,  provide  some  measure  of 
device  reliability,  and  (3)  identifying  the  origin  of  fracture  and  methods  of  mitigation. 
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CHAPTER  4 


JOHNS  HOPKINS  UNIVERSITY 


1.  INTRODUCTION 

This  chapter  summarizes  the  development  by  W.N.  Sharpe,  Jr.  and  coworkers,  at 
Johns  Hopkins  University,  of  a  testing  environment  for  measuring  the  Young’s  modulus, 
Poisson’s  ratio,  direct  tensile  strength,  and  fracture  toughness  of  poly-Si  thin  films.  This 
group’s  efforts  have  produced  the  following  publications  (the  list  was  kindly  provided  by 
W.N.  Sharpe,  Jr.). 

[1]  B.  Yuan,  and  W.N.  Sharpe,  Jr.,  Mechanical  Testing  of  Polysilicon 
Thin  Films  with  the  ISDG”,  Experimental  Techniques,  vol.21,  no.2,  p.  32-5, 1997. 

[2]  W.N.  Sharpe,  Jr.,  R.  Vaidyanathan,  B.  Yuan,  G.  Bao  and  R.L.  Edwards, 
“Effect  of  Etch  Holes  on  the  Mechanical  Properties  of  Polysilicon”  Journal  of 
Vacuum  Science  &  Technology,  vol.15,  no.5,  pp.  1599-603, 1997. 

[3]  W.N.  Sharpe,  Jr.,  B.  Yuan  and  R.L.  Edwards,  “ A  New  Technique  for 
Measuring  the  Mechanical  Properties  of  Thin  Films”,  Journal  of 
Microelectromechanical  Systems,  Vol.  6,  No.  3,  pp. 193-199, 1997. 

[4]  D.J.  Chang  and  W.N.  Sharpe,  Jr.,  Mechanical  Analysis  and  Properties 
of  MEMS  Materials,  Microengineering  for  Aerospace  Systems .  The  Aerospace  Press  of 
The  Aerospace  Corporation,  in  press,  1998. 

[5]  W.N.  Sharpe,  Jr.  and  B.  Yuan,  “Mechanical  Testing  of  MEMS 
Materials ”,  Proceedings  of  the  1995  SEM  Spring  Conference  on  Experimental 
Mechanics,  Grand  Rapids,  pp.  698  -  701, 1995. 

[6]  W.N.  Sharpe,  Jr.,  B.  Yuan,  R.  Vaidyanathan,  and  Edwards,  R.  L., 
“Direct  Measurements  of  Young’s  Modulus  and  Tensile  Strength  of  Polysilicon ”, 
Late  News  Proceedings  of  the  Solid-State  Sensor  and  Actuator  Workshop,  Hilton  Head 
Island,  SC,  pp.  27  -  28, 1996. 

[7]  W.N.  Sharpe,  Jr.,  B.  Yuan,  R.  Vaidyanathan  and  R.L.  Edwards,  “New 
Test  Structures  and  Techniques  for  Measurement  of  Mechanical  Properties  of  MEMS 
Materials  ”,  Proceedings  of  the  SPIE  Symposium  on  Microlithography  and  Metrology 
in  Micromachining  II,  Austin,  TX,  pp.  78-91, 1996. 

[8]  W.N.  Sharpe,  Jr.,  B.  Yuan,  R.L.  Edwards  and  R.  Vaidyanathan, 
“ Measurements  of  Young’s  modulus,  Poisson’s  ratio,  and  Tensile  Strength  of 
Polysilicon”,  Proceedings  of  the  Tenth  IEEE  International  Workshop  on 
Microelectromechanical  Systems,  Nagoya,  Japan,  pp.  424-429, 1997. 

[9]  W.N.  Sharpe,  Jr.,  B.  Yuan  and  R.L.  Edwards.,  “ Variations  in 
Mechanical  Properties  of  Polysilicon”,  Proceedings  of  43rd  International  Symposium 
of  the  Instrumentation  Society  of  America,  Orlando,  FL,  pp.  179-188, 1997. 

[10]  W.N.  Sharpe,  Jr.,  D.A.  LaVan  and  R.L.  Edwards.,  “Mechanical 
Properties  of  LIGA-Deposited  Nickel  for  MEMS  Transducers”,  Proceedings 
Transducers  ‘97,  Chicago,  IL,  pp.  607-610, 1997. 

[11]  W.N.  Sharpe,  Jr.  and  R.L.  Edwards,,  “Advances  in  Tensile  Testing  of 
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Polysilicon  Thin  Films”,  Applications  of  Experimental  Mechanics  to  Electronic 
Packaging”,  ASMEEEP-VOL.  22/AMD-VOL.  226,  pp.  113-116, 1997. 

[12]  W.N.  Sharpe,  Jr.,  D.A.  LaVan.  and  A.  McAleavey,  “Mechanical  Testing 
of  Thicker  MEMS  Materials ”,  Micro-Electro-Mechanical  Systems  (MEMS),  AS  ME 
DSC-Vol.  62/HTD-Vol.  354,  pp.  93-97, 1997. 

[13]  W.N.  Sharpe,  Jr.,  R.  Vaidyanathan,  B.  Yuan  and  R.L.  Edwards,  “A  New 
Technique  for  Measuring  Poisson’s  Ratio  of  MEMS  Materials”,  Materials  for 
Mechanical  and  Optical  Microsystems  Symposium,  Materials  Research  Society,  pp. 
185-190, 1998. 

[14]  W.N.  Sharpe,  Jr.  and  B.  Yuan,  “Fracture  Tests  of  Polysilicon  Film”, 
presented  at  Thin  Films  -  Stresses  and  Mechanical  Properties  -  Symposium  NN, 
Materials  Research  Society  and  submitted  to  the  proceedings  (December,  1997). 

[15]  W.N.  Sharpe,  Jr.,  K.  Turner  and  R.L.  Edwards,  “ Polysilicon  Tensile 
Testing  with  Electrostatic  Gripping”,  presented  at  Microelectromechanical  Structures 
for  Materials  Research  -  Symposium  N,  Materials  Research  Society  and  accepted  for 
publication  in  the  proceedings  (April,  1998). 

[16]  W.N.  Sharpe,  Jr.,  K.  Turner  and  R.L.  Edwards,  “Measurements  of  the 
Effect  of  Specimen  Size  on  Young’s  Modulus  and  Tensile  Strength  of  Polysilicon”, 
Late  News  Proceedings  of  the  Solid-State  Sensor  and  Actuator  Workshop,  Hilton  Head 
Island,  SC,  accepted  for  publication  (Junel998). 

What  makes  the  above  referenced  work  unique  is  the  use  of  laser  interferometry  to 
accurately  measure  the  relative  displacements  (of  gold  lines  deposited  on  the  specimen) 
needed  to  calculate  axial  and  lateral  strains.  As  a  result,  the  authors  have  measured  not 
only  the  tensile  strength  and  fracture  toughness  (using  a  finite  radius  center  notch)  of  poly- 
Si,  but  also  its  Young’s  modulus  and  Poisson’s  ratio. 

The  first  generation  testing  procedure,  described  in  References  [1]-[14],  uses  an 
adhesive  to  grip  the  ends  of  the  (3.5  pm  thick,  600  pm  wide,  and  4000  pm  long)  tensile 
specimen  to  a  testing  machine.  Very  reproducible  results  were  obtained,  suggesting  that 
the  experimental  and  material  processing  procedures  are  consistent  and  reproducible. 

Unfortunately,  the  fabrication  is  associated  with  only  one  specimen  per  die,  and 
involves  a  time  consuming  etching  procedure.  Moreover,  the  dimensions  of  the  specimen 
in  the  plane  of  the  film  are  relatively  large  compared  to  typical  components  of  MEMS 
devices,  and  therefore  the  scatter  in  the  experimental  measured  properties  is  not  expected 
to  represent  the  scatter  in  typical  device  components. 

For  these  reasons,  the  group  developed  a  second  generation  testing  system, 
described  in  References  [15]-[16],  which  involves  an  electrostatic  gripping  procedure 
similar  to  the  one  developed  by  Tsuchiya  and  co workers  (described  in  the  previous 
chapter).  This  significant  modification  makes  it  possible  to  produce  many  different  size 
specimens,  each  of  smaller  dimensions,  on  one  die.  This  author’s  opinion  is  that  when  the 
system  is  fine  tuned,  it  will  be  the  most  reliable  system  available  for  in-laboratory 
fundamental  studies  on  mechanical  properties  of  poly-Si. 

Because  the  above  referenced  work  is  associated  with  two  different  procedures, 
this  chapter  is  divided  into  two  parts.  The  first  describes  the  development  of  the  first 
generation  gripping  system  and  the  laser  interferometry  technique,  as  well  as '  the 
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measurements  made  with  this  system  of  Young’s  modulus  and  Poisson’s  ratio,  tensile 
strength,  and  fracture  toughness.  The  second  part  summarizes  the  group’s  adoption  of  the 
electrostatic  gripping  device  and  preliminary  results  of  Young’s  modulus,  Poisson’s  ratio, 
tensile  strength,  and  fracture  toughness. 

2.  FIRST  GENERATION  GRIPPING  SYSTEM 

General  description  of  experimental  technique 

A  schematic  of  the  testing  environment  is  shown  in  Figure  1.  The  tensile  specimen, 
which  is  described  subsequently,  is  placed  between  two  grips,  one  attached  to  a  fixed  base, 
the  other  to  a  linear  air  bearing  that  eliminates  friction  in  the  loading  mechanism.  The  force 
is  applied  through  a  piezoelectric  translator  (with  a  range  of  180/rm  and  resolution  of 
0.09/un,  as  reported  in  Reference  [7]),  and  measured  by  a  4.5N  load  cell  (with  resolution 
of  0.0045N,  as  reported  in  Reference  [7]).  The  test  is  run  under  displacement  control. 


Laser 


Figure  1-  Schematic  of  the  test  setup  (Reference  [8]). 

The  Interferometric  Strain/Displacement  Gage  (ISDG)  is  used  for  measuring 
strains.  This  optical  technique  measures  the  relative  displacement  between  reflective  gage 
markers,  as  shown  schematically  in  Figure  2.  In  this  work  the  markers  are  gold  lines, 
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Figure  2-  Schematic  of  the  ISDG  (Reference  [8]). 

10/zm  wide,  0.5^m  high,  200/mi  long,  and  250pm  apart,  as  shown  in  the  SEM  micrograph 
of  Figure  3.  In  this  figure  there  are  two  sets  of  lines,  one  set  perpendicular  to  the  other. 
One  set  is  used  to  measure  the  relative  displacement  in  the  direction  of  the  applied  force, 
the  other  in  the  direction  perpendicular  to  the  applied  force.  This  allows  measurement  of 
both  Young’s  modulus  and  Poisson’s  ratio.  As  explained  subsequently,  these  lines  are 
deposited  on  the  specimen  during  processing. 

Strain  is  measured  as  follows.  An  incident  laser  beam  is  used  to  illuminate  the  gold 
lines.  The  edges  of  the  lines  produce  diffracted  reflections  that  overlap  and  the  interfere. 
This  produces  fringe  patterns  that  move  as  the  lines  undergo  relative  displacement.  The 
motion  of  the  fringe  patterns  is  measured  with  fringe  detectors,  and  in  turn  the  strain  is 
calculated.  The  authors  report  in  Reference  [8]  that  their  ISDG  system  is  associated  with  a 
relative  uncertainty  of  ±  3%.  The  phenomena  that  contributes  to  this  uncertainty  is  stated 
in  Reference  [8],  The  optical  principles  and  practical  aspects  of  the  ISDG  technique  have 
been  reviewed  by  W.N.  Sharpe,  an  authority  on  these  systems,  in  W.N.  Sharpe ,  Jr.,  “An 
interferometric  strain  displacement  measurement  system,”  NASA  Technical 
Memorandum  101638, 1989. 


Figure  3-  SEM  photo  of  the  biaxial  strain  gage  (Reference  [8]). 
Description  of  specimen 

The  poly-Si  specimen  is  shown  mounted,  using  an  adhesive,  in  the  grips  of  the  test 
machine  in  the  center  of  Figure  4.  There  are  two  support  strips  parallel  to  the  specimen 
that  prevent  damage  of  the  specimen  during  fabrication  and  handling.  These  are  cut  (with 
a  diamond  saw),  as  shown  in  the  figure,  immediately  before  testing.  Details  of  the 
alignment  and  precautions  taken  during  handling,  gluing  and  curing  are  presented  in 
Reference  [7]. 


Figure  4-  A  polysilicon  specimen  mounted  in  the  grips  of  the  test  machine  with  the 

support  strips  cut  (Reference  [8]). 


Specimen  fabrication 

The  3.5pm  thick,  P-doped  poly-Si  thin  film  specimens  are  manufactured  at  MCNC 
using  the  so  called  MUMPs  surface  micromachining  process.  They  are  released  and 
prepared  for  testing  at  Johns  Hopkins  University.  The  steps  needed  to  fabricate  the 
specimen  are  quite  involved,  and  the  details  are  not  repeated  here.  They  can  be  recovered 
in  References  [7]-[8].  Only  a  brief  description  of  the  material  and  the  back-etching 
procedure  are  described  next. 

It  is  important  to  note  that  MEMS  devices  fabricated  at  MCNC  are  comprised  of 
two  layers  of  poly-Si,  deposited  by  LPCVD,  each  contained  within  sacrificial  layers  of 
PSG.  The  bottom  layer  of  PSG  is  separated  from  a  supporting  silicon  substrate  by  a  layer 
of  silicon  nitride.  The  as-deposited  films  are  amorphous;  a  fine  grain,  low  internal  stress, 
polycrystalline  microstructure  is  achieved  by  annealing  for  one  hour  at  1050°C  (each  layer 
is  annealed  separately).  As  shown  in  the  cross-sectional  TEM  micrograph  of  Figure  5,  the 
thickness  of  the  first  (second)  deposit  is  2pm  (1.5pm).  It  is  observed  that  this  final 
structure  is  nearly  columnar.  The  poly-Si  is  doped  during  the  annealing  from  the 
surrounding  PSG.  In  Reference  [8]  it  is  reported  that  the  PSG  layer  that  separates  the 
poly-Si  layers  was  removed  during  processing,  and  that  the  second  polysilicon  layer  was 
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deposited  directly  on  top  of  the  first. 


Figure  5-  TEM  photograph  of  the  two-layer  polysilicon  specimen.  The  bottom  layer 
is  2.0  pm  thick  and  the  darker  top  layer  is  1.5  pm  thick  (Reference  [8]). 

Figure  6a  shows  a  schematic  of  the  specimen  contained  within  a  one-centimer 
square  carrier  frame  that  includes  the  grip  ends  and  the  support  strips.  Figure  6b  shows  the 
evolution  of  section  A-A,  from  the  as  received  configuration  to  the  final  released 
configuration.  The  second  (from  the  top)  schematic  in  Figure  6b  shows  the  ceramic  carrier 
that  was  attached,  for  handling  purposes,  to  the  front  of  the  die  using  an  adhesive  wax. 
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Gold  Lines 


(d)  Cross-section  of  completed  tensile  specimen. 


(b) 


Figure  6-  Schematic  of;  (a)  poly-Si  pattern,  (b)  back-etch  process  as  viewed  at 


section  A-A  of  (a)  (Reference  [7]). 


An  SEM  micrograph  showing  an  oblique  view  of  the  specimen,  taken  from 
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Reference  [7],  is  shown  in  Figure  7. 


Results 

Young’s  modulus,  Poisson’s  ratio,  and  tensile  strength 

Typical  stress-strain  curves  are  shown  in  Figure  8.  As  stated  by  the  authors  in 
Reference  [8],  the  specimen  buckles  out-of  plane  when  the  support  strips  are  cut,  and 
therefore  a  preload  is  applied  before  initiation  of  displacement  controlled  loading.  The 
figure  indicates  a  tensile  strength  of  1.2GPa,  Young’s  modulus  of  169GPa  and  Poisson’s 
ratio  of  0.22.  This  latter  measurement  is  the  first  on  poly-Si. 


40 


BIAXIAL  TEST  #18 


Figure  8-  A  typical  stress  biaxial-strain  plot  (Reference  [8]). 

Table  1  presents  the  data  for  48  tests  on  specimens  produced  during  five  MUMPs 
production  runs.  Ten  specimens  were  tested  from  MUMPs  6,  10,  and  11,  and  nine  were 
tested  from  MUMPs  8  and  12.  Note  that  it  is  this  data  that  is  reported  in  the  previous 
chapter  describing  variations  in  mechanical  properties  of  polycrystalline  silicon. 


Table  1-  Mechanical  property  data  for  polysilicon  (Reference  [8]). 


MUMPs 

Young’s 

Modulus 

(GPa) 

Poisson’s 

Ratio 

Tensile 

Strength 

(GPa) 

6 

163  ±4.78 

- 

1.23  ±  .200 

8 

173  ±  6.94 

- 

1.27  ±  .107 

10 

173  ±  2.00 

- 

1.14  ±.155 

11 

168  ±7.15 

0.22  ±.015 

1.19  ±  .155 

12 

168  ±  1.59 

0.22  ±.003 

1.17  ±.100 

All  Data 

169  ±6.15 

0.22  ±  .011 

1.20  ±  .150 

The  statistical  analyses  performed  by  the  authors  showed  that  (1)  at  95% 
confidence  the  average  tensile  strength  of  1.2GPa  is  statistically  equal  to  the  average 
strength  of  the  individual  runs,  (2)  at  95%  confidence  the  average  Young’s  modulus  of 
169GPa  is  statistically  different  than  the  average  of  MUMPs  6  and  of  MUMPs  10,  (3)  the 
coefficient  of  variation  for  all  of  the  Young’s  modulus  data,  all  of  the  strength  data,  and  all 
of  the  Poisson’s  ratio  data  are,  respectively,  3.6%,  12.5%  and  5%.  All  and  all  the  results 
suggest  that  the  tests  and  the  fabrication  processes  are  very  consistent  and  reproducible. 


Fracture  toughness 

In  Reference  [14]  the  authors  used  the  experimental  procedure  to  measure  the 


41 


fracture  toughness  of  poly-Si,  from  the  MUMPs  13  run,  using  the  specimen  shown 
schematically  in  Figure  9.  The  fabrication  process  is  not  reported;  it  is  similar  to  the  one 
described  previously  for  the  tensile  specimens.  It  should  be  mentioned  that  this  group  tried 
to  use  the  first  generation  device  to  measure  fracture  toughness,  but  was  not  successful 
(Reference  [7]). 

This  fracture  mechanics  configuration  consists  of  a  300Qum  wide  by  4000//m  long 
strip  which  contains  a  lOO/rni  central  crack.  The  relevant  dimensions  of  the  crack  are 
shown  in  the  figure.  The  crack-tip  radius  is  reported  in  Reference  [14]  as  1/mi.  In  these 
experiments  the  ISGD  system  is  used,  with  the  gold  pads  shown  in  Figure  2,  to  measure 
the  COD  at  the  center  of  the  crack.  Note  that  because  the  pads  are  at  a  finite  distance 
from  the  crack  surfaces,  the  relative  displacement  that  is  measured  is  not  exactly  the  COD. 
A  photograph  of  the  center  crack  and  the  gold  pads  is  shown  in  Figure  10. 


Figure  10-  Close-up  photograph  of  the  center  crack  and  the  gold  pads  (Reference 

[14]). 

If  the  crack  in  this  configuration  were  mathematically  sharp,  the  stress  intensity 
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factor  and  plane  stress  COD  at  the  center  of  the  crack  would  be 


(la) 

(lb) 


A 


COD  =  2 


Pa 

AE 


where  P,  a  and  A  are,  respectively,  the  applied  force,  half-crack  length,  and  cross-sectional 
area.  Young’s  modulus  in  Equation  (lb)  was  assumed  to  be  169GPa;  this  value 
corresponds  to  the  average  value  measured  in  the  tests  described  in  the  previous  section. 

The  results  of  the  12  tests  are  presented  in  Figure  1 1  (COD)  and  Table  2  (fracture 
toughness).  It  is  observed  that  the  agreement  between  the  theoretical  and  measured  COD 
is  reasonable.  The  differences  are  due  to  the  fact  that  the  crack  is  of  finite  thickness  and 
not  sharp,  and  that  the  COD  is  measured  away  from  the  crack  surfaces. 


COD  -  microns 

Figure  11-  Measured  stress  versus  COD  and  comparison  with  theory  (Reference 

[14]). 


The  average  fracture  toughness,  its  standard  deviation,  and  its  coefficient  of 
variation  are,  respectively,  1.41MPa-m1/2,  ±0.15  MPa-m1/2,  and  11%.  The  fracture 
toughness  of  single  crystal  silicon  is  anisotropic,  its  range  typically  reported  as  0.7- 1.2 
MPa-mI/2  (ft.  Ballarini,  R.L.  Mullen,  Y.Yin,  H.  Kahn,  S.  Stemmer  and  A.H.  Heuer, 
“ The  fracture  toughness  of  polysilicon  microdevices:  a  first  report,”  Journal  of 
Materials  Research,  Vol.  12,  No.  4,  April  1997). 

The  values  in  Table  2  may  not  represent  the  true  fracture  toughness  because  the 
notch  tip  is  of  finite  radius.  However,  since  this  radius  is  much  smaller  than  the  other 
dimensions,  the  values  are  expected  to  be  closer  to  the  true  toughness  than  those  reported 
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by  Tsuchiya  (in  the  previous  chapter)  and  by  Ballarini  and  coworkers  (whose  work  is 
described  in  the  next  chapter). 


[  Test 

F2 

F3 

F9 

F10 

Fll 

F12 

F13 

F14 

F15 

F16 

F17 

F18 

1.14 

1.19 

1.32 

1.49 

1.51 

1.52 

1.51 

1.28 

1.53 

1.38 

1.42 

1.61 

Table  2-  Fracture  toughness  (MPa-m1/2)  (Reference  [14]). 

3.  SECOND  GENERATION  GRIPPING  SYSTEM 
General  description  of  experimental  technique 

Because  of  the  previously  discussed  limitations  of  the  first  generation  experimental 
procedure,  the  authors  replaced  the  adhesive  based  gripping  system  with  an  electrostatic 
one  similar  to  the  one  developed  at  TCRDL  and  described  in  the  previous  chapter.  A 
schematic  of  the  second  generation  environment  is  shown  in  Figure  12.  One  end  of  the 
specimen,  which  is  described  subsequently,  is  attached  to  a  550/zm  wide  probe  through 
electrostatic  forces.  This  probe,  which  is  cut  from  a  silicon  chip  and  insulated  by  a  600  nm 
thick  layer  of  nitride,  is  in  turn  mounted  to  a  load  cell  (as  reported  in  Reference  [15],  two 
load  cells  were  used,  with  ranges  of  5  and  100  grams).  Alignment  of  the  specimen  is 
achieved  by  mounting  the  load  cell  and  probe  on  a  manual  three-axis  translation  stage.  The 
other  end  of  the  specimen  is  attached  to  a  five-axis  piezoelectric  stage  (with  a  resolution 
of  20  nm,  as  reported  in  Reference  [15]),  which  allows  precise  alignment,  and  whose 
movement  produces  a  tensile  force  on  the  specimen. 

Two  methods  of  strain  measurement  are  applied;  both  correct  for  the  stiffness  of 
the  load  cell  On  the  narrower  specimens  the  movement  of  the  piezoelectric  stage  is 
monitored  with  a  100/xm  range  capacitance  gage.  On  the  wider  specimens,  the  ISGD 
system  is  also  used. 


Laser 


Base  PJate 
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Figure  12-  Schematic  of  testing  environment  (Reference  [15]). 

Specimen  fabrication 

The  specimens  were  fabricated  on  the  MUMPs  19  and  21  runs.  Figure  13  shows 
the  18  specimens  that  come  on  each  one  square  centimeter  die.  Figure  14  shows  a 
schematic  of  an  individual  specimen,  showing  the  dimensions  of  the  free  end  (which  is 
attached  to  the  probe)  and  the  fixed  end.  Note  that  the  600/an  width  of  the  free  end  is 
slightly  larger  than  the  550/im  probe  width.  Apparently  this  difference  made  it  more 
convenient  to  align  the  specimen.  Size  effects  were  studied  by  varying  the  thickness, 
width,  and  length  of  the  specimens.  One  row  of  specimens  was  made  1.5//m  thick,  the 
other  3.0/mi  thick.  Each  row  contained  one  50/mi  wide  (in  the  test  section)  specimen,  two 
that  were  20//m  wide,  three  5/im  wide,  and  three  2/zm  wide.  Three  specimens  are  500/im 
long  (in  the  straight  section),  three  are  200//m  long,  two  are  50/zm  long,  and  one  is  20/zm 
long.  The  dimensions  of  the  smallest  specimen,  which  was  1.5x2x20^m,  are  comparable  to 
those  of  a  typical  MEMS  component  (an  attachment  beam,  for  example). 


Figure  13-  Eighteen  specimens  on  a  one  square  centimeter  die  (Reference  [15]). 


Free  End  with 
Etch  Holes 


Fixed 

End 


Figure  14-  Schematic  of  the  specimen,  with  dimensions  in  fim.  The  length,  thickness, 
and  width  vary,  as  described  in  the  text  (Reference  [15]). 
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Results 


As  stated  by  the  leading  author  at  the  MRS  1998  Spring  meeting,  because  the 
second  generation  system  is  relatively  new,  the  results  presented  next  are  preliminary. 

Young’s  modulus  was  calculated  from  the  specimen’s  force-displacement  plot;  for 
the  wider  specimens  it  was  also  calculated  from  strain  calculated  from  the  relative 
displacements  measured  using  the  ISGD  system.  Using  these  two  methods  proved 
important,  as  explained  next.  Typical  plots  for  two  specimens  with  identical  dimensions  in 
the  plane  of  the  film  (20pm  x  500pm),  but  of  different  thicknesses  (1.5/un  versus  3.5pm), 
are  shown  in  Figure  15.  It  is  observed  that,  for  the  thicker  specimen,  the  values  of 
modulus  determined  using  the  two  methods  is  practically  the  same.  For  the  thinner 
specimen,  on  the  other  hand,  the  values  are  significantly  different.  This  discrepancy  is 
currently  being  investigated. 


FORCE-DISPLACEMENT 


STRESS-STRAIN 


Figure  15-  Force-displacement  (left)  and  stress-strain  (right)  plots  of  two  specimens 
identical  in  shape,  but  of  different  thickness  (Reference  [15]). 


Figure  16  presents  the  measured  Young’s  modulus  and  tensile  strength  of  all  the 
specimens.  Those  for  the  600pm  wide  specimens  were  taken  from  Reference  [8]  and 
included  for  comparison  purposes.  It  is  observed  that  the  scatter  in  the  Young’s  modulus 
measurements  is  significant.  At  this  point  it  may  be  premature,  but  it  seems  that  the 
strength  exhibits  the  size  effect  typical  of  brittle  materials;  the  widest  specimens  from 
Reference  [8]  are,  on  the  average,  the  weakest. 
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Figure  16-  Modulus  and  strength  data  for  specimens  of  various  sizes  (Reference 

[15]). 
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4.  CRITIQUE 


The  first  generation  testing  environment  introduced  the  use  of  laser  interferometry 
to  accurately  measure  strains  in  poly-Si  processed  using  MEMS  technology.  It  was  used 
to  produce  very  valuable  data  for  elastic  modulus,  Poisson’s  ratio,  and  the  direct  tensile 
strength  of  poly-Si  thin  films. 

The  second  generation  system,  which  combines  the  use  of  electrostatic  gripping 
with  laser  interferometry,  allows  measurements  on  specimens  whose  dimensions  are 
comparable  to  MEMS  devices.  In  the  scientist’s  opinion,  when  the  system  described  in  this 
chapter  is  fine-tuned,  it  will  be  the  most  robust  testing  environment  for  measuring  the 
elastic  modulus  and  direct  tensile  strength  of  poly-Si  films. 
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CHAPTER  5 


BERKELEY  SENSOR  AND  ACTUATOR  CENTER 
1.  INTRODUCTION 

This  section  describes  the  experiments  conducted  by  P.T.  Jones,  G.C.  Johnson  and 
R.T.  Howe  at  the  University  of  California  at  Berkeley  (UCB),  to  measure  the  tensile 
strength  and  Young’s  modulus  of  poly-Si  at  scales  relevant  to  microelectromechamcal 
devices.  Their  work  will  appear  as 

[1]  P.T.  Jones,  G.C.  Johnson  and  R.T.  Howe,  “Fracture  strength  of 
polycrystalline  silicon ,”  to  appear  in  Proceedings  of  the  1998  MRS  Spring  Meeting, 
San  Francisco. 

This  group  tested  P-doped  specimens  from  MCNC’s  MUMPsl9  and  MUMPs21 
production  runs,  and  undoped  specimens  produced  at  UCB.  Their  tests  on  the  MCNC 
material  represents  their  contribution  to  the  previously  mentioned  round-robin  program 
for  measuring  tensile  strength  and  Young’s  modulus  of  poly-Si  films.  A  comparison  of  the 
results  for  the  two  different  MCNC  MUMPs  runs  showed  that  there  were  no  run-to-run 
variations  of  mechanical  properties  using  the  nominally  identical  process.  As  expected,  a 
comparison  of  the  results  for  the  MCNC  material  and  the  Berkeley  material  showed  that 
differences  in  processing  can  result  in  variations  of  mechanical  properties. 


Figure  1:  SEM  micrograph  of  fracture  structure  test  device.  The  six  cantilever 
beams  to  be  tested  are  attached  to  a  shuttle  which  is  anchored  by  a  folded  flexure 

structure  (Reference  [1]). 
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This  group’s  device  determines  the  fracture  (tensile)  strain  of  poly-Si  films  by 
measuring  the  displacement  at  fracture  of  a  series  of  cantilever  beams  bending  in  the  plane 
of  the  film.  The  displacement  is  related  to  the  strain  through  the  analytical  solution  of  the 
elastica.  Concomitantly  a  different  device,  consisting  of  a  series  of  beams  of  different 
length  bending  out  of  the  plane  of  the  film,  is  used  to  measure  Young’s  modulus. 
Assuming  the  films  are  isotropic,  the  fracture  strain  and  Young’s  modulus  are  used  to 
determine  the  (bending)  tensile  strength. 

2.  GENERAL  DESCRIPTION  OF  EXPERIMENTAL  PROCEDURE 

Figure  1  shows  a  SEM  micrograph  of  the  500  /xm  x  320  yum  test  structure  used  to 
measure  the  fracture  strain.  Six  cantilever  beams  (three  different  length  pairs)  are  attached 
to  a  shuttle,  which  in  turn  is  anchored  to  the  substrate  by  a  folded-flexure  frame.  The 
beams  are  nominally  2  yum  thick,  3  pm  wide,  and  have  lengths  of  50  p m,  60  pm,  or  70  pm. 
The  authors  measured  the  thickness  and  width,  respectively,  as  1.9  yum  and  2.8  pm 
( private  communication  with  P.T.  Jones).  The  frame  is  designed  to  permit  a  relatively 
large  displacement  of  the  beams  in  the  plane  of  the  film,  and  relatively  small  in-plane 
rotation  and  out-of-plane  displacement  of  the  assembly. 

The  authors  mention  the  following  advantages  of  this  device  over  uniaxial 
structures. 

1.  Its  small  size  allows  a  large  number  of  devices  to  be  fabricated  on  a  single  die. 
This,  coupled  with  the  fact  that  as  many  as  six  data  points  could  be  obtained 
from  a  single  device,  allows  a  large  amount  of  statistical  data  to  be  obtained. 

2.  The  state  of  strain  in  the  beam  is  independent  of  the  thickness  of  the  film; 
therefore  the  thickness  need  not  be  measured.  Moreover,  this  allows 
application  of  the  test  device  to  the  study  of  bending  of  other  materials  such  as 
deep-etched  single  crystal  silicon,  or  even  hexsil  structures  ( C.G .  Keller  and 
R.T.  Howe,  “ Hexsil  tweezers  for  teleoperated  micro-assembly,”  Proceedings 
IEEE.  The  Tenth  Annual  International  Workshop  on  Micro  Electro 
Mechanical  Systems.  An  Investigation  of  Micro  Structures,  Sensors, 
Actuators,  Machines  and  Robots,  Nagoya,  Japan,  26-30  Jan.  1997.  New 
York,  NY,  USA:  IEEE,  1997.  p.  72-7). 

3.  The  beam  can  achieve  large  deformations  before  fracture.  This  is  an  advantage 
since  the  displacement  measurement  resolution  is  0.2 pm. 

4.  The  dimensions  of  the  specimens  are  comparable  to  those  of  typical  MEMS 
components. 

5.  The  structural  layer  of  the  device  can  be  defined  using  only  a  single  mask.  It  is 
noted  in  Reference  [1]  that  the  tested  devices  were  fabricated  using  a  minimum 
of  two  masks;  one  defines  the  anchor  region,  the  other  defines  the  structural 
layer. 

The  authors  also  note  that  the  force  on  the  beam  cannot  be  directly  measured.  This 
necessitates,  as  explained  subsequently,  a  different  device  to  determine  the  Young’s 
modulus. 
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Figure  2-  Scanned  video  image  of  fracture  structure  test  device.  An  off  chip  probe 
loads  the  shuttle  of  the  device.  Here  the  shortest  two  beams  have  already  broken  off 
as  other  beams  continue  to  be  bent  (Reference  [1]). 

As  shown  in  Figure  2,  the  beams  are  loaded  when  they  come  into  contact  with 
anchored  square  structures  as  the  shuttle  is  pushed,  at  a  continuous  rate  of  0.025  fim/s,  by 
an  off-chip  probe.  The  shuttle  is  marked  with  a  vernier  scale.  A  video  system  monitors  to  a 
0.2  resolution  the  shuttle’s  displacement,  and  the  vernier  is  used  to  determine  the 
fracture  displacement  of  each  beam.  This  latter  displacement  is  used  in  conjunction  with 
the  elastica  solution  to  determine  the  fracture  strain.  The  analysis,  which  is  standard  and 
not  presented  here,  appears  in  P.T.  Jones,  G.C.  Johnson,  “ Micromechanical  Structures 
for  Fracture  Testing  of  Brittle  Thin  Films,”  ASME  International  Mechanical 
Engineering  Congress  and  Exposition  Micro-Electro-Mechanical  Systems  (MEMS), 
DSC-Vol.  59,  p.  325-330, 1996.  A  close-up  view  of  one  beam  is  shown  in  Figure  3,  which 
clearly  demonstrates  that  the  end  of  the  beam  undergoes  a  relatively  large  displacement 

The  following  uncertainty  arises.  The  design  is  symmetric  with  respect  to  the 
direction  of  loading,  and  as  a  result,  two  beams  of  equal  length  are  tested  simultaneously. 
It  is  unlikely  that  both  fracture  at  the  same  exact  time.  After  the  first  beam  breaks, 
symmetry  is  lost,  and  fracture  of  the  second  beam  may  occur  at  a  different  strain  level  The 
authors  considered  this  issue  (private  communication  with  P.T.  Jones).  They  did  not 
observe  rotation  of  the  shuttle  during  testing  of  the  shortest  four  beams.  In  some  tests, 
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only  one  of  the  last  of  the  longest  beams  remains  intact,  and  the  shuttle  shifts  upwards;  the 
results  from  these  tests  are  discarded.  As  a  check,  they  compared  the  failure  strains  of  the 
longest  beams  with  those  of  the  shorter  beams,  and  found  a  close  match. 


Figure  3-  Scanned  video  image  of  a  2.8  //m  wide  by  1.9  fim  thick  beam  at  a  large 
displacement  The  distance  from  the  fixed  end  of  the  beam  to  the  contact  point  is  70 
fim.  The  beam  has  deflected  21.6  fim  (Reference  [1]). 

3.  SPECIMEN  FABRICATION 

The  above  referenced  paper  does  not  provide  the  details  of  the  fabrication 
processes.  These  could  easily  be  recovered  from  UCB  and  MCNC.  What  is  described  in 
the  following. 

The  MCNC  films  were  deposited  at  610°C,  and  then  sandwich-P-doped  (as 
described  in  the  JHU  summary)  during  a  1  hour  furnace  anneal  at  1050°C.  The  structures 
released  from  MUMPs  19  were  both  furnace  dried  and  super-critical  carbon  dioxide  dried 
(G.T  Mulhem,  D.S  Soane  and  R.T.  Howe,  “ Supercritical  carbon  dioxide  drying  of 
microstructures,”  Transducers  ‘93,  pp.  296-299, 1993). 

The  Berkeley  undoped  poly-Si  was  deposited  at  610°C  and  300  mTorr, 
sandwiched  by  LTO  and  furnace  annealed  at  950°C  for  Vi  hour. 

4.  FRACTURE  TENSILE  STRAIN  RESULTS 

The  measured  average  displacement  at  fracture  and  average  strain  at  fracture  are 
shown,  respectively,  in  Tables  I  and  II  (these  were  reproduced  from  Reference  [1]). 
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Table  I-  Average  Beam  Displacement  (pm)  at  Fracture  for  the  Materials  Tested. 
Uncertainties  are  one  standard  deviation  of  the  measured  displacement  (Reference 

[1]). 


mmjmm 

liiiSl 

MUMPS  19 
C02  Dry 

BSAC 

50  pm  beams 

13.5  ±  1.2 

13.1  ±  1.4 

13.2  ±  1.6 

10.4  ±0.9 

60  pm  beams 

18.6  ±2.3 

18.6  ±2.5 

17.7  ±  1.7 

14.7  ±  2.0 

70  pm  beams 

25.4  ±4.1 

27.1  ±3.6 

25.7  ±  2.4 

20.1  ±2.3 

Beams  Tested 

92 

64 

68 

76 

Table  II-  Estimated  Maximum  Strain  State  in  Beams  at  Failure  (%) 

(Reference  [1]). 


MmMim 

mBSSSBM 

BSAC 

Average 

1.5  ±0.2 

1.5  ±0.2 

1.5  ±0.2 

1.2  ±0.1 

Minimum 

0.81 

1.10 

1.18 

0.90 

Maximum 

1.94 

1.91 

1.89 

1.73 

The  authors’  results  and  conclusions  are  summarized  as  follows. 

1.  A  Student’s  t-test  of  the  data  suggests  that  the  fracture  strain  of  the  MUMPs 
poly-Si  is  statistically  different  than  that  of  the  Berkeley  poly-Si.  The  authors 
attribute  this  to  differences  in  processing  (including  the  fact  that  one  material  is 
doped,  the  other  undoped),  which  results,  among  other  things,  in  different 
levels  of  side-wall  roughness. 

2.  No  statistical  difference  exists  between  the  fracture  strains  of  (MCNC) 
specimens  that  were  super-critically  carbon  dioxide  dried  and  those  that  were 
furnace  dried. 

3.  All  MUMPs  specimens  have  the  same  average  fracture  strain. 

5.  YOUNG’S  MODULUS  DEVICE  AND  RESULTS 

The  force  applied  to  the  specimen  in  the  fracture  strain  device  cannot  be  measured 
directly.  Therefore  the  authors  used  the  device  suggested  by  Tai  (Y.C.  Tai,  “IC- 
Processed  Polysilicon  Micromechanics:  Technology,  Material,  and  Devices,” 
Dissertation  -  University  of  California  at  Berkeley,  1989 )  to  determine  Young’s  modulus 
of  the  MUMPs21  poly-Si  films.  As  shown  in  Figure  4,  it  is  comprised  of  an  array  of 
relatively  wide  cantilever  beams  15  pm  wide,  2  pm  thick.  Their  lengths  vried  from  25  to 
1000  pm  in  5  pm  steps.  However,  most  of  the  beams  were  too  compliant  for  testing 
( private  communication  with  P.T.  Jones).  A  Dektak  3030  profiler  is  used  to  measure  the 
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out-of-plane  displacement  of  the  free  ends  of  the  cantilever  beams  for  a  variety  of  applied 
loads.  The  Young’s  modulus  E  may  be  determined  from  the  equation 


(1) 


E  = 


4(1 -v2)  PL3 
Sbh 3 


where  v  is  Poisson’s  ratio  (assumed  to  be  0.25),  P  is  the  applied  end  load,  L  is  the  length 
of  the  beam,  8  is  the  end  displacement  of  the  beam,  b  is  the  width  of  the  plate,  and  h  is  the 
film  thickness.  The  displacement  of  the  end  of  the  beam  is  measured  relative  to  small 
anchored  plates  on  either  side  of  the  beam. 

It  is  important  to  note  that  the  modulus  determined  using  this  device  is  associated 
with  bending  out  the  plane  of  the  film,  while  the  fracture  strain  specimens  bend  in  the 
plane.  There  is  no  evidence  that  suggests  the  films  tested  are  isotropic.  Indeed,  the  MCNC 
films  are  typically  columnar,  as  described  in  the  work  at  JHU,  and  probably  transversely 
isotropic.  This  point  will  apparently  be  discussed  further  in  the  summary  of  the  round- 
robin  tests  (private  communication  with  P.T.  Jones). 


Figure  4-  SEM  of  cantilever  “beams”  used  for  determining  Young’s  modulus 

(Reference  [1]). 
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Plots  of  end  displacement  versus  beam  length,  for  three  different  applied  forces, 
are  shown  in  Figure  5.  These  plots  show  that  the  displacement  obeys  Equation  (1)  if 
Young’s  modulus  is  set  equal  to  173  ±  20  GPa,  where  the  uncertainty  is  a  single  standard 
deviation  of  the  measured  values. 

Using  these  values  of  the  modulus,  and  assuming  that  the  modulus  measured  is 
equal  to  the  in-plane  modulus,  the  average  tensile  strength  becomes  2.1GPa  and  2.6GPa 
for  the  MUMPS  and  Berkeley  films,  respectively.  However,  this  also  assumes  that  the 
material  remains  linear  elastic  up  to  the  maximum  strain.  This  point  was  brought  up  by 
Arthur  H.  Heuer  (of  CWRU)  after  the  authors  presented  this  work  at  the  MRS  Spring 
1998  meeting.  There  is  no  evidence  that  the  material  is  linear  elastic  at  the  strain  levels 
measured  in  the  experiments  described  above.  Silicon  whiskers  in  bending  do  remain 
elastic  at  room  temperature  up  to  failure  strains  approximately  equal  to  1%,  as  shown  in 
Figure  6  ( G.L .  Pearson,  W.T.  Read,  Jr.  And  W.L.  Feldmann,  “ Deformation  and 
fracture  of  small  silicon  crystals,”  Acta  Metallurgica,  Vol.  5,  pp.  181-191, 1957). 


Figure  5-  Plot  of  displacement  vs.  beam  length  for  three  different  force  settings  using 
the  Young’s  modulus  test  device  (Reference  [1]). 


Fig.  4.  Load -deflection  corvee  for  eilieon  whiskers 
at  several  temperatures* 


Figure  6-  Load-deflection  curves  for  silicon  whiskers  at  several  temperatures 

(Reference  provided  in  text). 


6.  CRITIQUE 

The  testing  environment  described  in  this  chapter  relies  on  relatively  simple  (in¬ 
plane  flexural)  loading  and  optical  observation  systems  to  measure  (in  a  very  reproducible 
way)  the  fracture  strain  of  thin  poly-Si  films.  For  measuring  fracture  strain,  the  system 
appears  to  be  sufficiently  robust.  However,  the  beam  bending  system  used  to  determine 
the  Young’s  modulus  suffers  from  the  fact  that  the  bending  is  out  of  the  plane  of  the  film, 
and  thus  measures  a  modulus  which  may  not  be  equal  to  the  one  associated  with  the  in¬ 
plane  bending  specimens. 
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CHAPTER  6 


EXPONENT  ANALYSIS 


1.  INTRODUCTON 

This  chapter  summarizes  the  experimental  research  performed  by  S.B.  Brown,  W. 
Van  Arsdell  and  C.L.  Muhlstein  at  Exponent  (formerly  Failure  Analysis  Associates)  on 
environmentally  assisted  fatigue  crack  growth  in  poly-Si.  Their  work  will  appear  as 

[1]  S.B.  Brown,  W.  Van  Arsdell  and  C.L.  Muhlstein,  “Materials  reliability 
in  MEMS  devices,”  Transducers  1997,  1997  International  Conference  on  Solid-State 
Sensors  and  Actuators,  Chicago,  Illinois,  June  16-19, 1997. 

This  group  has  demonstrated  that  poly-Si  (from  the  MCNC  MUMPs  13  run)  is 
associated  with  stress  corrosion  cracking  in  the  presence  of  moisture.  This  is  an  important 
result  that  has  direct  implications  on  the  long  term  reliability  of  MEMS  devices. 

2.  DESCRIPTION  OF  EXPERIMENTAL  PROCEDURE 

This  group  developed  the  in-plane  resonant  fatigue  device  shown  in  Figure  1  to 
study  slow  crack  growth  in  poly-Si.  The  notched  cantilever  specimen  is  attached  to  a 
wedge  shaped  plate  that  serves  as  the  resonant  mass.  Two  sets  of  comb  drive  structures 
are  formed  along  the  two  long  straight  sides  of  this  mass.  One  provides  electrostatic 
actuation,  the  other  capacitance  sensing  of  motion.  A  close-up  view  of  the  notched 
cantilever  specimen  is  shown  in  Figure  2. 

The  dimensions  of  the  device  are  not  provided  in  the  above  referenced  paper, 
because  the  paper  focuses  on  the  demonstration  that  poly-Si  is  associated  with  stress 
corrosion  cracking.  Quantitative  measurements  of  stress  or  strain  levels  during  the 
experiments  are  not  reported.  From  the  micrographs  it  is  observed  that  the  device  and 
specimen  are,  respectively,  on  the  order  of  hundreds  and  tens  of  microns. 

The  authors  do  state  that  the  geometry  of  the  notch  was  selected  to  ensure  that  at 
resonance  the  specimen  fractures  immediately  (a  very  short  amount  of  time  is  obviously 
needed  to  reach  resonance).  Subsequent  plots  of  the  data  are  normalized  with  respect  to 
the  excitation  voltage  required  to  break  the  specimen  immediately. 

Evidence  of  crack  growth  is  searched  for  by  measuring  variations  in  time  of  the 
natural  frequency  of  the  device.  (The  frequency  depends  on  the  stiffness  of  the  structure; 
as  the  crack  grows  the  stiffness  decreases).  S.  Brown  has  used  this  strategy  (with  different 
devices)  to  demonstrate  environmentally  assisted  fatigue  crack  growth  in  single  crystal 
silicon  and  nickel  (J.A.  Connally  and  S.B.  Brown,  “Slow  crack  growth  in  single-crystal 
silicon,”  Science,  Vol.  256,  pp.  1537-1539,  12  June,  1992;  J.A.  Connally  and  S.B. 
Brown,  “Micromechanical  fatigue  testing,”  Experimental  Mechanics,  pp.  81-90,  June 
1993;  S.B.  Brown,  G.  Povirk  and  J.  Connally,  “Measurement  of  slow  crack  growth  in 
silicon  and  nickel  micromechanical  devices,  ”  Proceedings,  MEMS-93,  Fort 
Lauderdale,  Florida,  pp.  99-104,  February  7-10,  1994).  The  base  line  (first  mode) 


frequency  is  found  by  sweeping  a  range  of  frequencies  around  the  expected  response,  and 
observing  the  displacement  amplitude.  A  second  order  polynomial  is  used  to  fit  the  peak, 
and  in  turn  the  maximum  amplitude  is  calculated.  Changes  in  resonant  frequency  are 
measured  as  follows.  For  a  period  of  time,  a  defined  voltage  is  used  to  excite  the 
specimen  at  the  peak  frequency.  Re-evaluation  of  the  frequency  response  by  sweeping  at 
the  excitation  frequency  permits  measuring  of  the  changes. 

3.  RESULTS 

The  experiments  were  conducted  at  75%  relative  humidity  and  27°C.  The  results 
are  presented  in  Figure  3.  This  plot  is  a  slight  modification  of  the  traditional  way  of 
presenting  fatigue  data.  The  vertical  axis  typically  presents  the  applied  stress.  To 
eliminate  the  effects  of  uncertainties  in  the  displacement  measurements,  the  modulus 
required  to  relate  displacement  to  stress,  and  other  parameters,  the  vertical  axis  was 
chosen  as  the  excitation,  normalized  with  respect  to  the  immediate  fracture  excitation. 

The  data  clearly  shows  that  slow  crack  growth  occurs.  The  only  question  that 
arises  is  whether  the  growth  is  due  to  stress,  or  to  humidity.  The  authors  claim  that  they 
obtained  preliminary  results  (not  presented  in  Reference  [1])  that  showed  crack  growth 


Figure  1-  In-plane  resonant  fatigue  device  (Reference  [1]). 
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Figure  2-  Close-up  view  of  notched  cantilever  specimen  (Reference  [1]). 


Time  to  Failure,  Cycles 
io9  1010 


Time  to  Failure,  minutes 

Figure  3-  Time  to  failure  versus  normalized  excitation  (Reference  [1]). 

requires  moisture.  Furthermore,  they  attribute  the  phenomena  to  stress  corrosion  of  the 
native  oxide  on  the  surface  of  the  specimen. 
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As  will  be  discussed  in  the  next  chapter  (which  summarizes  the  work  performed 
by  the  scientist  and  coworkers  at  Case  Western  Reserve  University),  there  is  some 
evidence  that  suggests  cyclic  loads  could  reduce  the  strength  of  poly-Si  in  the  absence  of 
moisture. 

4.  CRITIQUE 

The  experiments  described  above  are  very  important;  they’ve  shown  that  poly-Si 
is  susceptible  to  environmentally  assisted  subcritical  crack  growth.  This  has  significant 
implications  for  long  term  durability  of  MEMS  devices. 
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CHAPTER  7 


CASE  WESTERN  RESERVE  UNIVERSITY 
1.  INTRODUCTION 

This  chapter  summarizes  the  development,  by  the  author  and  his  coworkers  at 
Case  Western  Reserve  University,  of  two  specimens  for  measuring  the  toughness  of  poly- 
Si  films.  This  group’s  ultimate  goal  is  to  design  and  fabricate  a  robust  electrostatically 
loaded  test  specimens  that  will  allow  on-chip  testing  of  various  mechanical  properties  of 
MEMS  devices.  The  first  design  involved  electrostatically  loaded  deeply  cracked  double 
cantilever  beam  specimens  that  failed  to  produce  a  force  sufficient  for  fracture.  This  design 
is  not  being  used  anymore.  However,  in  order  to  generate  preliminary  toughness  data  that 
could  be  used  in  the  second  design,  these  specimens  were  wedge-loaded  to  fracture  by  an 
external  probe.  The  results,  which  will  be  summarized  in  the  next  section,  appear  in 

[1]  R.  Ballarini,  R.L.  Mullen,  Y.  Yin,  H.  Kahn,  S.  Stemmer  and  A.H. 

Heuer,  “ The  fracture  toughness  of  polycrystaUine  silicon  microdevices:  a  first 

report,”  Journal  of  Materials  Research,  Vol.  12,  No.  4,  pp.  915-922, 1997. 

The  second  electrostatic  actuator  design  was  much  more  successful  It  produced 
enough  force  to  fracture  films  associated  with  a  residual  compressive  stress,  and  produced 
enough  displacement  under  resonance  conditions  to  fracture  those  under  a  residual  tensile 
stress.  This  work  has  been  submitted  for  publication  as 

[2]  H.  Kahn,  R.  Ballarini,  R.L.  Mullen  and  A.H.  Heuer,  “Electrostatically 
actuated  failure  of  microfabricated  polysilicon  fracture  mechanics  specimens,” 
submitted  to  Proceedings  of  the  Royal  Society  of  London,  1998. 

The  geometry  of  this  second  device  is  currently  being  fine  tuned  to  increase  the 
force  transmitted  from  the  actuator  to  the  fracture  specimen,  and  to  mitigate  certain  out- 
of-plane  displacements  that  result  from  residual  stresses. 

Before  proceeding  with  a  summary  of  the  above  referenced  papers,  it  is  important 
to  note  that  the  polycrystalline  silicon  films  tested  in  Reference  [1]  are  not  the  same  as 
those  tested  in  Reference  [2].  In  Reference  [1]  the  deposition  temperature  of  the  undoped 
films  was  higher,  and  resulted  in  a  larger  average  grain  size.  Furthermore,  in  Reference  [1] 
the  doping  was  performed  during  the  crystal  growth,  while  in  Reference  [2]  it  was  done 
afterwards. 

1.  REFERENCE  [1] 

Summary 

In  this  work  the  authors  studied  the  fracture  characteristics  of  relatively  thick  (up 
to  7.5  /an)  poly- Si  films  fabricated  using  procedures  identical  to  those  used  for  current 
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MEMS  devices.  They  developed  a  fracture  specimen  whose  characteristic  dimensions, 
notch  tip  radius,  and  uncracked  ligament  length  are  comparable  to  the  grain  size  of  the 
polysilicon,  and  present  data  on  the  critical  value  of  the  nominal  7-integral  for  fracture 
initiation,  Jc.  Notwithstanding  the  inherent  difficulties  of  microfracture  mechanics 
experimentation,  the  values  determined  for  Jc  appear  to  be  four  times  higher  than  those 
obtained  using  macro  specimens  of  single  crystal  silicon. 

Fabrication  and  characterization  of  microfracture  specimens 

A  typical  microfracture  specimen  is  shown  in  Figure  la.  The  specimen  is 
essentially  a  deeply  cracked,  double  cantilever  beam,  suspended  from  an  anchor  through 
two  thin  beams.  The  authors’  initial  goal  was  to  load  the  specimen  electrostatically 
through  the  comb  actuators.  In  fact,  the  generated  force  for  this  design  was  not  sufficient 
to  initiate  crack  propagation,  and  the  fabricated  specimens  were  loaded  mechanically 
through  a  10  /mi  diameter  probe  inserted  along  the  crack  surfaces,  as  will  be  explained 
subsequently. 

The  fracture  devices  were  fabricated  by  standard  surface  micromachining 
techniques,  using  a  single-mask  process.  First,  wafers  underwent  wet  oxidation  at  1 100°C 
for  14  h,  to  grow  a  2.5  /an  thick  SiC>2  release  layer,  which  is  visible  in  (the  bottom  of) 
Figures  lb  and  lc.  This  was  followed  by  the  formation  of  the  polysilicon  films  by  CVD. 
The  undoped  films  were  grown  at  610°C  at  a  total  pressure  of  244  mTorr,  using  Silfr  gas; 
in  all  cases,  a  2.5  or  5  /an  thick  film  was  formed  in  one  run,  and  for  thicker  films,  a  second 
deposition  was  performed  after  exposure  to  the  laboratory  ambient.  The  B-doped  films 
were  grown  at  1100°C,  using  a  combination  of  SiCfeH  and  B2H6  gases,  following  a  brief 
nucleation  step  at  775°C  in  SiHU;  the  several  thicknesses  studied  were  each  grown  in  a 
single  deposition  run,  of  course  using  different  deposition  times  for  the  varying 
thicknesses.  Because  the  surfaces  of  the  as-deposited  poly-Si  films  were  very  rough,  some 
of  the  films  were  subjected  to  a  chemo-mechanical  polishing  with  a  colloidal  silica  slurry 
(Syton),  which  converted  the  rough  surface  of  the  polysilicon  films  to  a  mirror  finish.  The 
polysilicon  films  were  then  oxidized  to  form  a  masking  layer  for  a  later  polysilicon  etch. 
This  wet  oxidation  was  carried  out  at  1000°C,  for  100  to  300  min,  depending  on  the 
desired  final-masking  oxide  thickness. 

The  micromechanical  structures  were  patterned  using  standard  optical  litography 
techniques.  The  masking  oxide  layer  was  dry  etched  in  CHF3,  C2F6,  and  He,  flowing  at  20, 
8,  and  45  seem,  respectively,  at  2.4  Torr  and  350  W  rf,  and  an  etch  rate  of  300  nm/min. 
The  polysilicon  etch  recipe  contains  an  initial  native  oxide  removal  step,  to  remove  any 
Si02  on  the  surface.  This  etching  was  done  for  40  s  at  400  mTorr  and  200W,  with  SF6, 
Cl2,  and  He  flow  rates  of  50,  80,  and  120  seem,  respectively.  The  main  portion  of  the  etch 
was  performed  at  400  mTorr  and  200  W,  with  Cl2  and  He  flow  rates  of  80  and  120  seem, 
respectively,  using  an  optical  endpoint  detector.  The  etch  rate  was  400  nm/min,  and  a  1 
min  overetch  was  included.  Both  etch  steps  were  done  at  an  anode-to-wafer  gap  spacing 
of  6  mm.  The  selectivity  was  24:1  with  respect  to  SiQj.  This  etching  sequence  is  quite 
successful  in  achieving  smooth,  nearly  vertical  sidewalls  in  the  thick  polysilicon  films. 

The  final  step  was  to  release  the  structures  in  hydrofluoric  acid,  which  removed  the 
release  oxide,  as  well  as  any  remnant  masking  oxide.  This  etching  step  was  done  for  10 
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min,  so  that  some  oxide  remains  beneath  the  large  features  of  the  structures,  leaving  them 
anchored  to  the  substrate,  while  the  smaller  features  are  fully  released  and  free  to  move. 

In  both  doped  and  undoped  films,  the  individual  grains  have  an  average  diameter 
of  about  0.5  fm 1  (Reference  [1]  reports  1.0  fim),  as  seen  in  the  cross-sectional  TEM 
images  and  the  planar  view  micrographs  shown  in  Figures  lb  and  lc;  it  is  clear  that  the 
microstructure  of  the  fabricated  devices  is  nearly  columnar. 

Figure  lb  shows  the  microstructure  of  a  two-run  film  which  exhibits  a  second 
nucleation  layer  of  randomly  oriented  grains  which  formed  at  the  beginning  of  the  second 
run,  most  likely  due  to  a  thin  native  oxide  on  the  surface  of  the  first  layer.  The  most 
striking  feature  of  the  microstructure  is  the  density  of  vertically  oriented  {111}  twin 
boundaries,  which  are  known  to  be  the  major  intragranular  defect  in  CVD-grown 
polysilicon  (T.I.  Kamins  and  T.R.  Cass,  Thin  Solid  Films,  Vol.  16,  pp.  147,  1973).  The 
columnar  grains  form  by  “evolutionary  selection”  (A.  van  der  Drift,  Philips  Res.  Repts., 
Vol  22,  pp.  267, 1967)  and  result  in  an  overall  <1 10>  texture,  which  was  confirmed  by  x- 
ray  pole  figures.  The  boron-doped  films  (Figure  lc)  were  deposited  at  a  higher 
temperature  and  exhibit  a  lower  twin  boundary  density,  but  approximately  the  same 
average  grain  size  as  the  undoped  films  (cf.  Figures  lb  and  lc).  The  lower  twin  boundary 
density  in  the  doped  films  is  thought  to  be  due  to  the  higher  growth  temperature. 

Dislocations  were  also  present  in  both  the  doped  and  undoped  films,  but  the  high 
density  of  twin  boundaries  prevented  their  systematic  study.  These  dislocations  may  be 
associated  with  the  residual  compressive  stresses  present  in  as-deposited  and  oxidized 
films,  which  can  be  as  high  as  200  MPa  (H.  Kahn,  S.  Stemmer,  K.  Nandakumar,  A.H. 
Heuer,  R.L.  Mullen,  R.  Ballarini  and  M.A.  Huff,  “Mechanical  properties  of  relatively 
thick,  surface  micromachined  polysilicon  films,  ”  1996  Proceedings  of  the  9fh  IEEE 
International  Workshop  on  Micro  Electro  Mechanical  Systems,  MEMS  96,  pp.  343, 
San  Diego,  CA,  11-15  February  1996,  New  York:  IEEE). 


Figure  1-  Microstructure  of  polysilicon  microfracture  specimens:  (a)  shows  a  typical 
specimen;  (b)  shows,  for  an  undoped  CVD  film,  a  planar  view  (top  picture) 
micrograph  of  the  grain  and  twin  boundaries,  and  a  cross-sectional  TEM  (bottom 
picture);  (c)  shows,  for  a  boron  doped  CVD  film,  a  planar  view  (top  picture) 
micrograph  of  the  grain  structure,  and  a  cross-sectional  TEM  (bottom  picture) 

(Reference  [1]). 
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Figure  2-  (a)  Micrograph  of  wedge  loaded  microfracture  specimen;  (b)  schematic 
drawing  of  microfracture  specimen;  (c)  variables  and  nominal  dimensions  used  in 
the  mathematical  model  (Reference  [1]). 
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3  um 


Figure  3-  SEM  micrograph  of  a  fractured  6  /im  ligament  B-doped  specimen.  The 
crack  path  is  almost  entirely  transgranular  (Reference  [1]). 


Test  procedure 

A  schematic  drawing  of  the  loaded  fracture  specimen  is  shown  in  Figure  2b,  while 
the  nominal  dimensions  are  shown  in  Figure  2c.  The  thickness  t  of  the  specimen  and  the 
length  x  of  the  uncracked  ligament  were  either  2.5,  5,  7.2  and  7.5  jum,  and  6,  10,  and  20 
// m,  respectively.  The  comb  free-travel  distance  a  was  either  8  or  18  /im.  To  cause  failure 
of  these  microdevices,  a  rigid  probe  was  used  to  wedge  open  the  2  //m  wide  notch,  as 
shown  in  Figure  2a.  The  shaft  of  the  tungsten  carbide  probe  is  0.020”  in  diameter  and 
tapers  to  a  conical  point,  with  an  inclusive  angle  of  15°,  to  a  tip  radius  of  5  fim.  As  the 
probe  was  moved  toward  the  crack  tip,  its  position  at  crack  initiation  was  observed  and 
recorded  on  video.  This  information  in  turn  was  used  together  with  a  finite  element 
method  analysis  to  determine  Jc  ,  a  fracture  parameter  appropriate  for  these  specimens 
because  the  radius  of  the  notch  (~1  fim)  is  comparable  to  the  average  diameter  of  the 
grains  and  other  characteristic  dimensions.  A  fractured  sample  is  shown  in  Figure  3,  and 
demonstrates  that  fracture  is  largely  transgranular  (this  was  typical  of  all  specimens).  The 
finite  element  calculation  procedure  for  the  J-integral,  which  assumed  a  Young’s  modulus 
equal  to  160  GPa  and  a  Poisson’s  ration  equal  to  0.22,  is  not  presented  here  for  brevity. 
The  details  can  be  recovered  in  Reference  [1]. 

Results 
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As  shown  in  Table  1,  the  authors  present  the  critical  J-integral  as  a  measure  of  the 
fracture  resistance.  They  stated  that  the  tests  on  the  2.5  pm  thick  specimens  were 
associated  with  significant  out-of-plane  deflections.  Therefore  the  resulting  unphysically 
high  critical  energy  release  rate  values  were  discarded.  The  remaining  values  range  from 
16  to  62  N/m,  exhibit  a  large  scatter,  and  are  approximately  four  times  higher  than  those 
of  single  crystal  silicon. 


TABLE  I.  Critical  energy  release  rate  Qc )  data. 


Polysilicon 

Jc  (N/m) 

Thickness  (/im) 

Doped 

Syton  polished 

Uncracked  ligament  length  (^m) 

Average 

Standard  deviation 

Number  of  tests 

2.5 

Yes 

Yes 

6 

134 

48 

5 

2.5 

Yes 

Yes 

10 

126 

27 

6 

2.5 

Yes 

Yes 

20 

747 

412 

6 

5.0 

Yes 

Yes 

6 

53 

22 

7 

5.0 

Yes 

Yes 

10 

62 

20 

8 

5.0 

Yes 

Yes 

20 

58 

8 

6 

7.5 

Yes 

Yes 

6 

38 

16 

5 

7.5 

Yes 

Yes 

10 

19 

4 

6 

7.5 

Yes 

Yes 

20 

16 

8 

5 

5.0 

No 

No 

6 

26 

4 

4 

5.0 

No 

No 

10 

24 

8 

4 

5.0 

No 

No 

20 

28 

14 

6 

7.2 

No 

Yes 

6 

25 

11 

3 

7.2 

No 

Yes 

10 

24 

4 

3 

7.2 

No 

Yes 

20 

26 

6 

4 

Table  1-  Critical  energy  release  rate  data  (Reference  [1]). 


2.  REFERENCE  [2] 

Summary 

The  specimens  described  in  this  paper  are  fully  integrated  with  simultaneously 
fabricated  electrostatic  actuators  which  are  capable  of  providing  sufficient  force  to  ensure 
catastrophic  crack  propagation  from  notches  produced  using  micromachining.  Thus,  the 
entire  fracture  experiment  takes  place  on-chip,  without  the  need  of  an  external  loading 
device.  This  is  what  makes  this  work  unique  from  the  others  reviewed  in  this  report. 
Fracture  has  been  initiated  using  both  monotonic  and  cyclic  resonance  loading.  A 
reduction  in  the  nominal  toughness  under  cyclic  loading  is  attributed  to  subcritical  growth 
of  a  sharp  crack  from  the  micromachined  notches.  Fatigue  fracture  has  also  been  observed 
in  specimens  subjected  to  as  many  as  109  cycles.  Preliminary  results  suggest  that 
subcritical  crack  growth  can  result  from  the  applied  loading  (without  the  environment 
attacking  the  crack  front). 

The  critical  displacement  of  the  electrostatic  actuator  at  the  instant  of  fracture  is 
recorded,  and  subsequent  finite  element  modeling  is  used  to  determine  the  stress  intensity 
at  the  notch  tip.  Since  the  notch  tip  is  quite  blunt  (the  tip  radius  is  1  pm),  the  maximum 
principal  stress  and  the  critical  value  of  the  nominal  7-integral  for  fracture  initiation,  Jc,  are 
reported. 
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The  actuation  mechanism  exploited  to  create  the  necessary  far  field  force  for 
fracture  was  electrostatic  attraction  in  interdigitated  comb  structures.  This  mechanism  is 
quite  robust  and  is  the  basis  of  the  widely-used  accelerometers  that  comprise  most  current 
automotive  air  bag  sensors.  As  illustrated  in  Figure  4a,  when  a  voltage  is  applied  between 
two  sets  of  fingers,  the  capacitance,  C,  of  the  structure  will  increase  as  the  movable  fingers 
move  a  distance  x  and  increase  the  finger-finger  overlap.  The  force,  F,  generated  by  such 
electrostatic  actuation  is  given  by 


1  dC  , 
F  =  -—V7 

2  dx 


h  . 
ne-V 1 
8 


(1) 


where  n  is  the  number  of  fingers,  s  is  the  permittivity  of  the  fluid  between  the  fingers,  h  is 
the  height  of  the  fingers,  g  is  the  gap  spacing  between  the  fingers,  and  V  is  the  applied 
voltage  (W.C.  Tang,  T.-C.H.  Nguyen  and  R.T.  Howe,  “Laterally  driven  polysilicon 
resonant  microstructures,  ”  Sensors  and  Actuators,  Vol.  20,  pp.  25-32,  1989).  Such 
electrostatic  actuation  usually  involves  only  modest  numbers  of  interdigitated  combs. 
However,  400  pN  of  actuation  force  has  been  realized  in  bulk  micromachined  single 
crystal  silicon  devices  which  contained  2000  pairs  of  interdigitated  comb  fingers  (M.T.A. 
Saif  and  N.C.  MacDonald,  “A  millinewton  microloading  device ,”  Sensors  and 
Actuators,  A  52,  pp.  65-75,  1996).  The  authors  have  now  generated  similarly  large 
electrostatic  forces  in  surface  micromachined  polysilicon  devices. 

Experiment 

Device 

The  device  used  in  this  investigation  is  shown  schematically  in  Figure  4b.  The 
outer  frame  of  the  actuator  is  fully  released  from  the  substrate  in  the  final  stage  of  MEMS 
processing  and  is  free  to  move  when  subjected  to  the  electrostatic  attraction  of  the  comb 
drives.  The  frame  is  anchored  to  the  substrate  at  the  ends  of  eight  long  thin  beams,  which 
constrain  movement  in  the  vertical  direction,  but  allow  side-to-side  movement.  Two 
different  devices  were  fabricated,  containing  either  1456  or  2040  pairs  of  comb  fingers 
(rather  than  the  18  shown  in  Figure  4b),  which  pull  the  actuator  to  the  left  (in  the 
orientation  of  Figure  4b)  when  a  voltage  is  applied.  This  in  turn  pulls  on  the  connecting 
beam  and  the  top  of  the  fracture  mechanics  microdevice,  which  is  anchored  to  the 
substrate  on  the  opposite  side.  This  creates  a  stress  concentration  at  the  notch  of  the 
specimen,  which 
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fixed  fingers 


movable  fingers 


Actuator 


(b) 

Figure  4-  (a)  Sketch  of  two  sets  of  interdigitated  comb  fingers;  when  a  voltage  is 
applied  across  this  comb  drive,  the  movable  fingers  will  move  to  the  left;  (b) 
schematic  drawing  of  the  fracture  device  which  consists  of  an  actuator  (in  this 
sketch,  the  actuator  has  18  pairs  of  comb  fingers  acting  in  parallel)  and  a  fracture 

mechanics  specimen  (Reference  [2]). 
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Figure  5-  Cross-sectional  drawings  depicting  the  device  fabrication  sequence:  i)  a 
thick  Si02  film  is  grown  on  a  silicon  substrate,  followed  by  LPCVD  deposition  of  a 
polysilicon  film  and  a  Si02  film;  ii)  using  photolithography,  the  device  is  patterned, 
and  the  top  Si02  and  polysilicon  films  are  sequentially  etched;  iii)  the  top  Si02  film 
is  completely  removed,  and  the  polysilicon  film  is  boron-doped;  iv)  the  device  is 
time-released  in  hydrofluoric  acid,  such  that  some  Si02  remains  beneath  the  large 
anchor  pads,  while  the  smaller  features  are  fully  released  and  free  to  move 

(Reference  [2]). 
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Figure  6-  Transmission  electron  micrograph  showing  the  microstructure  of  the 

polysilicon  (Reference  [2]). 

mimics  a  conventional  single  edge-notch  fracture  mechanics  specimen.  When  the 
displacement  reaches  a  critical  value,  the  mechanical  forces  become  sufficiently  high  to 
initiate  fracture. 

Fabrication 

The  test  specimens  were  fabricated  using  standard  surface  micromachining 
techniques  and  a  one-mask  process,  illustrated  in  Figure  5  and  summarized  as  follows. 
(100)  silicon  wafer  substrates  are  oxidized  to  a  thickness  of  2.5  pm,  the  oxide  serving  as  a 
release  layer.  A  5.2  pm  thick  film  of  polysilicon  is  deposited  by  low  pressure  chemical 
vapor  deposition  (LPCVD)  at  580°C  using  SiR*  gas.  The  as-deposited  film  is  at  least 
partially  amorphous  (G.  Harbeke,  L.  Krausbauer ,  E.F.  Steigmeier,  A.E.  Widmer,  H.F. 
Kappert  and  G.  Neugebauer,  J.  Electrochem  Soc.,  Vol.  131,  pp.  675,  1984),  and  is 
annealed  at  1100°C,  during  which  crystallization  to  a  fine-grained  polycrystalline 
microstructure  occurs;  the  microstructure  is  shown  in  the  transmission  electron 
micrograph  in  Figure  6.  A  0.5  pm  film  of  Si02  is  then  deposited  by  LPCVD  at  450°C 
using  SiH4  and  O2  gases  for  use  as  a  masking  oxide.  Standard  photolithography  is  then 
used  to  pattern  the  devices  in  photoresist.  Subsequently,  the  masking  oxide  is  dry  etched 
in  a  CHF3/C2F6  plasma,  followed  by  dry  etching  of  the  polysilicon  in  a  Cl2  plasma.  The 
remaining  masking  oxide  is  chemically  removed  using  HF,  and  the  polysilicon  is  boron- 
doped  in  a  diffusion  furnace  at  1000°C.  (Such  boron-doped  polysilicon  specimens  have  a 
resistivity  of  3xl0'3  ohm-cm.)  The  devices  are  then  released  by  immersion  in  aqueous 
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hydrofluoric  acid  (HF)  to  dissolve  some  of  the  Si02  for  a  specific  time,  such  that  some 
oxide  remains  beneath  the  anchor  pads  leaving  them  attached  to  the  substrate,  while  the 
movable  portions  of  the  device  are  fully  released.  Following  HF  release,  the  devices  are 
placed  in  a  25%  NH4F  solution  before  rinsing,  to  help  prevent  stiction  of  the  moving  parts 
to  the  substrate.  Stiction  refers  to  the  immobilization  of  a  micromachined  specimen  during 
release.  The  NH4F  treatment  was  suggested  by  Houston  et  al.  ( M.R.  Houston,  R. 
Maboudian  and  R.T.  Howe,  Proceedings  of  the  8th  IEEE  International  Conference  on 
Solid-State  Sensors  and  Actuators  and  Eurosensors  IX,  Transducers  95,  pp.  210, 
Stockholm,  25-29  June  1995,  New  York:  IEEE,  also  appears  in  J.  Appl.  Phys.,  Vol.  81, 
No.  8,  pp.  3474-3483, 1997). 


Figure  7-  (a)  SEM  image  of  the  entire  fracture  device  (viewed  at  a  45°  angle);  (b)  a 
close-up  of  the  fracture  mechanics  specimen;  (c)  a  top  view  of  the  fracture 
mechanics  specimen  (rotated  90°)  showing  the  dimensions  (in  pm);  d)  one  of  the 
device’s  long  support  beams,  buckled  as-fabricated  due  to  compressive  residual 
stresses;  e)  a  high-magnification  image  showing  several  of  the  comb  fingers;  the 
fixed  and  movable  combs  are  not  co-planar  (Reference  [2]). 

A  scanning  electron  micrograph  (SEM)  of  a  fabricated  device,  along  with 
magnified  views  of  the  fracture  mechanics  specimen  are  shown  in  Figures  7a-c.  (The 
dimensions  of  the  particular  fracture  mechanics  specimen  shown  in  Figure  7c  are  also 
listed  in  Table  2  as  fracture  mechanics  specimen  A.)  The  1456  pairs  of  comb  fingers, 
situated  in  14  rows,  are  2  pm  wide,  20  pm  long,  and  are  too  small  to  be  distinguished  in 
the  SEM  image.  The  "release  holes"  visible  in  the  upper  portion  of  the  fracture  specimen 
and  the  outer  frame  of  the  actuator  allow  the  HF  etchant  to  flow  beneath  these  areas 
during  release.  The  boron-doped  polysilicon  contains  a  residual  compressive  stress  of 
32  ±5  MPa.  When  the  device  is  released,  these  residual  stresses  cause  the  anchored 
support  beams  to  buckle  slightly,  as  seen  in  Figure  7d.  In  addition,  the  movable  portion  of 
the  actuators  are  deflected  out-of-plane  by  these  residual  stresses,  as  seen  in  Figure  7e. 
Both  of  these  geometric  changes  in  the  device  structure  disappear  during  testing,  as 
described  below.  The  residual  stresses  were  measured  using  passive  micro-strain  gauges 
which  were  simultaneously  fabricated  on  the  same  wafers  (Kahn,  H.,  Stemmer,  S., 
Nandakumar,  K.,  Heuer,  A.H.,  Mullen,  R.L.,  Ballarini,  R.  &  Huff,  M.A.  1996 
Mechanical  properties  of  thick,  surface  micromachined  polysilicon  films.  In  Proc. 
IEEE  Int.  Workshop  on  Micro  Electro  Mechanical  Systems,  MEMS  96,  11-15  Feb 
1996,  San  Diego,  CA,  pp.  343-348.  New  York:  IEEE). 
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Table  2-  Fracture  mechanics  specimen  dimensions.  Specimen  A  is  shown  in 
Figure  7c.  All  dimensions  are  in  microns  (Reference  [2]). 


fracture 

mechanics 

specimen 

doped 

tip 

radius 

uncracked 

ligament 

length 

reduced  beam 

release  end 

width 

width 

A 

K9f 

2 

1 

4.2 

9.4 

15 

75 

156 

B 

no 

2 

1 

4.5 

10 

15 

75 

156 

C 

no 

2 

1 

3 

7 

15 

40 

100 

Experimental  procedure 

The  fracture  test  is  conducted  by  applying  a  dc  voltage  between  the  fixed  and 
movable  sets  of  comb  fingers.  This  was  done  on  a  probe  station,  with  one  probe  in  contact 
with  the  central  fixed  portion  of  the  device,  and  a  second  probe  in  contact  with  one  of  the 
outer  anchor  pads.  A  third  probe  was  positioned  in  contact  with  the  silicon  substrate.  This 
last  probe  was  shorted  to  the  probe  on  the  anchor  pad,  and  both  were  grounded.  At  the 
start  of  each  test,  the  applied  voltage  had  to  be  increased  to  a  certain  level  (~  40  V)  before 
the  movable  portion  of  the  actuator  would  be  pulled  flat  (to  counteract  the  processing- 
induced  residual  stresses),  so  that  the  movable  fingers  were  nearly  in-plane  with  the  fixed 
fingers  (to  within  approximately  1  yum,  the  resolution  of  the  focusing  stage  of  the  optical 
microscope  attached  to  the  probe  station).  After  that  point,  as  the  voltage  was  increased, 
the  actuator  (and  thus  the  released  side  of  the  fracture  specimen)  moved  by  increasing 
amounts  until  a  crack  emanated  from  the  micromachined  notch  and  fracture  ensued.  At 
this  time,  the  fingers  immediately  moved  into  contact  with  each  other,  shorting  the  device 
and  causing  some  melting  of  the  polysilicon.  During  the  test,  the  displacement  of  the 
actuator  was  recorded  visually,  with  an  accuracy  of  ±0.3  tun. 

Devices  were  also  fabricated  from  poly-Si  films  which  were  not  boron-doped,  but 
were  otherwise  processed  identically.  For  these  devices,  the  geometry  of  the  fracture 
specimens  (fracture  mechanics  specimen  B  in  Table  2)  was  essentially  identical  to  the 
doped  devices,  except  that  the  reduced  beam  width  was  10  tun  and  the  uncracked 
ligament  length  was  4.5  pm.  To  achieve  sufficient  conductivity  to  allow  electrostatic 
actuation,  these  devices  were  sputter-coated  with  ~10  nm  of  palladium  following  release. 
The  residual  stresses  in  these  devices  were  smaller,  19±5  MPa,  and  tensile.  This  increased 
the  stiffness  of  the  actuator  support  beams  to  such  an  extent  that  sufficient  displacement  of 
the  actuator  to  cause  fracture  could  not  be  attained.  (The  force  limitation  in  these 
specimens  arose  from  the  actuation  voltage  causing  the  outer  frame  of  the  actuator  to 
bend,  moving  the  comb  fingers  sideways,  and  shorting  the  device.)  However,  because  of 
the  tensile  residual  stress,  the  support  beams  were  not  buckled,  and  the  application  of  an 
ac  signal,  at  the  proper  frequency  (~20  kHz),  caused  the  device  to  resonate.  While  in 
resonance,  a  sufficiently  high  amplitude  displacement  could  be  achieved  with  a  relatively 
low  voltage  to  cause  catastrophic  crack  propagation  and  fracture  of  the  specimen.  To 
investigate  the  effects  of  the  laboratory  ambient,  in  particular  the  relative  humidity,  on 
fatigue  behavior,  some  resonance  tests  were  performed  in  a  vacuum  chamber  which  was 
evacuated  to  a  pressure  of  8  Pa. 
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Though  these  devices  could  not  be  fractured  with  a  dc  signal,  monotonic  fracture 
could  still  be  accomplished  by  mechanically  pushing  the  actuator  with  a  probe,  as  in 
Reference  [1].  The  probe  was  carefully  manipulated  to  ensure  strictly  lateral  stresses; 
though  for  this  experiment,  the  possibility  of  perpendicular  and  out-of-plane  stresses  could 
not  be  eliminated.  For  completeness,  boron-doped  devices  were  also  tested  with  the 
mechanical  probe  technique,  and  the  results  included  here. 

Modeling 

After  determining  the  displacement  at  which  fracture  occurred,  the  device  behavior 
under  electrostatic  loading  was  thoroughly  modeled  using  finite  element  analysis  (FEA)  in 
order  to  determine  the  relationship  between  the  observed  displacement  at  fracture  and  the 
maximum  tensile  stresses,  omax>  as  well  as  the  critical  energy  release  rates,  Jc.  The  finite 
element  model  included  both  the  fracture  mechanics  specimen  and  the  actuator.  The 
geometry  of  the  fracture  specimen  was  determined  from  measured  values  while  the 
dimensions  of  the  actuator  were  taken  from  nominal  design  values.  The  model  was 
constructed  from  approximately  16,000  eight-node  plane  stress  quadrilateral  elements. 
The  mesh  density  was  determined  by  a  study  of  the  convergence  of  the  strain  energy  with 
number  of  elements.  For  all  analyses,  Young’s  modulus  and  Poisson’s  ratio  of  polysilicon 
were  assumed  to  be  168  GPa  and  0.22,  respectively,  and  a  large  strain  formulation  was 
used.  Experimentally  measured  residual  stresses  were  included.  Figures  8  and  9  present 
typical  results  of  the  finite  element  calculations.  Figure  8  shows  the  magnified 
displacements  of  both  the  actuator  and  specimen  for  a  600  pN  load  on  a  boron-doped 
device.  The  effect  of  residual  stress  on  the  upper  right  support  beam  can  be  seen  in  the 
displaced  shape.  In  addition,  the  moment  generated  in  the  fracture  specimen  results  in  a 
small  rotation  of  the  actuator.  Figure  9a  shows  the  discretization  of  the  specimen 
geometry  and  the  contours  of  displacement  magnitude.  Contours  of  maximum  tensile 
stress  in  the  region  surrounding  the  notch  tip  are  shown  in  Figure  9b.  The  residual  stress 
(the  stiffness  of  the  support  beams)  has  a  significant  impact  on  the  actuator  behavior; 
however,  omax  in  the  fracture  specimens  was  independent  of  residual  stress  for  actuator 
displacements  resulting  in  fracture.  For  example,  at  a  7  pm  displacement  of  the  end  of  the 
fracture  specimen,  omax  was  calculated  to  be  4.07  MPa  in  a  linear  analysis  without 
residual  stress  and  4.09  MPa  using  large  strain  and  residual  stresses.  Reported  stress  and 
energy  release  values  presented  are  calculated  from  nonlinear  finite  element  results.  Linear 
calibrations  can  be  used  with  less  than  a  5  percent  difference  from  the  non-linear 
calculation.  For  the  boron-doped  devices,  the  linear  calibrations  indicated  that  Jc=1.27  8 
and  amax=0.58  5,  where  8  is  the  displacement  of  the  end  of  the  fracture  mechanics 
specimen  in  pm,  Jc  is  in  N/m,  and  amax  is  in  GPa.  For  the  undoped  devices,  Jc=1.38  8 
and  omax=0.61  8.  Since  displacement  was  experimentally  measured  at  a  specific  point  on 
the  actuator,  FEA  modeling  was  also  used  to  correlate  that  measurement  to  the 
displacement  of  the  specimen. 

The  FEA  modeling  also  determined  that  the  electrostatic  actuator  generated 
90  percent  of  its  theoretical  force,  as  calculated  by  equation  (1).  Most  likely,  the  modest 
loss  in  force  is  due  to  the  fixed  and  movable  fingers  not  being  exactly  co-planar  during  the 
experiment  and  to  the  slight  rotation  of  the  actuator. 


Figure  8-  Finite  element  calculations  showing  displacements  of  a  boron-doped 
device  under  a  600  pN  electrostatic  load.  The  sketch  greatly  exaggerates  the 
displacements;  the  original  position  of  the  device  is  outlined  in  blue  (Reference  [2]). 
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Figure  9-  Finite  element  calculations  showing  (a)  the  discretization  of  the  specimen 
geometry  and  the  contours  of  displacement  magnitude  and  (b)  the  contours  of 
maximum  tensile  stress  in  the  region  surrounding  the  notch  tip  (Reference  [2]). 
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Experimental  results  and  discussion 
Monotonic  and  cyclic  fracture 

The  deflection  of  a  typical  sample  during  monotonic  dc  loading  is  shown  in 
Figure  10.  The  specimen  fractured  at  a  displacement  of  8.3  pm.  The  movement  of  the 
actuator  increased  linearly  with  the  square  of  the  applied  voltage,  as  expected  from  the 
electrostatic  force  equation  (1). 


Applied  Voltage  Squared 

Figure  10-  The  displacement  of  the  fracture  device  versus  the  square  of  the  applied 
voltage;  the  device  fractured  at  8.3  fim  (Reference  [2]). 

The  displacements  at  the  point  of  fracture,  along  with  the  calculated  Jc  and  omax 
values,  for  the  five  sets  of  experiments  are  shown  in  Table  3.  The  results  for  the  two 
different  techniques  of  monotonically  loading  the  boron-doped  samples  (dc  electrostatic 
actuation  and  pushing  with  a  probe)  are  very  consistent.  This  suggests  that  the  two 
techniques  are  equivalent 
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Table  3-  Experimental  polysilicon  fracture  data  (Reference  [2]). 


monotonically 
(dc  electrostatically) 
loaded 
(specimen  A) 


average 


standard  deviation 


B-doped  polysilicon 

monotonically 
(pushed  with  probe) 
loaded 

(specimen  A) 


displacement  at 
fracture  (pm) 


5.0* 


6.3 


7.1 


7.3* 


7.8 


8.3 


7.0 


.2 


5.5* 


6.0* 


6.8* 


6.8* 


6.8* 


7.2* 


7.3 


7.3 


7.4* 


7.9* 


8.0 


8.5 


average 


standard  deviation 


undoped  polysilicon 

monotonically 
(pushed  with  probe) 
loaded 

(specimen  B) 


average 


standard  deviation 


undoped  polysilicon 


4.6 

29 

2.8 

4.6 

29 

2.8 

4.9 

33 

3.0 

5.0 

35 

3.1 

5.3 

39 

3.2 

5.3 

39 

3.2 

cyclically 

loaded 

(at  atmospheric  pressure) 
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(specimen  B) 

6.0 

50 

3.7 

6.6 

60 

4.0 

average 

5.3 

39 

3.2  | 

standard  deviation 

0.7 

11 

0.4 

undoped  polysilicon 

4.4 

26 

2.7 

5.3* 

39 

3.2 

cyclically 

5.3* 

39 

3.2 

loaded 

6.5* 

58 

4.0 

(at  low  pressure) 

7.2 

72 

4.4 

(specimen  B) 

8.7 

105 

5.3  1 

average 

6.2 

63 

4.0 

standard  deviation 

1.6 

27 

0.9  ! 

*These  devices  contained  actuators  with  2040  comb  fingers;  all  other  devices  had 
actuators  with  1456  comb  fingers. 

The  Jc  values  obtained  for  the  monotonically  loaded  boron-doped  and  undoped 
poly-Si  devices  are  comparable  to,  and  consistent  with  previously  reported  values  for 
specimens  with  similar  blunt  notches  (Reference  [1]).  The  small  differences  between  the 
boron-doped  and  undoped  samples  are  within  the  scatter  of  the  data;  therefore,  any  doping 
effects  on  toughness  must  be  small  or  nonexistent.  The  doping  procedure  also  induces 
some  changes  in  the  surface  morphology  of  the  poly-Si  devices.  Because  a  small  amount 
of  oxidation  occurs,  followed  by  the  removal  of  this  oxide  during  the  release  process,  any 
sharp  features  will  be  somewhat  rounded.  In  addition,  during  the  doping,  thermal 
grooving  occurs  where  the  grain  boundaries  intersect  the  free  surface.  Both  of  these 
features  can  be  seen  in  the  micrographs  presented  in  Figure  12. 

For  the  undoped  devices,  the  Jc  values  derived  under  cyclic  stress  are  lower  than 
under  monotonic  loading.  This  suggests  that  during  resonance  a  sharp  fatigue  crack  grows 
subcritically  to  a  critical  dimension.  (During  the  resonance  experiment,  the  specimens  were 
subjected  to  the  critical  displacement  amplitudes  for  approximately  0.1  to  5  seconds 
before  fracture,  corresponding  to  2000  to  100000  cycles.  This  is  referred  to  as  “fast 
fracture.”)  Additional  specimens  were  set  into  resonance  at  displacements  which  were  too 
low  to  initiate  fast  fracture.  The  specimens  were  allowed  to  continue  resonating  until 
fracture  occurred.  The  results,  including  those  from  the  fast  fracture  and  monotonic 
experiments,  are  plotted  in  Figure  11  as  the  maximum  tensile  stress  generated  at  the  notch, 
omax,  versus  the  number  of  cycles  required  for  fracture;  in  essence,  this  is  a  traditional  S- 
N  (stress  to  failure  versus  number  of  cycles  to  failure)  fatigue  curve.  All  of  these  tests 
were  carried  out  in  a  laboratory  ambient  -  approximately  23°C  and  40  percent  relative 
humidity.  The  results  represented  by  circles  were  obtained  using  the  fracture  mechanics 
specimen  described  above.  The  results  represented  by  triangles  were  obtained  using  a 
fracture  mechanics  specimen  with  the  same  shape  but  different  dimensions  -  these 
specimens  are  listed  as  geometry  C  in  Table  2.  The  notch  width  and  tip  radius  were 
identical  in  both  specimens.  Both  specimens  were  modeled  using  FEA  to  determine  amax 
as  a  function  of  displacement,  and  therefore  all  the  results  are  included  in  the  same  plot 
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Figure  11  illustrates  that  a  large  scatter  exists  in  the  data,  possibly  due  to 
differences  in  the  roughness  of  the  inside  notch  surface  caused  by  dry  etching,  and  that 
fatigue  fracture  can  be  achieved  at  sub-critical  stresses  after  as  many  as  109  cycles. 


Figure  11-  S-N  (stress  to  failure  versus  number  of  cycles  to  failure)  curve  for  fatigue 
fracture.  See  text  for  further  details  (Reference  [2]). 

For  cyclic  loading  at  low  pressures,  the  Jc  values  obtained  are  higher  than  for 
cyclic  loading  in  air,  but  are  lower  than  for  monotonic  loading.  This  suggests  that  fatigue 
crack  initiation  and  growth  occurs  during  cyclic  loading  in  both  environments,  but  the 
process  is  faster  in  air,  presumably  due  to  stress  corrosion  in  the  presence  of  humidity.  In 
the  low  pressure  ambient,  stress  corrosion  at  crack  tips  will  be  less  severe,  but  mechanical 
damage  can  still  occur  during  the  tensile-compression  fatigue  cycles.  This  phenomenon 
was  reported  by  Suresh  (S.  Suresh,  Fatigue  of  Materials,  Cambridge  University  Press, 
New  York,  1991)  for  zero  stress-compression  fatigue  cycles. 

Fractography 

Figure  12a  shows  the  top  view  of  a  fractured  specimen.  The  crack  path  is 
significantly  deflected;  this  was  observed  in  most  specimens.  The  fracture  surface  from  a 
B-doped  specimen,  which  fractured  under  monotonic  loading  at  a  deflection  of  7.9  pm,  is 
shown  in  Figures  12b  and  12c,  and  from  an  undoped  specimen,  which  fractured  under 
cyclic  loading  at  a  deflection  of  5.3  pm,  in  Figures  12d  and  12e.  (These  were  the  two  best 
examples  of  “classical”  fracture  features  which  were  observed.)  In  both  samples,  a  semi- 
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circular  “mirror”  is  present,  characteristic  of  brittle  fracture,  and  surrounds  the  crack- 
initiating  flaw.  Beyond  the  relatively  smooth  mirror  zone,  branching  cracks  can  be  seen 
radiating  from  the  crack  initiation  site.  In  brittle  ceramics,  the  radius  of  the  mirror  is 
known  to  be  6  to  10  times  the  radius  of  the  original  flaw  ( J.J.  Mecholsky,  Quantitative 
fractographic  analysis  of  fracture  origins  in  glass,  in  Fractography  of  Glass  (ed.  R.C. 
Bradt  &  R.E.  Tressler),  pp.  37-73.  New  York:  Plenum  Press,  1994).  For  a  semi-circular 
flaw,  the  stress  intensity  factor  is  given  by 

K  =  1.255 (2) 

where  S  is  the  nominal  stress,  and  a  is  the  flaw  radius  {Mecholsky,  1994).  Using  the 
values  of  amax  from  Table  2  for  S,  this  corresponds  to  a  fracture  toughness,  KIc,  of  2.1  to 
2.7  MPa-m1/2  for  boron-doped  polysilicon  and  1.1  to  1.4  MPa-m1/2  for  undoped 
polysilicon,  using  a  values  of  0.13  to  0.22  pm  and  0.07  to  0.12  pm,  respectively.  It  should 
be  emphasized  that  this  analysis  is  valid  only  for  perfectly  semi-circular  flaws  which  lie 
perpendicular  to  the  stress  axis  -  a  fact  that  could  not  be  determined  in  this  experiment. 
Therefore  the  differences  seen  between  the  doped  and  undoped  polysilicon  KJc  values  may 
not  be  significant 

These  values  are  still  higher  than  values  previously  reported  for  single  crystal 
silicon  (~0.9  MPa-m1/2).  The  higher  fracture  toughness  of  polysilicon  compared  to  single 
crystal  silicon  suggests  the  existence  of  a  “microstructural”  toughening  in  our  fine-grained 
polycrystalline  samples.  Current  research  is  now  being  devoted  to  producing 
microfabricated  fracture  mechanics  specimens  containing  atomically  sharp  cracks  to  study 
this  phenomenon  and  other  aspects  of  the  fatigue  of  polysilicon  (effects  of  humidity, 
temperature,  etc.). 
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Figure  12-  a)  SEM  image  of  a  fractured  specimen,  top  view;  b)  the  fracture  surface 
of  a  doped  specimen;  c)  higher  magnification  image  of  the  area  surrounding  the 
crack  initiation  site;  d)  the  fracture  surface  of  an  undoped  specimen;  e)  higher 

magnification  (Reference  [2]). 

3.  CRITIQUE 

The  electrostatically  actuated  fracture  specimen  described  in  this  chapter  is  the 
only  system  available  for  loading  a  thin  poly-Si  to  fracture  without  an  external  loading 
device.  Thus  represents  a  step  towards  the  development  of  on-chip  test  structures  for 
measuring  the  in-situ  mechanical  properties  of  MEMS. 
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CHAPTER  8 


CORNELL  UNIVERSITY  AND  UNIVERSITY  OF  ILLINOIS  AT  URBANA- 

CHAMPAIGN 


This  chapter  presents  a  very  brief  review  of  a  recent  paper  by  M.T.A.  Saif  and 
N.C.  MacDonald  that  describes  their  development  of  a  load  calibration  device,  and  its  use 
as  part  of  an  on-chip  strength  measurement  device  that  integrates  the  (electrostatic) 
actuator,  the  material  sample,  and  the  calibration  mechanism.  Their  system  relies  on  the 
knowledge  of  the  elastic  modulus  of  the  calibration  beam  described  next;  thus  it  was  used 
to  measure  the  strength  of  single  crystal  silicon.  While  it  cannot  be  directly  applied  to 
poly-Si,  it  is  included  in  this  report  because  it  (1)  reinforces  the  scientist’s  philosophy  of 
on-chip  test  structures,  and  (2)  illustrates  the  clever  procedures  that  are  continuously 
being  developed  by  the  MEMS  community.  The  paper  appeared  as 

[1]  M.T.A.  Saif  and  N.C.  MacDonald,  “ Measurement  of  forces  and  spring 
constants  of  microinstruments,”  Review  of  Scientific  Instruments,  Vol.  69,  No.  3,  pp. 
1410-1422,  March  1998. 

A  schematic  and  an  SEM  micrograph  of  the  left  side  and  part  of  the  right  side  of 
the  integrated  microinstrument  are  shown  in  Figure  1.  The  torsion  sample  is  a  single 
crystal  silicon  pillar  (in  the  direction  out  of  the  plane  of  the  figure);  as  shown  in  the 
schematic,  the  torque  is  achieved  through  the  forces  provided  to  the  lever  arms  by  the 
right  and  left  2000  comb  capacitor  actuators.  The  measured  buckled  shape  of  the  500  pm 
long  and  1  pm  wide  calibration  beam  is  used  in  conjunction  with  an  elaborate  structural 
analysis  of  the  instrument  is  performed  to  determine  the  forces.  The  purpose  of  this  short 
review  is  to  describe  the  authors’  concept.  Therefore  details  of  the  fabrication  procedure 
and  mechanical  analysis  (which  makes  up  much  of  the  above  referenced  paper)  are  not 
repeated  here.  In  fact,  the  authors  emphasize  in  their  paper  that  the  fracture  experiments 
are  intended  to  demonstrate  the  application  of  microinstruments  for  material 
characterization,  and  therefore  the  strength  data  should  not  be  considered  as  conclusive. 

Close-ups  of  the  calibration  beam  (and  the  Vernier  scale  used  to  determine  the 
maximum  displacement  of  the  beam)  and  the  torsion  sample  are  shown  in  Figures  2  and  3, 
respectively.  The  dimensions  of  the  two  pillars  tested  were  measured,  using  an  SEM,  as 
approximately  1  pm  x  1  pm  and  1.5  pm  x  1.6  pm. 

The  experiment  and  calibration  are  conducted  under  an  optical  microscope  using  a 
probe  station.  As  the  voltage  is  increased  on  each  instrument,  its  actuator  approaches  the 
calibrating  beam.  In  turn,  the  gap  St  decreases,  and  the  actuator  touches  the  lever  arm  of 
the  torsion  sample.  At  this  point  the  voltage  is  increased  until  the  calibration  beam  just 
buckles.  The  voltage  on  both  actuators  is  then  increased  equally.  As  the  transverse 
deflection  of  the  beam  rapidly  increases  and  is  recorded  (Figure  4),  the  actuators  slowly 
advance  along  the  axial  direction  and  twist  the  sample  (Figure  5).  At  a  critical  applied 
voltage  the  sample  fractures. 
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The  authors  present  a  very  detailed  and  informative  discussion  of  the  major 
assumptions  of  the  calibration  method,  its  limitations,  and  the  design  considerations  of  the 
calibration  beam.  They  also  suggest  that  the  calibration  procedure  is  currently  being 
verified  using  a  direct  force  sensor. 
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Figure  1-  (a)  Schematic  (not  to  scale)  of  microinstrument  and  torsion  sample,  (b) 

SEM  micrograph  (Reference  [1]). 
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Figure  2-  A  calibrating  beam  and  part  of  the  actuator  that  it  calibrates  (Reference 

[1]). 
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Figure  3-  Close-up  of  torsion  pillar  (Reference  [1]). 


v=D(co5(cos)-1  )/2 


Figure  4-  Optical  microscope  close-up  of  buckling  beam.  The  circles  denote  the 
deflection  predicted  by  theory  (Reference  [1]). 


Figure  5-  Twisting  of  a  torsion  sample  recorded  under  an  optical  microscope 

(Reference  [1]). 


91 


PART  II 


REVIEW  OF  THEORETICAL  WORK 


CHAPTER  9  STATISTICAL  METHODS  FOR  PREDICTING  ELASTIC 

CONSTANTS  OF  POLYCRYSTALLINE  FILMS 

CHAPTER  10  THE  EFFECTS  OF  HETEROGENEITY  AND  ANISOTROPY 
ON  CRACK-TIP  PARAMETERS  IN  CRACKED 
POLYCRYSTALLINE  FILMS 
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CHAPTER  9 


STATISTICAL  METHODS  FOR  PREDICTING 
ELASTIC  CONSTANTS  OF  POLYCRYSTALLINE  FILMS 

1.  INTRODUCTION 

Typical  stress  analyses  of  engineered  structures  rely  on  the  assumption  that  the 
constituent  materials  are  homogeneous,  and  in  most  cases,  isotropic.  These 
approximations  have  proved  to  be  appropriate  for  structures  whose  characteristic 
dimensions  are  much  larger  than  those  of  the  material  microstructure  (steel  and  concrete 
buildings,  for  example).  However,  they  may  be  inappropriate  for  MEMS  devices,  whose 
characteristic  dimensions  are  of  the  same  order  as  those  of  the  constituent  material’s 
microstructure.  Poly-Si  is  comprised  of  anisotropic  single  crystals  possessing  “cubic” 
symmetry.  As  shown  schematically  in  Figure  1,  the  {100}  and  {110}  are  the  most 
common  textures  that  result  from  MEMS  processing  techniques.  The  random  direction  of 
the  principal  material  axis  within  each  grain,  represented  in  this  figure  by  the  angle  y/ , 
results  in  an  inhomogeneous  polycrystalline  structure.  Predicting  the  scatter  in  the  nominal 
elastic  constants  of  structures  comprised  of  polycrystalline  silicon  motivated  the 
development  of  the  two  models  presented  in  this  chapter,  which  can  be  easily  modified  to 
analyze  general  anisotropy.  All  calculations  and  results  presented  in  this  chapter 
correspond  to  the  { 100}  texture.  It  is  again  noted  that  the  films  produced  at  CWRU  and 
MCNC  are  {110}.  Thus  the  numerical  results  cannot  be  directly  compared  to  the 
experimental  data. 


{100|  {100}  1100} 


Figure  1-  Schematic  diagram  of  crystal  orientations  in  (a)  {100},  and  (b)  {110}  thin 
film  textures,  showing  in-plane  rotations,  y/  (From  Reference  [3]  below). 
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Both  models  can  be  used  to  predict  the  average  values  and  scatter  of  the  nominal 
elastic  constants  of  a  film  aggregate  of  anisotropic  crystals.  The  first,  proposed  by  the 
author  and  his  coworkers,  is  a  finite  element  based  Monte  Carlo  method  model  described 
in  the  following  references. 

[1]  R.L.  Mullen,  R.  Ballarini,  Y.  Yin  and  A.H.  Heuer,  “ Monte  Carlo 
simulation  of  effective  elastic  constants  of  poly  crystalline  thin  films,”  Acta  mater.,  Vol. 
45,  No.  6,  pp.  2247-2255, 1997. 

[2]  Y.  Yin,  (iMonte  Carlo  simulation  of  effective  elastic  constants  of 
poly  crystalline  thin  films,”  M.S.  thesis,  Department  of  Civil  Engineering,  Case 
Western  Reserve  University,  Cleveland,  Ohio,  1996. 

In  the  above  referenced  work,  the  microstructure  is  approximated  as  a  Poisson- 
Voronoi  tessellation  and  modeled  explicitly  using  a  finite  element  discretization  of  all  the 
grains.  Universal  results  for  the  average  and  variation  of  nominal  plane  strain  Young’s 
modulus  and  Poisson’s  ratio  are  presented  as  functions  of  the  number  of  grains  within  a 
unit  volume,  and  two  parameters  that  quantify  the  level  of  the  crystalline  anisotropy.  The 
predictions  are  compared  with  formulae  that  are  derived  for  plane  strain  Voigt  and  Reuss 
bounds.  For  polycrystalline  silicon  films  comprised  of  a  large  number  of  grains  the  model 
predicts  an  average  plane  strain  Young’s  modulus  equal  to  168GPa,  and  an  average 
Poisson’s  ratio  of  0.22;  these  values  cannot  be  compared  directly  with  experimentally 
observed  values  of  178GPa  and  0.22  (W.  V.  Sharpe,  Jr.,  B.  Yuan,  R.L.  Edwards  and  R. 
Vaidyanathan,  Measurements  of  Young’s  modulus,  Poisson’s  ratio,  and  tensile 
strength  of  polysilicon,”  Proceedings  of  the  Tenth  IEEE  International  Workshop  on 
Microelectromechanical  Systems,  Nagoya,  Japan,  pp.  424-429, 1997).  The  scientist  and 
coworkers  are  currently  calculating  the  results  for  the  {110}  texture. 

This  model  can  be  improved  by  using  a  more  realistic  geometric  modeling  of  the 
microstructure.  There  are  several  ways  of  doing  this,  depending  on  the  objective.  If  one 
wanted  to  model  a  specific  experimental  specimen,  electron  backscattered  diffraction 
micrographs  of  the  specimen  can  be  used  to  determine  the  geometry  and  crystallographic 
directions  of  the  grains.  This  data  in  turn  can  be  translated  into  a  finite  element  mesh.  If 
the  objective  is  to  improve  on  the  statistical  description  of  the  microstructure  that  results 
from  a  specific  processing  procedure,  an  ensemble  of  micrographs  of  actual  films  can  be 
used  to  calibrate  the  Voronoi  model. 

The  model  suffers  from  a  relatively  serious  drawback.  Being  a  Monte  Carlo  based 
technique,  it  is  computationally  intensive,  and  the  results  are  dependent  on  load  and 
displacement  boundary  conditions,  grain  sizes,  and  structural  geometry.  Moreover,  while 
it  provides  average  values  and  standard  deviations  of  desired  parameters,  it  does  not 
provide  autocorrelation  functions. 

The  second  method  was  developed  to  (1)  eliminate  the  computational 
inefficiencies  associated  with  the  Monte  Carlo  method,  and  (2)  provide  an  autocorrelation 
function  for  the  process,  so  that  it  can  be  flexibly  used  in  stochastic  finite  element  codes.  It 
is  summarized  in  the  following  reference. 
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[3]  D.  Mirfendereski,  A.  Der  Kiureghian,  M.  Ferrari  and  G.  Johnson, 
“Probabilistic  characterization  and  response  prediction  of  micro-electro-mechanical 
sytems,”  Report  ( University  of  California,  Berkeley:  Structural  Engineering, 
Mechanics  and  Materials),  UCB/SEMM-96/04, 1996. 

This  procedure  employs  an  equivalent  continuous-parameter  random  field 
characterization  of  the  crystalline  microstructure,  in  which  the  auto-correlation  function 
for  equivalent  elastic  properties  of  the  field  are  generated  by  a  moving  window  approach 
in  fields  of  simulated  crystals.  As  in  the  first  model,  observed  geometries  can  be  used  if 
available. 

Determining  the  average  elastic  and  plastic  properties  of  polycrystalline  materials 
has  been  the  focus  of  much  research  during  the  past  fifty  years.  Before  the  advent  of  the 
computer,  analytical  models  relied  on  simplifying  assumptions  regarding  the  interactions  of 
individual  crystals  as  well  as  their  shape.  As  far  as  the  scientist  knows,  the  first  analytical 
model  for  an  aggregate  of  cubic  crystals  was  proposed  by  Hershey  (A.V.  Hershey,  “The 
elasticity  of  an  isotropic  aggregate  of  anisotropic  cubic  crystals,”  AS  ME  Journal  of 
Applied  mechanics,  Vol.  21,  pp.  236-240,  (1954)).  Christensen  (R.M.  Christensen,  “A 
critical  evaluation  for  a  class  of  micro-mechanics  models,”  Journal  of  the  Mechanics 
and  Physics  of  Solids,  Vol.  38,  pp.  379-404,  (1990))  provides  an  excellent  critique  of 
different  classes  of  analytical  models  for  predicting  the  elastic  and  inelastic  response  of 
composite  materials.  A  practical  limitation  of  these  models  is  that  they  do  not  provide 
estimates  for  the  scatter  of  the  properties  from  their  mean  value,  as  functions  of  the 
number  of  grains  within  a  unit  volume.  This  information  may  be  very  important  in  certain 
applications.  For  example,  specifications  for  disposable  blood  pressure  sensors  require  one 
percent  accuracy  in  pressure  readings  (/.  Bryzek,  R.  Mayer  and  P.  Barth,  “Disposable 
Blood  Pressure  Sensors  with  Digital  On-Chip  Laser  Trimming,”  IEEE  Solid-State 
Sensor  and  Actuator  Workshop,  Hilton  Head  Island,  South  Carolina,  June  6-9,  pp. 
121-122,  (1988)),  which  translates  into  at  most  one  percent  accuracy  in  predicted 
deflections,  and  hence  elastic  moduli.  Currently  this  performance  is  achieved  by  using  bulk 
micromachined  single  crystal  silicon.  Surface  micromachining  provides  more  flexibility  in 
geometric  design,  but  uses  polycrystalline  silicon.  If  this  material  is  to  be  used  in 
applications  which  require  precise  predictions  of  mechanical  response,  then  a  quantitative 
understanding  of  the  scatter  in  its  elastic  moduli  that  results  from  the  inhomogeneity  of  the 
microstructure  is  paramount.  To  estimate  the  dispersion  in  elastic  constants  accurately, 
one  needs  to  perform  analyses  that  model  the  microstructure  explicitly,  as  is  done  in  the 
above  referenced  papers. 

2.  REFERENCES  [1]  AND  [2] 

These  papers  present  a  Monte  Carlo  simulation  that  predicts  the  nominal  elastic  constants 
and  their  coefficient  of  variation  of  films  comprised  of  a  finite  assembly  of  cubic  crystals. 
Poly-Si  is  such  a  structures.  As  discussed  previously,  the  microstructure  is  typically 
(nearly)  columnar.  This  supports  the  two  dimensional  finite  element  analysis  which  is 
presented  in  the  next  section.  The  analysis  relies  on  a  Voronoi  tessellation  of  the 
microstructure,  and  follows  the  approach  taken  by  Kumar  ( S .  Kumar,  “Computer 
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Simulation  of  3D  Material  Microstructure  and  Its  Application  in  the  Determination  of 
Mechanical  Behavior  of  Crystalline  Materials  and  Engineering  Structures ”,  Ph.D. 
Thesis,  Department  of  Engineering  Science  and  Mechanics,  The  Pennsylvania  State 
University,  (1992)),  who  focused  on  comparing,  for  a  large  number  of  grains,  the 
predictions  from  this  model  with  those  derived  from  available  analytical  models. 

Description  of  numerical  model 


With  respect  to  the  global  coordinate  system  x-y-z  shown  in  Figure  2,  the  thin  film 

is  analyzed  under  plane  strain  conditions,  u.  =-^-^  =  0,  where  uz  is  the  displacement  in 

oz 


the  z  direction,  and  ( )  denotes  any  physical  quantity.  For  a  given  number  of  grains  n,  the 
topology  of  the  microstructure  within  the  unit  volume  shown  in  Figure  2  is  approximated, 
for  each  of  m  simulations,  as  a  Voronoi  tessellation.  The  tessellation  is  constructed  by  first 
generating,  using  a  Poisson  process,  a  set  of  random  nucleus  points  rt  ( i=l,n ).  The  set  of 
all  points  which  are  closer  to  rt  than  to  any  other  nucleus  r-  (i*j)  comprise  the  i*h  Voronoi 


cell.  The  assembly  of  the  n  convex,  planar  edged  cells,  defines  the  tessellation.  Each  grain 
is  discretized  with  a  sufficient  number  of  quadratic  displacement  finite  elements.  The 
details  of  the  algorithm  used  to  construct  the  tessellation  and  the  finite  element  method 
solution  procedure  can  be  found  in  Reference  [2]. 


Figure  2-Poisson-Voronoi  tessellation  of  microstructure  for  /i-grains  (n=5-1000) 

(Reference  [1]). 
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Each  grain  is  assumed  to  possess  cubic  symmetry.  With  respect  to  the  local 
coordinate  axes  1-2  shown  in  Figure  2,  which  represent  the  principal  material  directions  of 
a  typical  grain,  the  stress-strain  relations  are  given,  in  matrix  and  component  form,  by 
Hooke's  law  for  a  cubic  material  in  plane  strain 


cn 

C12 

0  ‘ 

£u 

{<?}  =  < 

g22 

=  [C]{£}S 

C\2 

C11 

0  • 

£ 2,  ’ 

0 

0 

2C44_ 

where  (Jtj  and  ei}  ( i ,  j  =  1,2)  are,  respectively,  components  of  the  symmetric  stress  and 

strain  tensors.  System  1-2  is  rotated  through  angle  8  from  the  global  x-axis.  The 
requirement  that  the  strain  energy  be  positive  definite  places  the  following  restrictions  on 
the  elastic  constants 


C11  >  0  cn  ci2  >  0  C44  >  0 


(2) 


The  level  of  anisotropy  is  quantified  by  the  parameters  A  and  v  ,  defined  as 

2c. 


A  =■ 


^44 


—  c12 
V  = - - — 

CU  C12 


(3a) 

(3b) 


The  combination  (cu  -c12)/  2  represents  the  shear  modulus  associated  with  the  shear 

strain  of  an  element  whose  normals  are  rotated  45°  with  respect  to  the  principal 
directions,  while  v"  is  the  apparent  plane  strain  Poisson  ratio  associated  with  uniaxial 
loading  along  the  principal  directions.  For  isotropic  materials  A-l  and  v  =  v ,  where  v  is 
Poisson’s  ratio.  The  randomness  of  the  orientation  of  the  principal  directions  is  introduced 
by  assigning,  to  each  cell,  a  uniformly  distributed  random  angle  6. 

The  thin  film  structure  is  loaded,  as  shown  in  Figure  2,  through  the  following 
boundary  conditions 


Uy  =  A, 

y=l 

(4a) 

Ux  =CTXy 

=  0,  x  =  0,  x  =  1 

(4b) 

uy=G*y 

=  0,  y=0 

(4c) 

These 

are 

associated  with 

nominal  (average) 

strain  components 

e'„  =  A  / 1,  £^=0,  £nxy  =0.  The  finite  element  method  is  used  to  calculate  the  total 
reaction  forces  (per  unit  thickness)  Fx  and  Fy  shown  in  Figure  3,  which  for  the  unit 

length  edges,  correspond  to  nominal  stress  components  and  a”y ,  respectively.  The 
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associated  nominal  values  of  the  Poisson  ratio  v"  and  plane  strain  Young's  modulus 
£"s£"/[l-(v")2]  are  defined  by  invoking  Hooke's  law  for  the  isotropic  continuum 
under  plane  strain  conditions 


II 

Fn 

1-v" 

vn  0 

\on  -  F 

i-/ 

v” 

l-vn  0 

el  =  a  ► 

yy  y 

,<=0j 

~  (l  +  v")(l-2vn) 

0 

0  (l-2v")_ 

yy 

Equations  (5)  provide 

(6a) 

(6b) 


v"  = 


1 


1+£ 

F. 


E  = 


E”  (l  +  vn)(l-2vn)  Fx 
[l-(vn)2]"  vn[l-(v")2]  A 


The  steps  described  above  are  repeated  m  times  for  each  value  of  n. 

Voigt  and  Reuss  Bounds  for  Plane  Strain 

The  plane  strain  elastic  constants  defined  as  Reuss  bounds  are  isotropic  within  the 
plane  and  are  calculated  by  averaging  the  single  crystal  stiffness  matrix  over  all  possible 
orientations  of  the  random  angle  6.  The  average  stiffness  matrix  is  calculated  as 

[C]Rcuss  ~  —  fK[T]T[C][T]d6  (7) 

2  TL  ® 

where  matrix  [T]  and  its  transpose  [T]t,  which  transform  the  stiffness  matrix  from 
principle  axes  to  those  rotated  by  angle  #,  are  defined  by 


[T]m 


cos2  6  sin2#  cos#sin# 

sin2#  cos2#  -cos#sin# 

2cos#sin#  2cos#sin#  cos2#-sin2# 


(8) 


Equating  the  result  of  the  integration  in  (7)  with  the  stiffness  matrix  in  (5)  provides 

„Reuw  _  C11  +3C12  ~  2°44 


rReujj 


4(cu  +c12) 

?  Reuss  _|_  y  Reuss  ^ 


[l_(vR^)2]  2[l-(vRe““)2] 


(Cn  0,2+2044) 


(9a) 


(9b) 
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The  plane  strain  elastic  constants  defined  as  Voigt  bounds  are  isotropic  within  the 
plane  and  are  calculated  by  averaging  the  single  crystal  compliance  matrix  over  all  possible 
orientations  of  the  random  angle  0.  The  average  compliance  matrix  is  calculated  as 


[sr^\yncrme 


(10) 


where  the  superscript  -1  denotes  matrix  inverse.  Equating  the  result  of  the  integration  in 
(10)  with  the  inverse  of  the  stiffness  matrix  in  (5)  provides 


Voigt  _ 


V  5  = 


2s, ,  +  6s,2 


13  44 


2(2 s,,  -2,s12  +$44) 


(11a) 


in  which 


-pVoigt 
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Figure  3-  Scatter  of  silicon’s  nominal  plane  strain  Young’s  modulus  E"  as  a 
function  of  number  of  grains  n  (Reference  [1]). 


Results 

Typical  results  are  presented  in  Figure  3,  which  show  the  scatter  in  silicon’s 
(cu  =  235cI2  =2 cu  =160 GPa,  (A  =  1.65,  v  =  0.28) )  E"  for  m=100  realizations.  It  is 
observed  that  even  for  1000  grains  the  scatter  is  on  the  order  of  a  few  percent.  This 
implies  that  polysilicon  microdevices  whose  characteristic  dimensions  are  less  than  ~  1mm 
cannot  be  used  in  designs  that  rely  on  knowledge  of  the  elastic  modulus  to  within  a  few 
percent. 

The  coefficient  of  variation  ( COV)  of  E"  of  eight  cubic  materials,  defined  as 


COV  = 


(13) 


is  presented  in  Table  1  as  a  function  of  number  of  grains  in  the  film.  In  Equation  (13)  the 
subscript  “avg”  denotes  the  average  value  of  m  simulations. 

Table  2  (m=100)  shows  that  the  plane  strain  Young’s  modulus  and  Poisson’s  ratio 
predicted  by  the  model  for  1000  grains  he  withing  the  Voigt  and  Reuss  bounds.  Figures 

4  and  5  show,  for  v  =0.1,0.2,03,0.4,  the  effects  of  A  on  E"  and  v”  for 
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Table  1-  Coefficient  of  variation  of  nominal  plane  strain  Young’s  modulus  E"  of 
various  cubic  materials  versus  number  of  grains;  c^x  1CT6  Pa  shown  (Reference  [1]). 


Character  Properties 

Voigt 

Reuss 

Model  (1000  Grains, 

100  Simulations) 

Materials 

ell 

cl2 

c44 

Young’s 

Modulus 

Poisson 

Ratio 

Young’s 

Modulus 

Poisson 

Ratio 

Young’s 

Modulus 

Poisson 

Ratio 

Aluminum 

108.24 

62.16 

28.41 

79.0 

0.349 

78.4 

0.351 

78.7 

0.350 

Ammonium 

Bromide 

29.60 

5.90 

0.258 

20.8 

0.294 

Argon 

5.29 

1.35 

1.59 

4.64 

0.232 

4.60 

0.235 

4.62 

0.234 

Calcium 

Fluoride 

162.80 

43,30 

33.40 

128 

0.274 

121 

0.292 

125 

0.284 

Cesium 

2.59 

0.310 

1.29 

1.98 

0.369 

Diamond 

950.00 

390.00 

430.00 

928 

901 

0.247 

913 

0.242 

Gold 

166.00 

40.00 

94.0 

74.3 

0.442 

86.3 

0.432 

Silicon 

168.00 

66.00 

84.00 

172 

0.212 

165 

0.229 

168 

0.221 

Table  2-  Comparison  of  predicted  mean  values  of  plane  strain  modulus  E"  and 
Poisson’s  ratio  v"  with  Voigt  and  Reuss  bounds  (#*=1000)  (Reference  [1]). 

In  these  plots  E"  is  normalized  with  the  "isotropic"  plane  strain  Young's  modulus 

*■“;  isotropic  .  , 

E  obtained  by  equating  the  first  entry  in  the  stiffness  matrix  for  cubic  symmetry,  for 
A=1  (for  fixed  v")  with  that  of  an  isotropic  material,  i.e. 

—isotropic  _  C, ,  (l  +  V  )(l  —  2V  )  ,, 

E  "  <l-v)J  04) 
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It  is  observed  that  the  mean  value  of  the  elastic  constants  is  relatively  insensitive  to 
the  number  of  grains  for  the  practical  range  A<10.  The  predicted  values  lie  between  the 
Voigt  and  Reuss  bounds.  It  is  also  interesting  to  note  that  negative  values  of  the  effective 
Poisson  ratio  are  realized  as  a  result  of  the  inherent  anisotropy  of  the  constituent  grains. 

Figure  6  shows,  for  v*  =  0.1,0.2,0.3,0.4,  the  COV  of  E"  as  functions  of  A  for 
n=l 0,1 00, 1000.  These  plots  clearly  show  that  the  COV  of  the  nominal  Young’s  modulus 
is  more  sensitive  to  v  than  to  A. 


A 

Figure  4-  Mean  value  of  normalized  nominal  plane  strain  Young’s  modulus 
E" / Elso,rop'c  as  functions  of  anisotropy  parameter  A  (Reference  [1]). 


Figure  7  shows  the  standard  deviation  (SD)  of  v"  ,  defined  as 


SD  = 


(15) 


Figures  6  and  7  can  be  used  to  predict  the  inherent  scatter  in  the  nominal  plane 
strain  elastic  constants  of  a  thin  film  comprised  of  a  cubic  material. 
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Mean  Value  of  Nominal  Poisson  Ratio 


Figure  5-  Mean  value  of  nominal  Poisson’s  ratio  v "  as  functions  of  anisotropy 

parameter  A  (Reference  [1]). 
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Figure  6-  Coefficient  of  variation  of  nominal  plane  strain  Young’s  modulus  E "  as 
functions  of  anisotropy  parameter  A  (Reference  [1]). 
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Figure  7-  Standard  deviation  of  nominal  Poisson’s  ratio  v"  as  functions  of 
anisotropy  parameter  A  (Reference  [1]). 
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3.  REFERENCE  [3] 


In  this  report  the  authors  developed  a  random  field  approach  for  characterizing  the 
flexural  and  axial  behavior  of  polycrystalline  beams.  The  authors  also  suggest  a 
methodology  for  extending  the  method  to  two  and  three  dimensional  grain  structures.  An 
equivalent  continuous  parameter  random  field  characterization  of  the  crystalline 
microstructure  is  developed,  whereby  the  auto-correlation  function  for  equivalent  elastic 
properties  of  the  field  are  generated  by  a  moving  window  approach  in  fields  of  simulated 
crystals.  The  advantages  of  this  method  over  the  previously  described  Monte  Carlo 
approach  (which  they  also  used  to  calculate  benchmark  solutions)  are  that  (1)  the 
randomness  in  the  elastic  modulus  at  a  point  along  the  beam  length  is  incorporated  in  one 
equation,  written  in  terms  of  random  variables  and  the  deterministic  spatial  coordinate, 
and  (2)  length  scale  effects  are  included  through  an  autocorrelation  function  and 
correlation  length.  Therefore  the  same  random  field  model  can  be  re-used  within  the 
context  of  conventional  stochastic  finite  element  methods  to  calculate  the  mean  and 
standard  deviation  of  the  response  statistics. 

The  response  of  a  beam  is  essentially  a  one  dimension  problem  that  depends  on  the 
bending  stiffness  EI(x),  where  E  is  Young’s  modulus,  I  is  the  moment  of  inertia,  and  x  is  a 
point  along  the  beam.  In  this  section  only  the  results  for  constant  inertia  are  summarized. 
The  authors  did  calculate  examples  with  random  variations  in  inertia;  these  can  be 
recovered  in  Reference  [3]. 

Calculation  of  random  field  using  the  finite  element  method 

The  original  predictions  of  the  model  relied  on  a  finite  element  stress  analysis  of 
the  discretized  polycrystalline  structure.  The  first  step  consists  of  constructing  Poisson- 
Voronoi  tessellations  (for  a  given  value  of  the  parameter  that  quantifies  the  average 
number  of  grains  in  a  unit  area)  of  an  ensemble  of  polycrystalline  beams,  as  shown  in 
Figure  9.  Next  the  geometry  is  discretized  using  finite  elements,  and  a  crystal  orientation 
distribution  function  is  assigned.  Consider  the  axial  loading  case  of  Figure  9(a),  which 
shows  a  window  containing  three  grains.  For  each  geometry,  this  window  of  width  W  and 
depth  d  at  position  x  is  subjected  to  a  prescribed  stress  oT ,  and  the  average  displacement 
of  the  right  side  of  the  window  relative  to  the  left  is  calculated.  Dividing  this  displacement 
by  the  width  of  the  window  provides  the  average  strain  £avg ,  and  in  turn  the  effective 

(tension)  Young’s  modulus  within  the  window,  ET  =  oT/£avg .  Repeating  these 

calculations  for  a  sufficient  number  of  points  that  span  the  beam  provides  the  continuous 
effective  modulus  Eeff  (x),  which  represents  the  limit  of  ET  as  W  approaches  zero. 

Typical  realizations  are  shown  in  Figure  8.  This  ensemble  of  Eeff  (x)  realizations  is  then 

used  to  find  the  mean,  standard  deviation,  and  autocorrelation  function  of  the  random  field 
E (x) .  Because  this  random  field  is  to  be  used  as  an  input  into  stochastic  finite  element 

codes,  it  is  discretized  into  a  set  of  deterministic  functions  and  a  finite  number  of  random 
variables  using  the  Expansion  Optimal  Linear  Estimation  (EOLE)  method  ( C.C .  Li  and  A. 
Der  Kiureghian,  “ Optimal  Discretization  of  Random  Fields,”  Journal  of 
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Figure  8-  Random  field  model  (Reference  [3]). 
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Engineering  Mechanics,  ASCE,  Vol.  119,  No.  6,  pp.  1136-1154,  1993).  The  flow  chart 
of  the  procedure  is  shown  in  Figure  8. 

The  authors  also  calculated  the  modulus  functions  for  flexural  loading  using  the 
configuration  shown  in  Figure  9(b).  The  details  of  these  slightly  more  involved 
calculations  are  not  presented  here. 

Because  the  finite  element  calculations  require  repeated  remeshing  of  the  windows, 
the  procedure  is  cumbersome  and  computationally  demanding.  The  authors  subsequently 
developed  explicit  equations  for  the  effective  Young’s  modulus  through  an  analytic- 
geometric  approximation  that  eliminated  the  need  of  finite  element  calculations,  as 
described  next 


Analytic-geometric  derivation 

Consider  a  vertical  slice  of  length  <5r  along  the  polycrystalline  beam  shown  in 
Figure  8.  This  slice  is  shown  schematically  in  Figure  10.  As  the  length  approaches  zero, 
this  slice  reduces  to  a  section  of  the  beam,  which  is  composed  of  n  crystals,  each  with  a 
distinct  and  random  orientation  angle  yr,  and  projected  area  Ar  The  different  shadings 
indicate  the  different  crystals.  The  distance  from  the  axis  to  the  top  surface  of  each  crystal 
is  denoted  by  d{ .  The  contribution  of  each  crystal  to  the  effective  Young’s  modulus  in  the 
direction  of  the  beam’s  axis  is  approximated  using  standard  (stiffness)  transformation 
formulas,  and  can  be  written  in  terms  of  the  angle  between  the  principal  direction  and  the 
axis,  and  the  stiffness  coefficients  defined  in  Equation  (1),  as 

E-t  ,Cjj  ,Cj2  *£44)  (16) 

For  uniform  tension,  the  effective  modulus  is  approximated  as  a  (relative  area)  weighted 
average 


(17) 

^total  i=l 

For  flexural  loading,  the  standard  transformed  cross-sectional  area  approach  is 
used.  The  resulting  analytic  expressions  for  the  effective  modulus  are  not  reproduced  here. 
However,  it  can  easily  be  shown  that  for  both  tension  and  flexural  loading,  in  the  limit  as 
the  slice  approaches  zero,  the  following  expression  results 

n 

£tf(*)  =  2>,£(V',)  (18) 

1=1 

where  the  weights  w,  are  functions  of  the  di . 

Using  Equation  (18),  realizations  of  the  random  field  of  the  effective  Young’s 
modulus  along  a  multicrystalline  beam  are  computed.  One  such  realization,  for  flexural 
loading  and  uniformly  distributed  y/ ,  is  shown  in  Figure  1 1.  By  generating  an  ensemble  of 
such  realizations,  statistical  characteristics  of  the  random  field  (mean,  variance,  and  auto¬ 
correlation  function)  were  obtained  using  a  class  of  models  introduced  by  Vanmarcke  (E. 
Vanmarcke,  Random  Fields:  Analysis  and  Synthesis,  The  MIT  Press,  Cambridge, 
Massachusetts,  1983).  The  details  are  not  reproduced  here.  A  comparison  of  beam 
deflection  using  this  approximate  method  with  the  finite  element  based  method  is  shown  in 
Figure  12.  It  is  observed  that  the  approximate  model  does  a  good  job  at  predicting  the 
deflection. 
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Figure  10-  Cross-section  of  vertical  slice  (Reference  [3]) 


Figure  11-  One  realization  of  the  (Young’s  modulus)  random  field  (Reference  [3]). 


Figure  12-  Comparison  of  beam  deflection  using  multicrystalline  finite  element 
calculations  and  using  analytic-geometric  closed  form  solutions. 

4.  CRITIQUE 

The  theoretical  models  presented  in  this  chapter  demonstrate  the  powerful 
computer-based  stochastic  methods  available  to  predict  the  elastic  properties  of 
heterogeneous  materials  and  structures.  While  they  require  lots  of  time  on  the  computer, 
they  are  superior  to  classic  closed-form  solutions  that  only  provide  average  values  and 
wide  bounds;  they  can  be  used  to  generate  as  many  moments  of  elastic  properties  as 
desired. 
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CHAPTER  10 


THE  EFFECTS  OF  HETEROGENEITY  AND  ANISOTROPY  ON  CRACK-TIP 
PARAMETERS  IN  CRACKED  POLYCRYSTALLINE  FILMS 

1.  INTRODUCTION 

This  section  summarizes  two  papers  that  investigate  the  effects  of  heterogeneity  on 
crack-tip  parameters  in  cracked  polycrystalline  films.  Specifically,  the  models  provide 
quantitative  estimates  of  how  the  local  stress  intensity  factors  and  energy  release  rate 
differ  from  those  associated  with  the  assumption  that  the  material  is  homogeneous  and 
isotropic.  This  information  is  needed  to  determine  whether  fracture  mechanics  analyses  of 
microelectromechanical  system  devices  require  discrete  modeling  of  the  material 
microstructure,  or  whether  the  assumption  of  homogeneity  (and  isotropy)  will  suffice.  The 
first  model,  developed  by  the  author  and  coworkers  as  part  of  their  research  on  poly-Si, 
has  been  submitted  for  publication  as 

[1]  R.  Ballarini,  R.L.  Mullen  and  A.H.  Heuer,  “ The  effects  of  heterogeneity 
and  anisotropy  on  the  size  effect  in  cracked  polycrystalline  films,”  submitted  for 
publication. 

The  second  model,  developed  Abdel-Tawab  and  Rodin  in  their  study  of  columnar 
sea-ice,  is  described  in 

[2]  K.  Abdel-Tawab  and  G.J.  Rodin,  “Fracture  size  effects  and 
polycrystalline  inhomogeneity,”  MSS&M  Report  No.  97/12,  Research  Center 
for  Mechanics  of  Solids,  Structures  &  Materials,  The  University  of  Texas  at 
Austin,  July  1997. 


In  linear  elastic  fracture  mechanics  (LEFM)  theory  of  homogeneous  materials,  the 
nominal  stress  on  is  expressed  in  terms  of  the  stress  intensity  factor  Kj  as 

K, 


°n  = 


,>/2; 


(1) 


where  a  is  the  crack  length  and  F  is  a  dimensionless  function  which  depends  on  the 
loading  and  specimen  geometry.  In  this  one  parameter  system,  the  critical  stress  <Tnc  is 
determined,  for  the  desired  loading,  geometry,  and  crack  length,  by  equating  the  stress 
intensity  factor  with  the  material’s  fracture  toughness,  Klc.  The  associated  crack  stress 
analysis  for  the  homogeneous  material  produces  a  macroscopic  (nominal  or  far-field) 
stress  intensity  factor,  which  will  henceforth  be  denoted  K~ .  The  critical  stress  intensity 
factor  criterion  is  equivalent  to  the  condition  that  the  energy  release  rate,  defined  for  self¬ 


similar  extension  as  G“  = 


(*r)2 


reaches  a  critical  value  Glc 


K 


Ic 


E' 


—  2  y ,  where 
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E'=  E  for  plane  stress,  £"=  £7 (1  - v2)  for  plane  strain,  y  is  the  surface  energy,  E  is 
Young’s  modulus,  and  v  is  Poisson’s  ratio.  LEFM  is  based  on  the  assumption  that  the 
“process  zone”  in  the  immediate  vicinity  of  the  crack  tip,  within  which  irreversible 
deformation  occurs  and  linear  elasticity  breaks  down,  is  sufficiently  small  compared  to 
other  relevant  length  scales,  so  that  there  exists  an  annular  region  surrounding  it  within 
which  the  stresses  are  dominated  by  the  elastic  stress  intensity  factor.  Moreover,  the 
assumption  of  homogeneity  implicitly  assumes  that  the  process  zone  extends  from  the 
crack  tip  to  a  distance  much  larger  than  the  crystal  (grain)  size,  and  thus  engulfs  and  in 
turn  is  engulfed  by  a  sufficiently  large  number  of  grains  that  behave  as  a  homogeneous 
medium.  Abdel-Tawab  and  Rodin  (Reference  [2])  suggest  that  a  material’s  tendency  to 
behave  homogeneously  can  be  screened  by  the  grain  brittleness,  GB,  defined  as 

d 

GB= —  (2) 

rp 

Here  d  is  the  grain  size,  and  rp  is  the  size  of  the  process  zone,  which  can  be  estimated  as 


where  au  represents  a  limiting  stress  appropriate  for  a  given  material  system  (yield  stress 
for  metals,  tensile  strength  for  brittle  materials),  and  c  is  a  constant  close  to  unity.  The 
assumption  of  homogeneity  is  appropriate  for  GB»1\  if  GB  is  not  sufficiently  large,  then 
the  irreversible  deformations  at  the  crack  tip  do  not  overwhelm  microstructure  induced 
inhomogeneity,  and  the  details  of  the  microstructure  surrounding  the  process  zone  may 
have  a  significant  effect  on  the  fracture  process.  Indeed,  because  of  the  randomness  of  the 
distribution  and  orientation  of  the  grains,  these  effects  will  be  of  a  statistical  nature.  It  has 
been  established  experimentally  that  the  fracture  toughness  of  most  engineering  materials 
depends  on  specimen  dimensions.  The  amount  of  this  so  called  “size  effect”,  which  results 
from  differences  in  the  relative  size  of  the  process  zone  in  geometrically  similar  but 
different  sized  specimens,  varies  from  material  to  material.  Much  research  has  been 
conducted  by  the  fracture  mechanics  community  in  the  past  twenty  years  aimed  at 
describing  and  predicting  the  size  effect.  This  effort  has  produced  scaling  laws  through 
cohesive  zone  models  (B.N.  Cox  and  D.B.  Marshall,  Acta  Metall.  Mater.  42,  pp.  341- 
363,  1994),  dimensional  analysis  (Z.P.  Bazant,  ASCE  Journal  of  Engineering 
Mechanics  119,  pp.  1828-1844,  1993),  and  fractal  geometry  concepts  (A.  Carpinteri, 
Mechanics  of  Materials  18,  pp.  89-101, 1994).  The  main  advantage  of  these  formulations 
is  that  they  can  be  easily  implemented  into  analysis/design  tools  for  predicting  crack 
propagation  in  engineered  structures.  However,  by  smearing  out  the  details  of  the 
deformation  mechanisms  involved  in  the  fracture  process,  and  by  treating  the  phenomena 
deterministically  and/or  phenomenologically,  they  lose  information  that  can  be  valuable  in 
the  interpretation  of  scatter  in  experimental  data. 

From  the  limited  available  experimental  data,  it  appears  that  polycrystalline  silicon 
does  not  develop  microcracking  in  the  vicinity  of  a  dominant  crack.  Therefore  there  won’t 
be  any  contributions  from  this  phenomena  to  the  size  effect.  It  remains  to  be  determined 
whether  heterogeneity  produces  significant  size  effects. 
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The  work  summarized  in  this  chapter  focuses  on  characterization  of  the  fracture 
mechanics  parameters  of  a  crack  in  a  columnar  aggregate  of  randomly  orientated,  perfectly 
bonded  crystals  under  plane  deformation.  The  following  question  is  addressed:  given  the 
parameters  that  quantify  the  level  of  crystalline  anisotropy,  the  number  of  grains  in  a  unit 
volume,  and  a  far  field  loading  defined  by  KJ ,  what  are  the  average  values  and  standard 
deviations  of  the  microscopic  (local)  stress  intensity  factors  Kl,  and  Kl„  and  energy 

release  rate  Gl  ?  This  information,  which  should  shed  some  light  on  why  the  fracture 
toughness  of  brittle  polycrystalline  structures  is  prone  to  scatter  and  size  effects,  can  be 
obtained  only  by  using  micromechanical  models  that  treat  the  microstructural  details  and 
associated  deformation  mechanisms  explicitly.  As  stated  in  Reference  [2],  the  ultimate 
goal  of  such  models  is  to  develop  size  effect  maps  analogous  to  those  for  deformation 
mechanisms  (H.J.  Frost  and  M.F.  Ashby,  Deformation  Mechanism  Maps,  Pergamon 
Press,  Oxford,  1982),  and  that  characterization  of  the  individual  size  effects  produced  by 
inhomogeneity  is  the  first  step  towards  this  end.  The  next  logical  step  towards  this  end,  is 
to  include  the  effects  of  the  microcracking  (within  the  grains  or  along  grain  boundaries) 
that  often  develops  in  the  vicinity  of  the  crack  tip. 

In  Reference  [2],  which  focused  on  fracture  mechanics  size  effects  in  sea-ice,  an 
approximate  analytical  model  and  a  finite  element  plane  strain  model  calculated  the 
difference  between  KJ  and  K\  at  the  tip  of  a  long  crack  surrounded  by  a  large  number 
of  brittle  hexagonal  crystals  with  perfectly  bonded  interfaces.  This  local  parameter  is 
needed  because,  for  a  heterogeneous  material,  GJ  and  y  in  the  fracture  criterion  should  be 


replaced  by  Gl  and  /,  where  f  is  the  local  value  of  critical  energy  release  rate.  Abdal- 
Taweb  and  Rodin’s  contribution  showed  that,  in  an  aggregate  comprised  of  a  very  large 
number  of  grains,  the  grains  remote  from  the  crack  tip  do  not  induce  significant 
differences  between  the  macroscopic  and  microscopic  stress  intensity  factors.  For 
example,  their  calculations  showed  that  even  for  highly  anisotropic  crystals,  K\  is  only 
twenty  percent  less  than  KJ  ( Kln  was  negligible  for  the  nominal  mode-I  loading).  They 


concluded  that  for  such  materials  the  contribution  of  inhomogeneity  to  the  size  effect  is 
relatively  weak.  Moreover,  their  results  suggest  that  the  crack-tip  parameters  of  materials 
with  large  values  of  GB  are  controlled  by  only  a  few  grains  near  the  crack  tip,  and  provide 
a  rational  explanation  of  why  fracture  toughness  data  exhibits  significant  statistical  effects. 
As  they  pointed  out,  the  results  are  expected  to  change  dramatically  if  microcracking 
within  grains  or  along  grain  boundaries  is  included  in  the  models.  As  will  be  discussed 
later  on  in  this  chapter,  the  results  are  also  expected  to  change  dramatically  if  the  structure 
is  comprised  of  a  very  few  grains. 

One  of  the  issues  that  was  not  discussed  in  Reference  [2]  is  fatigue.  For  brittle 
polycrystalline  materials,  crack  growth  under  cyclic  KJ ,  or  under  constant  KJ  and 
corrosive  environments,  has  been  quantified  with  power  law  functions  of  the  type 

and  — —  =  D{fCJ  )p ,  where  N  and  t  represent  number  of  cycles  and  time, 


dt 


respectively,  and  C,D,m  and  p  are  experimentally  measured  constants.  Because  m  and  p 
are  typically  larger  than,  say  10,  small  heterogeneity  induced  fluctuations  in  local  crack-tip 
parameters  may  suddenly  become  important  in  lifetime  predictions. 
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In  Reference  [1]  the  authors  developed  a  Monte  Carlo  finite  element  based  model 
for  quantifying  the  size  effect  due  to  heterogeneity  and  anisotropy  in  polycrystalline  films. 
The  calculations  predict  the  average  and  standard  deviation  of  the  microscopic  (local) 
stress  intensity  factors  and  energy  release  rate  of  a  crack  in  a  columnar  aggregate  of 
randomly  orientated,  perfectly  bonded,  orthotropic  crystals  (grains)  under  plane 
deformation.  The  boundary  of  the  near-tip  region  is  subjected  to  displacement  boundary 
conditions  associated  with  a  macroscopic  (far  field  or  nominal)  Mode-I  stress  intensity 
factor  and  average  elastic  constants  calculated  for  the  uncracked  film  with  a  large  number 
of  grains.  The  average  and  standard  deviation  of  the  microscopic  stress  intensity  factors 
and  energy  release  rate,  normalized  with  respect  to  the  macroscopic  parameters,  are 
presented  as  functions  of  the  number  of  grains  within  the  near-tip  region,  and  the 
parameters  that  quantify  the  level  of  crystalline  anisotropy.  It  is  shown  that  for  a  given 
level  of  anisotropy,  as  long  as  the  crack  tip  is  surrounded  by  at  least  ten  grains,  then  the 
expected  value  and  standard  deviation  of  the  crack  tip  parameters  are  insensitive  to  the 
number  of  crystals.  For  selected  values  of  crystalline  anisotropy,  the  probability 
distributions  of  Mode-I  stress  intensity  factor  and  stress  ahead  of  the  crack  are  also 
presented.  The  results  suggest  that  the  size  effect  due  to  heterogeneity  and  anisotropy  is 
weak;  crack  initiation  load  and  direction  are  governed  only  by  the  details  of  the  grains  in 
the  immediate  vicinity  of  the  crack  tip. 

This  model  extends  the  calculations  presented  in  Reference  [2]  to  include 
orthotropic  crystals  and  an  arbitrary  number  of  crystals  within  a  unit  volume. 
Furthermore,  the  model  can  be  generalized  to  quantify  the  effects  of  microcracking  in  the 
process  zone  that  often  surrounds  the  crack  tip  in  a  brittle  system.  It  was  developed  to 
quantify  the  effects  of  heterogeneity  and  anisotropy  on  the  scatter  and  potential  size 
effects  in  the  nominal  fracture  toughness  (critical  value  of  K~)  and  fatigue  data  of 
microelectromechanical  systems.  Using  a  first  generation  microfracture  specimen,  they 
reported  that  the  critical  J  integral  (which  provided  a  nominal  fracture  toughness)  of 
notched  polycrystalline  silicon,  which  is  one  of  the  principal  materials  currently  used  in  the 
fabrication  of  MEMS  devices,  exhibited  significant  scatter,  but  was  not  sensitive  to  the 
number  of  crystals  in  the  uncracked  ligament.  The  number  of  crystals  in  the  uncracked 
ligament  was  larger  than  ten. 

Silicon  possesses  cubic  symmetry,  a  K /c  =0.8  MPa  Jin,  =  \2GPa ,  and  a 

crystal  size  d  ~  05fim.  The  resulting  GB  ~  10  suggests  that  heterogeneity  and  anisotropy 
may  play  a  significant  role  in  the  fracture  behavior  of  polycrystalline  silicon  devices  that 
are  comprised  of  a  small  number  of  grains  (connecting  beams  are  typically  a  few  microns 
wide,  and  may  be  comprised  of  less  than  ten  grains).  The  microstructure  of  typical  devices 
is  (nearly)  columnar.  This  supports  the  assumption  of  plane  strain  in  subsequent  finite 
element  analyses. 

2.  REFERENCE  [1] 

Review  of  near-tip  mechanics 

This  section  defines  the  mechanics  problem  and  a  review  of  the  deterministic  dual 
length  scale  mechanics  near  the  tip  of  a  crack  in  a  polycrystalline  material.  The  average 
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properties  of  the  film,  which  enter  into  subsequent  calculations,  were  calculated  using  the 
Monte  Carlo  method  described  in  the  previous  chapter  (Refer  to  Figure  3  of  the  previous 
chapter,  and  define  l  as  the  average  crystal  dimension).  Figure  1  shows  the  same  aggregate 
with  a  crack,  loaded,  as  explained  subsequently,  with  prescribed  displacements  at  the 
boundary  consistent  with  a  macroscopic  stress  intensity  factor  K"  (or  G“ ).  The  problem 
consists  of  calculating  the  average  and  standard  deviation  of  the  microscopic  crack-tip 
parameters  K\ ,  Kl„  and  Gl  as  functions  of  the  level  of  crystalline  anisotropy  and  the 
number  of  grains  in  a  unit  volume.  This  information  is  needed  because  the  crack  will 
extend  at  a  critical  combination  of  the  microscopic  stress  intensity  factors  K\  and  Kl„ , 
which  are  linearly  related  to  K“ ,  but  otherwise  depend  on  the  details  of  the 
microstructure  in  the  immediate  vicinity  of  the  crack  tip.  Alternatively,  the  crack  can  be 
assumed  to  propagate  at  an  energy  release  rate  Gl  (which  is  defined  subsequently)  equal 
to  a  critical  value  yl .  The  propagation  criterion  in  a  single  crystal,  which  is  beyond  the 
scope  of  this  paper,  is  an  important  issue  that  is  receiving  attention  (A.  Azhdari  and  S. 
Nemat-Nasser,  J.  Mech.  Phys.  Solids  44,  pp.  929-951, 1996;  A.  Azhdari  and  S.  Nemat- 
Nasser,  to  appear  in  Mechanics  of  Materials,  1998). 

The  problem  of  predicting  crack  propagation  in  the  polycrystalline  medium  is 
complicated  by  the  fact  that  function  F  in  equation  (1)  now  depends  on  details  such  as  the 
number  of  crystals  in  a  unit  area,  the  relative  position  of  the  crack  tip  within  the  crystal, 
and  the  orientations  of  the  principal  material  directions  within  each  crystal. 

The  crystalline  structure  is  approximated  as  a  Poisson-Voronoi  tessellation, 
consisting  of  n  grains,  using  the  approach  described  in  the  previous  chapter.  However, 
each  grain  is  now  assumed  to  possess  either  cubic  or  orthotropic  symmetry.  This  is 
accomplished  by  relaxing  the  condition  c22  =  cu  in  Equation  (1)  of  the  previous  chapter, 
which  results  in  three  anisotropy  parameters 


cu  c12 


(4) 


(5) 


v  = 


-12 


c„  +C, 


12 


(6) 
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A  Nominal  energy  release  rate 


Figure  1-  Poisson-Voronoi  tessellation  of  a  crack  in  an  aggregate  of  crystals 

(Reference  [1]). 
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Figure  2-  Average  values  of  normalized  plane  strain  Young’s  modulus 
and  Poisson’s  ratio  as  functions  of  anisotropy  parameters;  v~  =  0,0.1, 0.2, 0.3, 0.4 ,0-5 


(Reference  [1]). 


For  isotropic  materials  A=l,  v  =  v ,  and  A,=l.  For  clarity,  the  following  definitions  from 
the  previous  chapter  are  repeated.  The  nominal  Poisson’s  ratio  and  nominal  plane  strain 
Young’s  modulus  are  defined  as  v"  =v  and  En  /  ( l-(v"))2,  respectively.  The 
normalized  plane  strain  modulus,  E,  is  obtained  by  normalizing  the  nominal  plane  strain 
modulus  with  respect  to  the  “isotropic”  plane  strain  modulus  £lsotrop,c  ^  which  is  obtained 
by  equating  the  first  entry  in  the  stiffness  matrix  C,  for  A=1  (for  fixed  v),  with  that  of  an 
isotropic  material,  i.e.  Ebotropic  =  c, ,  (1  -  2v) /(I  —  v2 ) . 

For  an  aggregate  comprised  of  a  very  large  number  of  grains  (an  average  of  n=750 
grains  within  a  unit  area  of  film).  Figures  2  show,  as  functions  of  A  and  v ,  the  average 
values  of  the  nominal  Poisson’s  ratio  and  normalized  plane  strain  Young’s  modulus. 
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0.12 


Average  plane  strain  modulus 


Figure  3-  Distributions  of  average  plane  strain  modulus  (Reference  [1]). 

Figures  3  present,  for  highly  anisotropic  cubic  and  orthotropic  crystals,  the  fraction 
of  occurrences  of  the  average  plane  strain  modulus.  It  is  observed  that  these  are  nearly 
symmetric  about  the  mean,  and  are  associated  with  a  well  defined  peak.  These 
distributions  will  be  revisited  later  on  as  part  of  the  discussion  of  the  distributions  of 
crack-tip  parameters. 
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Figure  4-  A  schematic  representation  of  the  regions  dominated  by  distinct 
asymptotics  (Reference  [1]). 

The  near-tip  mechanics  is  summarized  through  the  local-global  representation  of  the 
cracked  polycrystalline  plane  region  shown  in  Figure  4.  Let  L  represent  the  macroscopic 
length  scale  (crack  length,  uncracked  ligament  length,  or  any  other  relevant  characteristic 
length),  and  recall  that  l  represents  the  average  grain  size.  For  a  given  (deterministic) 
aggregate,  it  is  expected  that  for  large  values  of  LA,  a  near-tip,  small  scale  heterogeneous 
zone  exists,  and  arguments  analogous  to  those  presented  for  small  scale  process  zones  can 
be  used  to  study  the  evolution  of  the  near-tip  stress  fields  (R.  Ballarini,  P.G. 
Charalambides  and  S.  Islam,  International  Journal  of  Fracture  70,  pp.  275-304,1995; 
M.  Jha,  P.G.  Charalambides  and  R.  Ballarini,  Int.  J.  Solids  Structures  34,  pp.  1849- 
1871,  1997;  M.  Jha,  Ph.D.  Dissertation,  Department  of  Mechanical  Engineering, 
University  of  Maryland  at  Baltimore  County,  1997).  For  example,  in  small  scale  yielding, 
the  near-tip  plastic  zones  that  develop  with  increasing  load  disturb  the  singular  elastic 
solution  giving  rise  to  other  near-tip  asymptotics  such  as  the  HRR  fields  in  deformation 
plasticity  and  power  law  hardening  materials  (J.R.  Rice  and  G.F.  Rosengren,  Journal  of 
the  Mechanics  and  Physics  of  Solids  16  pp.  1-12,  1968;  J.W.  Hutchinson,  Journal  of 
the  Mechanics  and  Physics  of  Solids  16  pp.  13-31 ,  1968).  Similarly,  it  is  expected  that 
for  the  polycrystalline  material,  the  heterogeneity  perturbs  the  singular  fields  predicted  by 
the  homogeneized  aysmptotic  solution  giving  rise  to  dual  asymptotics  consistent  with  the 
dual  length  scales  in  the  near-tip  region.  Moreover,  for  large  Ul,  the  mechanical  response 
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of  the  polycrystalline  material  at  distances  sufficiently  away  from  the  physical  crack  tip  can 
be  studied  by  ignoring  the  microstructure  while  employing  the  properties  of  a 
homogeneous  isotropic  medium  instead.  The  latter  properties  can  be  calculated  using  the 
properties  of  the  individual  crystals  through  homogeneization  techniques;  this  paper  uses 
the  results  presented  in  Figures  2.  The  stresses  associated  with  the  response  of  the 
homogeneous  medium  reflect  average  stress  measures  acting  over  the  length  of  an 
appropriate  microstructural  unit  cell.  Under  such  conditions,  the  average  unit  cell  stresses 
and  displacements  in  the  neighborhood  of  the  crack  tip  are  square  root  singular,  and  can 
be  characterized  via  an  isotropic  macroscopic  stress  intensity  factor  KJ .  These  are  given 

by  (T.L.  Anderson ,  Fracture  Mechanics  Fundamentals  and  Applications,  2nd  Edition, 
CRC  Press,  1995) 
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where  the  r-6  polar  coordinate  is  centered  at  the  crack  tip,  G  =  £72(1  + v)  is  the  shear 
modulus,  k  =  3  -  4v  for  plane  strain,  K  =  (3-v)/(l  +  v)  for  plane  stress. 

As  shown  in  Figure  4,  the  outermost  radius  rs  from  the  physical  crack  tip  of  the 
zone  dominated  by  K~  is  expected  to  depend  on  the  macroscopic  characteristic  length  L, 
and  in  turn,  r3  is  some  fraction  of  L.  Thus,  in  a  homogeneized  medium,  region  III  as 
shown  in  Figure  4  is  defined  as  being  dominated  by  the  isotropic  singular  stress  field  given 
by  Equations  (7)  and  (8).  While  this  may  be  a  realistic  approximation  for  material  systems 
with  large  LA  ratios,  it  is  underscored  that  this  approach  may  only  provide  useful  stress 
field  approximations  at  material  points  located  within  region  III  at  distances  from  the 
physical  crack  tip  which  are  much  larger  than  the  microstructural  characteristic  length  /. 
As  a  result  of  the  microstructure,  it  is  clear  that  distinctly  different  asymptotics  will 
dominate  the  mechanics  at  material  points  in  region  I  (which  contains  the  physical  crack 
tip  and  extends  to  a  distance  r, ),  and  that  these  will  determine  the  crack  propagation  load 
and  direction.  The  extent  r,  of  region  I  from  the  crack  tip  depends  on  the  position  of  the 
crack  tip  relative  to  the  adjacent  crystal,  and  the  degree  of  crystalline  anisotropy.  The 
distinct  asymptotics  dominating  regions  I  and  III  discussed  above  are  assumed  to  develop 
in  heterogeneous  systems  with  large  Ul.  In  polycrystalline  systems  with  nearly  isotropic 
crystals,  the  dominance  of  region  I  would  clearly  prevail  over  region  III,  with  r,  extending 
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to  r3 .  For  high  levels  of  anisotropy,  the  asymptotic  fields  in  regions  I  and  III  are  expected 
to  be  matched  through  other  fields  dominating  the  transition  region  II,  extending  from  r{ 
to  r2. 

In  the  immediate  vicinity  of  the  crack  tip  (region  I),  which  is  anisotropic,  the  stress 
and  displacement  components  are  given  by  (G.C.  Sih,  P.C.  Paris  and  G.R.  Irwin, 
International  Journal  of  Fracture  Mechanics  1,  pp.  189-203, 1965 ) 
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where  Re  (Im)  denotes  real  (imaginary)  part 
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the  Cy  are  the  flexibility  coefficients  (the  components  of  the  inverse  of  matrix  C  in 
Equation  (4)),  transformed  to  the  global  axes,  and  /ij  are  the  roots  (with  positive 
imaginary  parts)  of  the  characteristic  equation 

Cut1  —  ^ci6^  +  (2c12  +  )/i  ~  2c26// +  c22  =  0 

(11) 

The  energy  release  rate  for  self-similar  extension  is  given  as 


G  —  2  Krc22 


Im 


^(/X1+/4)  +  ^/ 


+  (12) 


Because  the  crystals  in  a  real  material  are  randomly  orientated,  the  behavior  of  the  stress 
intensity  factors  and  of  the  stresses  in  the  vicinity  of  the  crack  is  more  complicated;  these 
physical  parameters  can  only  be  discussed  in  terms  of  their  probability  distributions. 
Indeed,  the  scatter  in  the  stresses  ahead  of  the  crack  may  be  large  enough  to  render  the 
previous  near-tip  discussion  irrelevant.  In  other  words,  distinct  asymptotics  may  not  be 
observed,  as  will  be  subsequently  shown. 

Near-tip  finite  element  model 

Consider  the  near-tip  plane  strain  region  shown  in  Figure  1.  The  physical 
dimensions  of  this  zone  are  RxR  (with  R=l)  such  that  R  <  r3,  where  r3  signifies  the 
boundary  of  region  III  consistent  with  Figure  4.  In  this  model,  the  Poisson-Voronoi 
tessellation  approximation  of  the  polycrystalline  structure  is  retained.  The  region  under 
consideration  includes  the  crack  tip,  whose  position  within  the  crystal  is  a  random 
variable.  The  crack  surfaces  are  traction-free.  Mode-I  loading  is  applied  by  prescribing, 
along  the  boundary  of  the  region,  the  displacements  in  Equations  (8).  It  should  be 
understood  that  this  model  does  not  incorporate  finite  geometry  characteristics.  In  other 
words,  R  does  not  directly  represent  a  macroscopic  length,  rather  it  provides  an  indirect 
measure  of  a  macroscopic  characteristic  length.  The  only  requirement  on  R  is  that  stated 
earlier,  i.e.  R<r3.  Various  morphologies  corresponding  to  different  LA  are  modeled  by 
changing  the  number  of  crystals,  n,  within  the  modeled  region.  It  is  expected  that  this 
approach  will  capture  the  various  near-tip  asymptotics  and  the  transitional  characteristics 
from  region  I  to  region  III  through  region  II.  The  relation  between  R  and  L  could  be 
quantified  only  by  modeling  specific  finite  geometries. 

The  microscopic  stress  intensity  factors  are  calculated  using  the  displacement 
correlation  technique  (R.S.  Barsoum,  International  Journal  of  Numerical  Methods  in 
Engineering  11,  p.  85,  1977;  A.R.  Ingraffea  and  C.  Manu,  International  Journal  of 
Numerical  Methods  in  Engineering  15,  pp.  1427-1445,  1980),  which  involves  equating 
the  square  root  portion  of  the  finite  element  displacement  interpolation  of  the  crack 
opening  and  crack  sliding  displacements,  with  the  analytical  values  given  by  Equation  (9b). 
These  in  turn  are  used  to  calculate  the  local  energy  release  rate. 
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Local  stress  intensity  factor 

Figure  5-  Fraction  of  occurrences  of  local  stress  intensity  factor  for  highly 
anisotropic  cubic  crystals  (Reference  [1]). 


Results 


The  average  and  standard  deviation  (denoted  by  subscripts  avg  and  sd, 
respectively)  of  K  s  Kl,  /K™  and  G  =  Gl/GJ‘  are  summarized  in  Table  1.  For  m 
simulations  (the  results  presented  in  Table  1  correspond  to  m=200-400 ),  the  standard 


deviation  of  x  is  defined  as 


I*’ 


mx 


avg 


m  —  1 


.  The  results  for  the  isotropic  crystals 


(A=1=I)  show  that  the  errors  in  the  calculated  stress  intensity  factors  and  energy  release 
rates  are  on  the  order  of  a  few  percent. 

For  cubic  crystals  (1=1),  the  effects  of  anisotropy  and  heterogeneity  on  the  crack- 
tip  parameters  are  relatively  small.  Consider  the  results  for  n=500;  for  the  highest  levels  of 
anisotropy  considered  (A=9,  v  =0.4),  the  average  stress  intensity  factor  is  slightly  reduced 
to  ~  0.87  ± 0.06 (±x  represents  ±  one  standard  deviation),  and  the  average  energy 
release  rate  is  reduced  to  *0.85  ±0.12.  Moreover,  the  results  are  not  sensitive  to  n. 
Polycrystalline  silicon  is  defined  by  A=1.65,  v  =0.282;  for  all  values  of  n,  the  average 
stress  intensity  factor  is  ~  0.93 ± 0.02 ,  and  the  average  energy  release  rate  is  *  0.91  +  0.03 . 

The  results  for  orthotropic  crystals  show  that  anisotropy  and  heterogeneity 
produce  significant  effects  for  highly  anisotropic  crystals.  For  example,  if  n=500,  A=9, 
1=10,  and  v  =0.4,  the  stress  intensity  factor  is  =  0.98  ±0.29,  and  the  energy  release  rate 
is  =  0.89  ±  0.30 .  Moreover,  for  highly  anisotropic  crystals,  the  results  are  significantly 
influenced  by  the  number  of  grains;  for  n=10,  the  standard  deviations  of  the  stress 
intensity  factor  and  energy  release  rate,  for  the  previous  example,  increase  to  0.39  and 
0.49,  respectively. 


Local  stress  intensity  factor 


Figure  6-  Fraction  of  occurrences  of  local  stress  intensity  factor 
for  highly  anisotropic,  slightly  orthotropic  crystals  (Reference  [1]). 
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Figures  5-7  show,  respectively,  the  fraction  of  occurrences  of  the  local  Mode-I 
stress  intensity  factor  for  highly  anisotropic  cubic  crystals,  highly  anisotropic  “slightly” 
orthotropic  crystals,  and  high  anisotropic  orthotropic  crystals.  It  is  observed  that  for  cubic 
crystals,  the  distribution  is  relatively  narrow  and  nearly  symmetric  about  the  mean,  while 
for  orthotropic  crystals  the  distribution  is  wider  and  not  symmetric.  Moreover,  for 
orthotropic  crystals,  a  small  population  of  realizations  is  associated  with  stress  intensity 
factors  that  are  more  than  double  of  the  nominal.  Recall  that  the  distributions  of  the 
nominal  elastic  constants  are  symmetric  for  both  cubic  and  orthotropic  crystals. 
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Figure  7-  Fraction  of  occurrences  of  local  stress  intensity  factor 
for  highly  anisotropic  orthotropic  crystals  (Reference  [1]). 
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The  scatter  of  the  normalized  stress  c  = 


V2  nRo^ 

— — —  ahead  of  the  crack  for  highly 
K, 

anisotropic  cubic  crystals,  and  highly  anisotropic  orthotropic  crystals,  are  shown  in 
Figures  8  and  9,  respectively.  In  these  figures  the  distance  from  the  crack  tip,  r,  actually 
represents  the  nondimensional  distance  r/R,  with  R=l.  It  is  observed  that  for  cubic  and 
orthotropic  crystals,  even  at  high  levels  of  anisotropy,  the  dual  length  scale  asymptotics 
disappear  within  the  scatter;  the  stress  is  square  root  singular,  with  upper  and  lower 
bounds  modulated  by  the  deviation  of  the  local  stress  intensity  factor  from  the  mean. 
Furthermore,  the  distribution  of  stress  is  nearly  independent  of  the  number  of  crystals  n. 

It  is  important  to  note  that  the  effects  of  anisotropy  are  expected  to  increase 
drastically  if  (1)  cracking  within  crystals  or  along  grain  boundaries  is  taken  into  account, 
and  (2)  if  the  structure  is  comprised  of  a  very  few  grains.  The  authors  are  currently 
calculating  these  effects. 
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Table  1-  Calculated  crack  tip  parameters  (Reference  [1]). 


135 


3.  REFERENCE  [2] 


This  paper,  as  far  as  the  scientist  knows,  is  the  first  to  show  that  heterogeneity 
produces  a  weak  size  effect  in  polycrystalline  films,  and  instead  induces  statistical  effects 
controlled  by  a  small  number  of  grains  surrounding  the  crack  tip.  The  authors  accomplish 
this  through  several  models,  the  first  analytic,  the  second  finite  element  based. 

The  material  that  was  considered  is  S2  ice,  which  is  formed  by  columnar  grains,  as 
shown  in  Figure  10.  The  global  cartesian  axes  are  defined  such  that  the  x3  axis  is  along 
the  thickness  direction,  and  the  negative  x,  axis  is  along  the  crack  line.  The  grains  are 
assumed  to  be  hexagonal  (transversely  isotropic)  crystals  whose  c-axes  are  randomly 
oriented  in  the  x3  =  0  plane.  Thus  a  sufficiently  large  sheet  is  transversely  isotropic  with 

the  symmetry  axis  along  the  x3  axis.  The  local  coordinates  within  a  grain  are  x3 ,  which 
coincides  with  the  c-axis,  and  the  x, ,  which  coincides  with  the  global  x3  axis. 

Asymptotic  analysis 

The  authors  first  formulate  an  asymptotic  statistical  estimate  of  AK ,  which 
represents  the  stress  intensity  factor  induced  by  grains  of  diameter  d  that  lie  in  the  region 
between  two  quasi-circular  contours  of  radius  Nd  and  N'd  centered  at  the  crack  tip.  The 
estimate  is  valid  for  sufficiently  large  N  and  involves  the  limit  as  N'  approaches  infinity. 
The  material  outside  this  region  is  assumed  homogeneous. 

AK  is  computed  as  follows.  Each  grain  is  replaced  with  an  equivalent 
transformation  strain,  and  the  stress  intensity  factor  it  induces,  8K  is  calculated.  Then, 
using  methods  of  probability  theory,  AK  is  calculated  as  the  sum  over  SK  ’s.  The 
equivalent  transformation  strain  y  (tensor)  of  a  reference  grain  with  stiffness  tensor  C , 
embedded  in  a  material  with  stiffness  tensor  C,  is  determined  from  Eshelby’s  equivalence 
condition  (T.  Mura,  Micromechanics  of  Defects  in  Solids,  Martinus  Nijhoff  Publishers, 
The  Hague,  1982),  which  is  written  in  terms  of  the  ambient  strain  field  £  and  Eshelby’s 
tensor  5,  as 


(13) 


C(£  +  Sy)  =  C(£  +  Sy-y) 


Because  the  stress  field  induced  by  y  at  a  distance  significantly  larger  than  d  is  of  interest, 
the  authors  did  not  calculate  the  exact  solution  of  Equation  (13).  Rather  they  developed 
an  approximate  solution  by  assuming  that  (1)  the  reference  grain  is  circular,  (2)  the  elastic 
interaction  between  grains  is  neglected,  and  (3)  the  variation  of  £  within  the  reference 
grain  is  neglected.  These  assumptions,  together  with  the  assumption  that  the  ambient  field 
is  fully  described  by  a  Mode-I  loading  with  stress  intensity  factor  K,  provide  the  closed 
form  solution 
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(14) 


7  = 


K 


f2nr 


0 C-CS  +  csr  ( c  -  C)C  t(6) 


where  S  is  Eshelby’s  tensor  of  a  long  cylinder,  the  components  of  r  are  the  angular 
functions  in  Equations  (7)  associated  with  a  Mode-I  plane  strain  solution,  and  r  and  9  are 
the  polar  coordinates  of  the  grain  center  (Figure  12(b)).  As  pointed  out  by  the  authors, 
strictly  speaking,  in  Equation  (6)  both  S  and  r  must  be  derived  as  solutions  of  anisotropic 
elasticity.  One  can  prove  that  the  symmetry  of  C  is  such  that  8K  can  be  calculated  using 
solutions  of  isotropic  elasticity.  Therefore  C  is  redefined  as  an  isotropic  stiffness  tensor 
whose  constants  E  and  v  coincide  with  the  corresponding  in-plane  constants  of  the  original 
transversely  isotropic  C  (and  S  and  r  are  redefined  as  the  solutions  corresponding  to  the 
new  C).  The  constants  E  and  v  were  calculated  using  the  self-consistent  approach  outlined 
by  Wang  and  Schapery  (L.  Wang  and  R.A.  Schapery,  “ Prediction  of  elastic  and 
viscoelastic  properties  of  columnar  ice,”  in  J.P.  Dempsey  and  Y.D.S.  Rajapakse  (eds.), 
Ice  Mechanics  1995,  AMD-Vol.  207,  The  American  Society  of  Mechanical  Engineers, 
New  York,  pp.  33-48).  The  expressions  for  all  tensors  can  be  recovered  in  Reference  [2]. 

As  shown  by  Hutchinson  (J.W.  Hutchinson,  “ On  steady  quasi-static  crack 
growth ,”  Report  AFOSR-TR7 4-1042,  1974)  and  Lambropoulos  (J.C.  Lambropoulos, 
“Shear,  shape  and  orientation  effects  in  transformation  toughening,”  International 
Journal  of  Solids  and  Structures,  Vol.  22,  pp.  1083-1106, 1986) 


8K  = 


4n 


16>/2  1 — v' 


■d2r~V2n»y 


(15) 


where  /1  is  given  in  Reference  [2]. 

The  authors  sum  over  all  crystals,  assume  that  A  is  a  uniformly  distributed  random 
variable,  and  employ  principles  of  probability  theory  to  derive  a  simplified  relation  for  the 
probability  that  AK/K  is  less  than  a  certain  number.  The  details  are  not  reproduced  here. 
For  S2  ice,  which  is  slightly  anisotropic,  their  model  predicts  that  for  N=6,  the  probability 
that  AK/K<0.1  is  above  0.99;  for  /V=9,  it  is  above  0.9999.  These  calculations  show  that  in 
this  material  remote  grains  cannot  induce  significant  stress  intensity  factors,  and 
consequently  the  heterogeneity  produces  a  weak  size  effect.  To  check  whether  highly 
anisotropic  crystals  produce  larger  perturbations  of  the  applied  stress  intensity  factor, 
they  calculated  their  model  for  a  ^-fold  reduction  in  one  of  their  elastic  constants.  Table  2 
shows,  for  selected  values  of  if/,  the  number  of  grains  N  needed  to  ensure  that  AK/K<0. 1 
with  a  probability  equal  to  0.99. 


137 


Figure  10-  Fracture  specimen  with  columnar  microstructure  (Reference  [2]). 
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Crack 


Figure  11-  Analysis  of  A K  induced  by  remote  grains:  (a)  specimen  (b)  reference 

grain  (Reference^]). 
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Numerical  estimates 


In  Reference  [2]  the  authors  also  conducted  finite  element  calculations  that 
reinforce  the  notion  that  grains  remote  from  the  crack  tip  do  not  induce  significant  stress 
intensity  factors.  The  polycrystalline  microstructure  was  approximated  and  analyzed  as  a 
circular  cylindrical  grain  centered  at  the  crack  tip,  surrounded  by  quadrilateral  grains 
(Figure  12).  In  each  numerical  experiment,  the  c-axis  of  the  grain  that  contains  the  crack 
tip  was  fixed,  while  the  c-axes  of  the  remaining  grains  were  either  aligned  with  it,  or 
randomly  oriented. 

The  results  are  presented  in  Figure  13.  The  solid  curve  corresponds  to  S2  ice, 
while  the  dashed  and  dotted  lines  correspond  to  <A=2  and  iff =5,  respectively.  Remarkably, 
each  of  these  curves  represents  four  realizations,  one  of  them  being  the  experiment  for 
which  the  c-axes  of  all  grains  are  aligned,  and  the  other  three  the  experiments  with  random 
directions  of  the  grains  surrounding  the  one  centered  at  the  crack  tip.  These  results 
provide  additional  evidence  that  the  grains  other  than  the  one  that  contains  the  crack  tip 
cannot  significantly  perturb  the  applied  stress  intensity  factor. 
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Figure  12-  Finite  element  mesh:  (a)  grain  that  contains  the  crack  tip  (b)  remaining 

grains  (Reference  [2]). 


K°/K 


N 

l 

6  s 

2 

8 

5 

13 

10 

33 

20 

266 

100 

588 

Table  2-  Evaluation  of  the  size  effect  for  various  levels  of  anisotropy  (Reference  [2]). 


Figure  13-  Normalized  local  stress  intensity  factor  as  a  function  of  the  angle  A 

(Reference  [2]). 
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PART  IV 


SUGGESTIONS  FOR  FUTURE  RESEARCH 
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The  work  described  in  this  report  is  representative  of  the  state-of-the-art  of 
experimental  characterization  and  theoretical  modeling  of  the  mechanical  properties  of 
poly-Si  used  for  MEMS.  It  is  the  opinion  of  the  scientist  that  the  quality  of  the  reviewed 
work  is  excellent.  However,  much  more  needs  to  be  done  if  mechanical  reliability  of 
MEMS  devices  is  to  be  predicted  (and  hopefully  guaranteed). 

The  following  research  is  suggested. 

1.  Performance  criteria  should  be  defined.  At  this  point  of  the  game,  it  is  not 
known  what  potential  loading  and  environmental  conditions  MEMS  devices 
are  expected  to  survive.  Grant  it  these  should  soon  be  understood  as  new 
applications  are  envisioned.  Continuous  communication  between  device 
developers  (who  will  dictate  performance  requirements)  and  mechanical 
characterization  types  is  paramount. 

2.  Robust  tests  should  be  developed  for  measuring  residual  stresses,  elastic 
constants,  tensile  strength,  fracture  toughness,  fatigue  crack  growth  rates,  etc.. 
Different  experimental  techniques  should  be  used  for  quantifying  the  effects  of 
specimen  geometry  and  loading.  Those  techniques  that  are  deemed  the  most 
robust  should  be  considered  as  candidates  for  standards.  To  this  end,  the  use  of 
round-robin  testing  should  be  encouraged.  As  mentioned  in  Chapter  1,  one 
such  exercise  (funded  by  the  Air  Force  Office  of  Scientific  Research)  has 
recently  been  completed.  This  involved  measuring  the  Young’s  modulus  and 
tensile  strength  of  poly-Si  films  fabricated  at  MCNC.  The  participants  included 
W.N.  Sharpe  and  coworkers  at  JHU,  George  Johnson  and  coworkers  at 
BSAC,  S.  Brown  and  coworkers  at  Exponent,  and  W.  Knauss  and  coworkers 
at  Caltech.  The  results  are  presented  in  Table  1  (private  communication  with 
W.N.  Sharpe).  Note  that  the  Caltech  results  are  not  yet  available,  and  that  the 
test  methods  used  to  generate  the  results  have  been  described  in  this  report. 
The  scientist  does  not  have  information  about  the  resonant  bending  tests 
performed  by  the  group  at  Exponent.  The  importance  of  this  program  cannot 
be  overstated.  For  example,  it  resulted  in  a  discussion  at  the  Atlanta  Spring 
1998  ASTM  meeting  in  of  Committee  E08.01  (an  education-focused 
committee  that  serves  at  the  current  home  base  of  MEMS  materials  testing) 
that  launched  a  round  robin  program  for  measuring  residual  stresses  in  films 
produced  at  MCNC.  This  program  is  underway. 

Table  1-  Results  of  round-robin  testing  program  (courtesy  of  W.N.  Sharpe) 


BSAC 

Exponent 

JHU  (using 
second 
generation 
system) 

JHU  (using 
second 
generation 
system) 

Film  thickness 
(microns) 

2.0 

2.0 

•  1.5 

3.5 

E  (GPa) 

174 

137 

136 

142 

Strength  (GPa) 

2.8 

2.7 

1.3 

1.3 
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3.  A  sufficiently  large  data  base  for  mechanical  properties  as  functions  of 
processing  techniques  should  be  developed  using  the  most  robust  procedures. 

4.  Mechanisms  of  fatigue,  stress  corrosion  and  fracture  toughness  have  to  be 
studied  and  understood.  It’s  been  convincingly  demonstrated  that  poly-Si  is 
prone  to  static  fatigue.  Crack  growth  rates  should  be  measured  for  materials 
processed  under  different  conditions  and  operating  under  different 
environments.  Some  evidence  was  presented  in  Chapter  7  suggesting  that  poly- 
Si  may  fatigue  in  the  absence  of  humidity.  This  is  but  one  example  of  an 
unexpected  mechanism;  others  may  appear  in  micron  size  structures. 

5.  Theoretical  models  should  be  developed  for  characterizing  and  predicting 
crack  growth  rates. 

6.  Improved  statistical  models  should  be  developed  for  predicting  the  scatter  in 
elastic  constants  of  poly-Si  and  other  material  systems. 

7.  Statistical  fracture  models  should  be  explored  for  predicting  the  scatter  in 
toughness  of  poly-Si  components  comprised  of  a  few  grains. 

8.  Tasks  1-7  should  be  performed  for  other  materials  that  are  currently  used  to 
fabricate  MEMS,  and  for  the  new  material  systems  that  will  be  developed  for 
applications  that  current  materials  cannot  be  used  (high  temperature,  for 
example). 

9.  Experiments  should  be  designed  to  explore  potential  failure  mechanisms  in 
these  new  material  systems,  including  interfacial  fracture  and  creep  in  the  SiC- 
poly-Si  structural  system  developed  by  M.  Mehregany  and  coworkers  (K.N. 
Vinod,  C.A.  Zorrnan,  A.  Yasseen  and  M.  Mehregany,  “Fabrication  of  low 
defect  density  3C-SiC  on  Si02  structures  using  wafer  bonding  techniques ,” 
Journal  of  Electronic  Materials,  in  press )  at  CWRU  and  shown  in  the  Figure. 
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Figure  1-  XTEM  image  of  SiC  (for  high  temperature  applications)  grown  on  as- 
deposited  poly-Si  (reference  provided  in  text). 


